24/12/17 10:43:34.51 W/3j2qe4.net
# # 行列 mat の各行の共通要素を求める
common_elements_in_matrix <- function(mat) {
# ベクトルの場合
if (is.vector(mat)) {
return(mat)
}
# 行列でない場合はエラーを返す
if (!is.matrix(mat)) {
stop("Input must be a matrix or a vector.")
}
# 行数を取得
num_rows <- nrow(mat)
# 行数に応じて処理を分ける
if (num_rows == 1) {
common_elements <- mat[1, ] # 1行の場合はそのまま返す
} else {
common_elements <- Reduce(intersect, split(mat, row(mat))) # 複数行の場合は共通要素
}
# 共通要素がない場合は NULL を返す
if (length(common_elements) == 0) {
return(NULL)
} else {
return(common_elements)
}
}
310:卵の名無しさん
24/12/17 10:44:39.65 W/3j2qe4.net
d4=t(combn(8,4)) ; tail(d4)
r3=t(combn(70,3)) ; tail(r3)
gr=expand.grid(0:1,0:1,0:1)
f=\(mat) all(apply(gr,1,\(y) length(setdiff(common_elements_in_matrix(mat[y==1,]),mat[y==0,])))<2)
calc=\(x){
rbind(d4[x[1],],d4[x[2],],d4[x[3],]) |> f()
}
r3[12345,] |> calc()
apply(r3,1,calc) |> which() -> idx
re=r3[idx,]
tail(re)
g=\(x) rbind(d4[x[1],],d4[x[2],],d4[x[3],])
nrow(re)
(x=g(re[sample(6720,1),]))
f(x)
matrix(c('赤','青','黄','緑','白','黒','紫','茶')[x],ncol=4) |> noquote()
311:卵の名無しさん
24/12/19 20:34:36.80 cvKGuYd1.net
Wolframによるカウント
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= s4=Subsets[Range[8],{4}];
In[2]:= grp=Subsets[s4,{3}];
In[3]:= Short[grp];
In[4]:= rat=Subsets[Range[3],3];
In[5]:= f[x_]:=(
z1=Length@Complement[Range[8],Flatten@x];
z28=Table[
dead = Intersection@@x[[rat[[i]]]];
alive= Flatten@x[[Complement[Range[3],rat[[i]]]]];
Length@Complement[dead,alive],{i,2,8}];
AllTrue[Flatten@{z1,z28},#==1&]
)
In[6]:= ans=Select[grp,f];
In[7]:= Short[ans]
Out[7]//Short= {{{1, 2, 3, 4}, {1, 2, 5, 6}, {1, 3, 5, 7}}, {<<3>>}, <<6717>>, {{2, 6, 7, 8}, <<2>>}}
In[8]:= Length[ans]
Out[8]= 6720
312:卵の名無しさん
24/12/19 20:34:45.09 cvKGuYd1.net
問題の意味は小学生でもわかる問題（手作業で数えるのは無理だと思う）
８本の注射液があり、赤,青,黄,緑,白,黒,紫,茶のラベルが貼られている。
いずれも十分な量があり、そのうち１本が致死的な毒薬であることがわかっている。
３匹のネズミに注射液を混合して注射してどれが毒薬かを特定したい。
複数の注射液を混合しても毒性は変わらないとする。
各々のネズミに注射できるのは１回だけである。
赤 青 黄 黒 の混合液
赤 白 黒 茶 の混合液
青 緑 白 黒 の混合液
を３匹に１つずつ注射すれば特定できる。
ネズミを区別しないとして、このような組み合わせは何通りあるか。
313:卵の名無しさん
24/12/21 01:57:46.63 GFvlmbSb.net
(* 点Pの座標 *)
pointP = {Latitude -> φ_P, Longitude -> λ_P};
(* 大円を定義する2点 *)
pointA = {Latitude -> φ_A, Longitude -> λ_A};
pointB = {Latitude -> φ_B, Longitude -> λ_B};
(* 大円の距離を計算する関数 *)
GreatCircleDistance[pointP, GeoPath[{"GreatCircle", {pointA, pointB}}]]
314:卵の名無しさん
24/12/21 02:04:29.20 GFvlmbSb.net
(* 点の座標（緯度と経度） *)
point = {30, 45}; (* 例: 緯度30°, 経度45° *)
(* 大円を定義する2点の座標（緯度と経度） *)
linePoint1 = {0, 0};
linePoint2 = {0, 90};
ToCartesian[{lat_, lon_}] := {
Cos[lat Degree] Cos[lon Degree],
Cos[lat Degree] Sin[lon Degree],
Sin[lat Degree]
}
p = ToCartesian[point];
lp1 = ToCartesian[linePoint1];
lp2 = ToCartesian[linePoint2];
greatCircleNormal = Cross[lp1, lp2];
angle = ArcSin[Norm[Cross[greatCircleNormal, p]]/Norm[greatCircleNormal]];
distance = angle (* 球の半径が1の場合 *)
315:卵の名無しさん
24/12/21 08:39:16.68 GFvlmbSb.net
rm(list=ls())
# URLﾘﾝｸ(i.imgur.com)
source('toolmini.R')
Plot(-20,5,axes=F,zero=F)
pt(0i,'◯',col='red',cex=0.75)
A=3+3i
pta(A)
B=-17+3i
pta(B)
C=Conj(B)
pta(C)
D=Conj(A)
pta(D)
Polygon(A,B,C,D)
P=-3+3i
pta(P)
Q=-3-17i
pta(Q)
R=3-17i
pta(R)
Polygon(A,P,Q,R)
θ=33
theta=pi/180*θ
ro=exp(1i*theta)
a=A*ro
pta(a)
b=B*ro
pta(b)
c=C*ro
pta(c)
d=D*ro
pta(d)
Polygon(a,b,c,d,Col='red')
I=intsect(P,Q,c,d)
pta(I)
J=intsect(C,D,a,b)
pta(J)
seg(P,I,col='green',lwd=2)
seg(J,D,col='green',lwd=2)
abs(P-I)
abs(J-D)
316:卵の名無しさん
24/12/21 08:39:33.21 GFvlmbSb.net
rm(list=ls())
source('toolmini.R')
solve=\(L=20,S=6,θ=33,verbose=TRUE){
A=S/2+1i*S/2
B=-(L-S)+1i*S/2
C=Conj(B)
D=Conj(A)
P=-S/2+1i*S/2
Q=-S/2-(L-S)*1i
R=S/2-(L-S)*1i
theta=pi/180*θ
ro=exp(1i*theta)
a=A*ro
b=B*ro
c=C*ro
d=D*ro
I=intsect(P,Q,c,d)
J=intsect(C,D,a,b)
if(verbose){
Plot(-L,S)
pt(0i,'*',col='red')
Polygon(A,B,C,D)
Polygon(A,P,Q,R)
Polygon(a,b,c,d,Col='red')
seg(P,I,col='green',lwd=2)
seg(J,D,col='green',lwd=2)
}
c(縦緑=abs(P-I),横緑=abs(J-D))
}
solve()
317:卵の名無しさん
24/12/21 19:38:34.87 2ZERFniW.net
NSolve[R Cos[x]/4==10 && R (Pi/2-x)/(2 Pi)==6,{R,x},Reals]
318:卵の名無しさん
24/12/23 05:02:42.03 jcrTO2R5.net
t=seq(-pi,pi,le=100)
library(rgl)
mat=cbind(cos(t),cos(t)^2,sin(t))
plot3d(mat)
319:卵の名無しさん
24/12/23 08:22:05.35 Odl9ursf.net
寄生虫トキソプラズマは感染したヒトを魅力的にして、モテやすくする⁈
URLﾘﾝｸ(nazology.kusuguru.co.jp)
人間においては寄生虫の感染が統合失調症に影響するという有意な調査報告があり、
320:卵の名無しさん
24/12/26 03:50:20.36 TtLzjWgm.net
library(RcppAlgos)
a=2:(-10)
b=(-3)^a
e=exp(1)
#abs(d-e)
f=permuteGeneral(0:2,length(a),rep=T)
#tail(f)
fn=\(c) abs(e-sum(b*c)) < 1e-5
ans=f[apply(f,1,fn),]
ans
if(is.matrix(ans)) {apply(ans,1,\(c) sum(b*c))
}else sum(b*ans)
321:卵の名無しさん
24/12/26 04:34:29.84 TtLzjWgm.net
library(RcppAlgos)
a=3:(-11)
b=exp(a)
e=exp(1)
f=permuteGeneral(0:2,length(a),rep=T)
fn=\(c) abs(pi-sum(b*c)) < 1e-5
ans=f[apply(f,1,fn),]
ans
if(is.matrix(ans)) {apply(ans,1,\(c) sum(b*c))
}else sum(b*ans)
322:卵の名無しさん
24/12/26 04:45:36.46 TtLzjWgm.net
(* 初期設定 *)
eValue = N[E, 50]; (* e の高精度値 *)
piValue = N[Pi, 50]; (* π の高精度値 *)
base = piValue;
(* 桁数の設定 *)
digits = 10;
(* 結果を格納するリスト *)
coefficients = {};
(* 変換のメインループ *)
For[i = 1, i <= digits, i++,
integerPart = Floor[eValue];
AppendTo[coefficients, integerPart];
eValue = (eValue - integerPart) * base;
]
(* 結果の表示 *)
coefficients
323:卵の名無しさん
24/12/29 20:50:22.37 wMWbX3gA.net
Worth to mention that Tag Times ... is Protected message is in 90% cases the result of a missing semicolon.
324:卵の名無しさん
24/12/29 20:51:14.03 wMWbX3gA.net
sim[] := Module[(* 1:横綱 2:大関 3:関脇*)
{h12,h23,h13,bout,lastWinner,flg,count,winner},
h12={10,9}; (* 通算成績 1 vs 2 *)
h23={8,7}; (* 通算成績 2 vs 3 *)
h13={6,5}; (* 通算成績 1 vs 3 *)
bout=RandomSample[Range[3],1][[1]];(* 初戦の試合番号=抜け番 *)
lastWinner=0;(*前勝者*)
flg=False;(*連勝フラッグ*)
count=0;(*試合数*)
While[!flg,
count++;
If[bout==1,(* 2 vs 3 *)
i=2-RandomVariate[BinomialDistribution[1,h23[[1]]/Total@h23],1][[1]];(* 勝者のindex *)
h23[[i]]=h23[[i]]++;(*通算成績更新*)
winner={2,3}[[i]];(*勝者番号*)
bout={3,2}[[i]](*次の試合番号*)
];
If[bout==2,(* 1 vs 3 *)
i=2-RandomVariate[BinomialDistribution[1,h13[[1]]/Total@h13],1][[1]];
h13[[i]]=h13[[i]]++;
winner={1,3}[[i]];
bout={3,1}[[i]]
];
If[bout==3,(* 1 vs 2 *)
i=2-RandomVariate[BinomialDistribution[1,h12[[1]]/Total@h12],1][[1]];
h12[[i]]=h12[[i]]++;
winner={1,2}[[i]];
bout={2,1}[[i]]
];
flg = winner==lastWinner;(*連勝フラッグ更新*)
lastWinner=winner(*前勝者更新*)
];
{winner,count}(*優勝者,試合数*)
k=10^5;
res=Table[sim[],k];
Winner=res[[All,1]];
Games=res[[All,2]];
Counts[Winner]/k//N//ReverseSort
Histogram[Games]
Mean@Games//N
]
TableForm[Transpose@Table[sim[],20],TableHeadings->{{"優勝者","試合数"},None}]
325:卵の名無しさん
25/01/18 20:02:47.91 u0+3XfyQ.net
solve=Vectorize(function(z) integrate(function(x) dlnorm(x,0.7,0.5)*dlnorm(x-z,0.9,0.4),0,Inf)$value)
curve(solve,0,5,xlab='days after onset',ylab='false negative rate')
solve(0:4)
plot(solve(1:24/24))
326:卵の名無しさん
25/01/18 20:17:15.74 u0+3XfyQ.net
pdfZ=Vectorize(function(z) integrate(function(x) dlnorm(x,0.7,0.5)*dlnorm(x-z,0.9,0.4),0,Inf)$value) # wind - incu
curve(pdfZ,-5,5,xlab='days after onset',ylab='false negative rate')
solve=Vectorize(function(d) integrate(pdfZ,d,Inf)$value)
solve(0)
solve(1:4)
curve(solve,0,5)
327:卵の名無しさん
25/02/09 06:28:59.93 rE+LhOlF.net
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= gr=Tuples[{"+","-","*","/",""},8];
In[2]:= join[x_] := (
{x1,x2,x3,x4,x5,x6,x7,x8}=x;
"1"<>x1<>"2"<>x2<>"3"<>x3<>"4"<>x4<>"5"<>x5<>"6"<>x6<>"7"<>x7<>"8"<>x8<>"9"
)
In[3]:= op[x_] := ToExpression@join@x
In[4]:= re = op /@ gr;
In[5]:= solve[n_] := join /@ Select[gr,op@#==n&]
In[6]:= solve[1234]
Out[6]= {1+234*5-6+78-9, 1234+5+67-8*9, 1234-5-67+8*9}
In[7]:=
328:卵の名無しさん
25/02/25 16:20:23.53 q+PMcb4A.net
ステップ 2. 「HKEY_LOCAL_MACHINE > SOFTWARE > Microsoft > Windows NT > CurrentVersion > Winlogon」に移動します。
ステップ 3. 「DefaultPassword」と「DefaultUserName」にスクロールダウンし、ダブルクリックします。ウィンドウが表示され、保存されたパスワードが明らかになります。
329:卵の名無しさん
25/02/26 05:25:02.59 RrSkJf5m.net
. 拡張機能やプライバシー設定の影響
広告ブロッカーやプライバシー系のアドオンがログインを妨げることがあります。
対処法:
about:addons で拡張機能を無効化して試す。
about:preferences#privacy で「強化型トラッキング防止機能」を「標準」にする。