宝の埋め方の組み合わせを列挙して分類したら
>95のようになるのは同意？
探検
分からない問題はここに書いてね448
レス数が1000を超えています。これ以上書き込みはできません。
宝の埋め方の組み合わせを列挙して分類したら
>95のようになるのは同意？
この結果面白いね
問題が２つ見つけるまでやって多く獲った方どっち？だったらイーブンだけど、１つ目を先に獲った方が勝ち、とすると差が出る
この場合、２番目を先に見つける確率にもきっと差があるのだろう
納得できないのは直観的に納得できないだけで、そういうことになるよなぁとはわかっていると思います
今ちょっと考えているのが、遅く見つけたほうが勝ちというルールで行うなら
Ｑ：ABCDEFGHIJKL
Ｐ：AEIBFJCGKDHL
では、Ｐ君の方が勝率は高いということ。
じゃあ、Ｑに対してＰ以上に勝率の高い文字列（検索順序）は存在するはずだけど
それらを具体的に求める方法は？
とか考えてしまう。
で、頭がぐーるぐーるるるるるる
先に２つの宝を発見した方を勝者とするのでやってみた。
n=2
ABC
DEF
の場合
> t232=treasure2(2,3,2)
P1st Q1st even
5 4 6
短軸方向探索Pが先に２つの宝を発見する埋め方：５通り
> print(matrix(LETTERS[t232$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5]
[1,] A A B B D
[2,] D E D E E
長軸方向探索Qが先に２つの宝を発見する埋め方：４通り
> print(matrix(LETTERS[t232$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4]
[1,] A A B C
[2,] B C C D
同時に２つめの宝を発見する埋め方：６通り
> print(matrix(LETTERS[t232$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] A B C C D E
[2,] F F E F F F
n=3
ABCD
EFGH
IJKL
の場合
> t342=treasure2(3,4,2)
P1st Q1st even
27 26 13
> #短軸方向探索Pが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B B C C C D E E E E F
[2,] E F I J K E F I J K I J K K F I J K I
[,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27]
[1,] F F G G G I I J
[2,] J K I J K J K K
> #長軸方向探索Qが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B C C C C C D D D D D
[2,] B C D G H C D G H D E F G H E F G H I
[,20] [,21] [,22] [,23] [,24] [,25] [,26]
[1,] E E F F G H H
[2,] G H G H H I J
> #同時に２つめの宝を発見する埋め方
> print(matrix(LETTERS[t342$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,] A B C D D E F G H H I J K
[2,] L L L J L L L L K L L L L
R の コードは http://tpcg.io/RZo4hd
EFGH
を一般化するとこうなるかな？
横一列に並んだ ABCDEFGH の中からランダムに2つを選んで宝を隠しておく
ABCDEFGH の中から宝を探す順番は8！通りある
最初の宝を見つけた時点で終わるものとするとき、
8！通りの探し方の中で最も有利な探し方はどのような探し方か
じゃんけんと一緒。
どんな列取ってきてもその先頭文字を末尾に回した探索には負ける。
A■■■□
E■■■■
I ■□■■
けど、ある特定の列に対して最も勝率が高いのはどれだろうとは気になる。
けど、先頭文字を末尾に回した奴が一番勝率高くなるのかな。
ABCDEFGに対してなら
BCDEFGAが一番勝率高い気がする
二つ目の当たりとの組み合わせは考慮しなくていい
当たりがどの座標のマスに置かれても
要素の個数は変化しないので
どの方向からの探査によっても確率は変化しない
やっと構造がなんとなくわかってきた
自分の脳みその弱さが悲しくなってくる
部屋の数についての帰納法でいけるんじゃね？
主張は
部屋の数が n の時 P:A[1]A[2]…A[n] に引き分けないという条件下で勝つ確率最大なのは Q:A[2]A[3]…A[n]A[1]。
以下Qの探索順をB[i]とする。
n=3では多分成立。
n<k で成立として n=k のとき。
P が Q に勝つのはA[1]とA[2]以外に宝が配置されるときでその確率は (n-2)/C[n,2]。
引き分けるのはA[1]、A[2]に配置されるときで確率1/C[n,2]。
よってQがPに勝つ確率は (C[n] - 1 - (n-2)/(C[n,2] - 1)。
容易にA[1]≠B[1]の場合はコレより確率は大きくならないとわかる。
A[1] = B[1]の場合を考えればよい。
このとき引き分けないという条件下では宝箱はA[1]以外の２つに配置される場合でその場合Qの勝つ条件付き確率の最大値は (C[n-1] - 1 - (n-3)/(C[n-1,2] - 1)。
多分 (C[n] - 1 - (n-2)/(C[n,2] - 1) ≧ (C[n-1] - 1 - (n-3)/(C[n-1,2] - 1)より成立。
おお、そういう風に片付くのか
帰納法で出来ないかとも考えたけど、自分の頭では無理だったのです
これで自分はスッキリしました！
一列じゃなくて長方形型 n×n+1の配置じゃないの？
nの次は (n+1)(n+2)では
縦nマス、横n+1マスのn(n+1)マスのとき
n=1でイーブン
n=2,3で長軸方向探索が有利
n=4以上で短軸方向探索が有利となっているので
数学的帰納法はn=3で適応できないと思う。
考えてる問題が違う。
# ABCDEFGに対してなら
# BCDEFGAが一番勝率高い気がする
library(gtools)
n=7
k=2
perm=permutations(n,n)
Q=perm[1,]
np=nrow(perm)
p1st=numeric(np)
for(i in 1:np){
P=perm[i,]
tre=combn(n,k)
nt=ncol(tre)
re=numeric()
for(j in 1:nt){
re[j]=min(which(tre[1,j]==P),which(tre[2,j]==P))-
min(which(tre[1,j]==Q),which(tre[2,j]==Q))
}
p1st[i]=sum(re<0)
}
plot(p1st)
p1st[which.max(p1st)]
(p.max=which(p1st==15))
print(matrix(LETTERS[perm[p.max,]],ncol=7),quote=F)
#
# ABCDEFGに対してなら
# BCDEFGAが一番勝率高い気がする
一番勝率高い探索順は４通りあった
> print(matrix(LETTERS[perm[p.max,]],ncol=7),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] B C D E F A G
[2,] B C D E F G A
[3,] B C D E G A F
[4,] B C D E G F A
なるほどね
先回り側がEの次の部屋へ進むってことは当たりはFGだからどっちに進んでも同じか
先回り側の順序の最後の2つは決して実行されない
その４つの順序のどれでも最後から3番目のFかGまでで決着が付くから
ABCDEFGに対して一番勝率高い探索順は？
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] A C D E F G B
[2,] B C D E F G A
[3,] C A D E F G B
[4,] C B D E F G A
俺もそっちに興味があるが、証明できる頭脳はない。
整数苦手だからよくわかんない
そうは問屋が卸さないみたいだよ。
縦４マス、横５マスで宝箱を１から７まで増やしてみると
宝が６個になると短軸有意から長軸有意に逆転した。
処理速度の制約であまり大きな数字で検証できないのだが。
> sapply(1:7,function(k) treasure(4,5,k))
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
短軸有利 9 84 463 1776 5076 11249 19797
長軸有利 9 83 453 1753 5075 11353 20057
同等 2 23 224 1316 5353 16158 37666
気長にやってみた。
> sapply(1:20,function(k) treasure0(4,5,k))
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11]
短軸有利 9 84 463 1776 5076 11249 19797 28057 32243 30095 22749
長軸有利 9 83 453 1753 5075 11353 20057 28400 32528 30250 22803
同等 2 23 224 1316 5353 16158 37666 69513 103189 124411 122408
[,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]
短軸有利 13820 6656 2486 695 137 17 1 0 0
長軸有利 13831 6657 2486 695 137 17 1 0 0
同等 98319 64207 33788 14114 4571 1106 188 20 1
４×５の場合
宝：１個 同等
宝：２～５個 短軸有利
宝：６～１３個 長軸有利
宝：１４～２０個 同等
と推移した。
えー
なかなか簡単にはいかせてもらえないね
考えながら仕事片付けるとしよう
それに宝箱１のときは、イーブンだろ。
どこかで派手な勘違いをしたままでやってたぽい
そりゃあっちこっちぐるぐるうううになるわ…
>>131の同じになるってのもとりあえず保留
宝２個でn=30まで計算させてみた。４以上で短軸有利は不変だった。
> t(sapply(1:30,treasure1))
P1st Q1st even
[1,] 0 0 1
[2,] 4 5 6
[3,] 26 27 13
[4,] 84 83 23
[5,] 203 197 35
[6,] 413 398 50
[7,] 751 722 67
[8,] 1259 1210 87
[9,] 1986 1910 109
[10,] 2986 2875 134
[11,] 4320 4165 161
[12,] 6054 5845 191
[13,] 8261 7987 223
[14,] 11019 10668 258
[15,] 14413 13972 295
[16,] 18533 17988 335
[17,] 23476 22812 377
[18,] 29344 28545 422
[19,] 36246 35295 469
[20,] 44296 43175 519
[21,] 53615 52305 571
[22,] 64329 62810 626
[23,] 76571 74822 683
[24,] 90479 88478 743
[25,] 106198 103922 805
[26,] 123878 121303 870
[27,] 143676 140777 937
[28,] 165754 162505 1007
[29,] 190281 186655 1079
[30,] 217431 213400 1154
の解説をして頂けませんか？
答えは2,2 2,3だと思うのですが解答が無くて
よろしくお願いします
お願いします
xは偶数しかありえないのでx=2mとおけば
y≦x-1のとき,
x!+y!=y!(x!/y!)+y!=y!((x!/y!)+1),
(3≦)(x!/y!)+1=x･(x-1)!/y!+1とxは互いに素だから, x!+y!≠x^y.
すなわちx≦y.
3≦xのとき,
x!+y!=x!(1+(y!/x!))は(x-1)(≧2)の倍数.
x-1とxは互いに素であり, x!+y!≠x^y.
すなわちx≦2.
1)x=1のとき, 与式を満足させるyはない.
2)x=2のとき, 2+y!=2^y.
y≧4とすれば,
2+y!=2+24･(y!/4!)>2+3･2^(k-1)>2^k.
すなわちy≦3.
よって求める組は(x,y)=(2,2), (2,3).
できました！
もとの 4x3 型の12部屋で宝箱の数を変えてみると
1と8以上で有利不利無し、それ以外は長軸優先有利となるな
初見での印象よりも随分奥深いなこれ
計算式お願いする
もう少し正確にいうと、
全単射関数 n : N -> N で 数列 a[ i ] を n で並べ替えた数列b[ i ]を
b[ i ] = a[ n(i) ] で定義する。
b[ i ] の収束性、極限値はどうなるでしょう?
部屋が
ABCD
EFGH
IJKL
として
宝物10個のときはABが空きなら縦の勝ち、
AEが空きなら横の勝ち
縦勝ちの宝物9個の配置
CDEFGHIJK
CDEFGHIJL
CDEFGHIKL
CDEFGHJKL
CDEFGIJKL
CDEFHIJKL
CDEGHIJKL
CEFGHIJKL
DEFGHIJKL
横勝ち
BCDFGHIJK
BCDFGHIJL
BCDFGHIKL
BCDFGHJKL
BCDFGIJKL
BCDFHIJKL
BCDGHIJKL
BCFGHIJKL
BDFGHIJKL
以下各個数での勝敗の数
treasures 1: p win 5 q win 5 even 2
treasures 2: p win 26 q win 27 even 13
treasures 3: p win 73 q win 76 even 71
treasures 4: p win 133 q win 140 even 222
treasures 5: p win 167 q win 176 even 449
treasures 6: p win 148 q win 153 even 623
treasures 7: p win 91 q win 92 even 609
treasures 8: p win 37 q win 37 even 421
treasures 9: p win 9 q win 9 even 202
treasures 10: p win 1 q win 1 even 64
treasures 11: p win 0 q win 0 even 12
treasures 12: p win 0 q win 0 even 1
Rでよければこんな感じ
# 宝の数を変化させる
treasure0 <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
return(c(短軸有利=sum(re<0),長軸有利=sum(re>0),同等=sum(re==0)))
}
sapply(1:12,function(k) treasure0(3,4,k))
> sapply(1:12,function(k) treasure0(3,4,k))
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
短軸有利 5 26 73 133 167 148 91 37 9 1 0 0
長軸有利 5 27 76 140 176 153 92 37 9 1 0 0
同等 2 13 71 222 449 623 609 421 202 64 12 1
>128に書いたけど
４x５だと宝箱を増やすと途中で短軸有利から長軸有利に変わっちゃうので自分でもびっくりした。
> 141
よろしければプログラムコードをアップしていただけませんか？Pythonでしょうか？
ちょっと整理してました。
NB. n comb n returns all n length set from 0..m-1
comb =: dyad define
if. x=1 do.
(1,~y)$i.y
elseif. x=y do.
(1,y)$i.y
elseif. do.
((y-1) ,/"0 1 (x-1) comb y-1 ), x comb y-1
end.
)
NB. usage: 3 4 game 2
game =: dyad define
p =. ,/ |: x $ i. */x
q =. i. */x
g =. y comb */x
d =. (<./"1)@(g &((i."1 0)~))
r =.(d p)-(d q)
y, (+/ r<0), (+/ r>0), (+/ r=0)
)
NB. run 3 4 games n for n in 1..12
smoutput 'tre p q even'
smoutput 3 4 game "1 0 (1+i. 12)
ここで実際に動かしてみることができます
https://goo.gl/znRTwf
U[j=1, n]A_j=ℝ となるn個の集合 A_j (*1) について,
j=1,2,...,n で a_j∈A_j となるような変数 a_j を取り,
lim[a_j→α] f(a_j) =k (*2) が全ての j について言えたならば,
lim[x→α] f(x) =k (*2) が言えますか。
例えば, p∈ℚ, q∈ℝ\ℚ とすると, p と q を合わせれば全実数を取ります。このとき,
lim[p→α] f(p) =lim[q→α] f(q) =k
かつ
lim[x→α] f(x) ≠k
となる f(x) は存在しますか。
(*1)αに十分近い要素も含む
(*2)離散的極限
お手数かけました。
残念ながら自分の知識ではアラビア文字のように理解不能でした。
一つ質問ですが
スタート地点Aに宝があるとゲームスタートと同時に
同着でゲーム終了になるけど、ポイントAに宝は設置されるのですか？
R^nからR^nへの写像で、(x_1,...,x_n)を(x_σ(1),...,x_σ(n))にうつすものは連続であることを示して下さい。
示してください
条件収束する級数を考えればa[i]とb[i]の収束性に関係がないことは明らか
その場合は引き分けで終了。
宝の置き方はランダム。
12C2=66通りに等確率で配置。
数字を増やしたらサイトの時間制限を超えて結果がでなくて残念。
尚、>142のRはメモリ不足で停止しました。
NB. usage: 5 6 game 2
NB. run 5 6 games n for n in 1..30
smoutput 'tre p q even'
smoutput 5 6 game "1 0 (1+i. 30)
分散分析でF分布の値の比に F-ratio というのが出てくるの知ってた？
a[i] → c とする。
e>0 とする。
|a[i] - c| ≧ e である i は有限個。
∴ |b[i] - c| ≧ e である i は有限個。
∴ b[i] → c。
離散的極限って離散位相での極限？
だったら Aj が disjoint な集合なら
a[1]→α、a[1]∈A[1]、a[2]→α、a[2]∈A[2] 自体が起こりえないやろ？
誘導位相？
12部屋から6部屋選ぶ組み合わせは924通りしかないのに
20部屋から10部屋だと184756通り、
30部屋から15部屋だと155117520通り、
という感じなのでどうしても時間やメモリを食いますよね
場合分けなどが面倒くさくて疲れ果てたけど、計算結果は>>133と一致。
P1st(n)-Q1st(n) が（偶奇によらず） (n^2-2n-6)(n-1)/6 になったので、n=2,3でＱが、n≧4でＰが有利。
コードはSagemath。
from sage.calculus.calculus import symbolic_sum
,var m,l,k,a,n
P1 = (symbolic_sum((m-1)*(m)-2*l-1, l,1,m-2)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,2*k) + symbolic_sum((m-1-k)*(m-k)+3*k-2*l-1, l,2*k+1,m-2), k,1,a-2)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,m-2), k,a-1,m-2)
).substitute({a:m/2}).substitute({m:n+1})
P2 = (symbolic_sum((m-1)*(m)-2*l-1, l,1,m-2)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,2*k) + symbolic_sum((m-1-k)*(m-k)+3*k-2*l-1, l,2*k+1,m-2), k,1,a-1)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,m-2), k,a,m-2)
).substitute({a:(m-1)/2}).substitute({m:n+1})
def P1st(x):
return P1.substitute({n:x}) if mod(x,2) == 1 else P2.substitute({n:x})
Q1 = (symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,l-1), l,0,a)
+ symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,m-l-1) + symbolic_sum((m-k)^2+k-l,k,m-l,l-1), l,a+1,m-1)
+ symbolic_sum(k^2-2*m*k+m^2+k-m,k,0,m-2)
).substitute({a:(m-1)/2}).substitute({m:n})
Q2 = (symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,l-1), l,0,a-1)
+ symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,m-l-1) + symbolic_sum((m-k)^2+k-l,k,m-l,l-1), l,a,m-1)
+ symbolic_sum(k^2-2*m*k+m^2+k-m,k,0,m-2)
).substitute({a:m/2}).substitute({m:n})
def Q1st(x):
return Q1.substitute({n:x}) if mod(x,2) == 1 else Q2.substitute({n:x})
P1 == 1/24*(6*n^3 + 20*n^2 - n - 27)*(n - 1) # nが奇数のとき
P2 == 1/4*n^4 + 7/12*n^3 - 7/8*n^2 - 13/12*n + 1 # nが偶数のとき
Q1 == 1/24*(6*n^2 + 10*n - 3)*(n + 1)*(n - 1) # nが奇数のとき
Q2 == 1/24*(6*n^2 - 2*n - 5)*(n + 2)*n # nが偶数のとき
>133です。労作ありがとうございます。
コードは全く読めないのですが、宝の数を増やしての計算はこのコードで可能なのでしょうか？
４×５の場合で宝を増やすと
宝：１個 同等
宝：２～５個 短軸有利
宝：６～１３個 長軸有利
宝：１４～２０個 同等
に変化したので差分はどんな関数なのだろうかとか、
５×６ではどうなるのか（メモリ不足で実行できませんでした）とか興味があります。
(n^2-2n-6)(n-1)/6 をグラフ表示してみました。
http://i.imgur.com/Qel6ZAy.png
>>161は多項式にまでするために、部屋をn x (n+1)、宝を２個と特殊化したものです。
#nloc(m,n,k,l)は縦m、横nの部屋で横優先が部屋(k,l)で初めて宝を発見する場合で
#宝が置かれても縦優先に先を越されない部屋の数。
def nloc(m,n,k,l):
q,r = divmod(n*k+l,m)
return (n-q)*(m-k)+q-1-l + ((k-r) if r > k else 0)
#nwin(m,n,c)は部屋が縦m、横n、宝がc個で横優先が勝つ宝の配置の数
def nwin(m,n,c):
return sum(binomial(nloc(m,n,k,l),c-1) for k in range(m) for l in range(n) if k*(n-1)<l*(m-1))
縦優先は縦横を替える。
レスありがとうございます。
コードは読めないのですが、
部屋数から宝部屋の組合せを列挙してどちらが縦横どちらが先にみつけるかを探る手続きで必要な計算式をプログラムが絞りだしてくれるという理解でいいのでしょうか？
いえ、計算式そのものです。数式で書けば
nwin(m,n,c) := Σ[(k,l)∈{0,…,m-1}×{0,…,n-1}, k*(n-1)<l*(m-1)] binomial((n-q)*(m-k)+q-1-l + (k-r)δ(r > k), c-1)、
ただし、n*k+lをmで割った商をq、余りをrとし、δ(P)をPが真なら1、偽なら0である関数とする。
例えば→√2を考えたい時、qの近づけ方は問題ないんでしょうが、pの近付き方を、p_n→√2になるような有理数列p_n上で考えることは出来ないんでしょうか。
→0なんかも、実際に関数に0を入れるわけではなくギリギリまで近付けるように、p自身が√2を取れないのは、定義できないほどの大問題でしょうか。
スタート地点のポイントAに宝があると
ゲーム開始とともに同着でゲーム終了になるので除外する
宝がいくつあったとしても、P君とQ君のどちらかが先に
一つでも宝を見つけるとそこでゲーム終了となる
縦方向の探査をｎ、横方向の探査をｎ＋１として
宝の個数をｋと置くと、調査する全範囲は
｛ｎ（ｎ＋１）－１｝－（ｋ－１）＝ｎ（ｎ＋１）－ｋと考えられる
Ω＝｛ｎ（ｎ＋１）－ｋ）|ｎ≧２，ｎ（ｎ＋１）－１＞ｋ≧１｝
■縦方向に探査をするP君の確率空間は
Ω＝｛（i，j）|１≦i≦ｎ，１≦j≦ｎ（ｎ＋１）－ｋ｝から
#A＝ｎ｛ｎ（ｎ＋１）－ｋ｝－｛ｎ（ｎ＋１）－ｋ－１｝（ｎ－１）
＝ｎ（ｎ^2＋ｎ－ｋ）－｛ｎ（ｎ^2－１）－ｋ（ｎ－１）－（ｎ－１）｝
＝ｎ^3＋ｎ^2－ｋｎ－ｎ^3＋ｎ＋ｋｎ－ｋ＋ｎ－１
＝ｎ^2＋２ｎ－ｋ－１
#Aは事象Aに含まれる要素の個数
■横方向に探査をするQ君の確率空間は
Ω＝｛（i，j）|１≦i≦ｎ＋１，１≦j≦ｎ（ｎ＋１）－ｋ｝から
#B＝（ｎ＋１）｛ｎ（ｎ＋１）－ｋ｝－ｎ｛ｎ（ｎ＋１）－ｋ－１｝
＝｛ｎ（ｎ＋１）^2－ｋ（ｎ＋１）｝－｛ｎ^2（ｎ＋１）－ｋｎ－ｎ｝
＝｛ｎ^3＋２ｎ^2＋ｎ－ｋｎ－ｋ｝－｛ｎ^3＋ｎ^2－ｋｎ－ｎ｝
＝ｎ^2＋２ｎ－ｋ＝ｎ（ｎ＋２）－ｋ
#Bは事象Bに含まれる要素の個数
■［ｎ≧２，ｎ（ｎ＋１）－１＞ｋ≧１］の条件下で以下の式が成立する
∴P(A)＝｛ｎ（ｎ＋２）－ｋ－１｝／｛ｎ^2（ｎ＋１）－ｋｎ｝
∴P(B)＝｛ｎ（ｎ＋２）－ｋ｝／｛ｎ（ｎ＋１）^2－ｋ（ｎ＋１）｝
Ωの部分集合を事象と言う
Ω自身は全事象と言う
最初に探す方向を i
行または列が変わる時を ｊ として
P君とQ君のうちどちらが先に宝を見つけるのかという
事象Aと事象Bを考える.
A＝｛（i，j）| i または j が宝｝
B＝｛（i，j）| i または j が宝｝
別スレでは等確率とデタラメ書いてたよなぁ。
https://rio2016.5ch.net/test/read.cgi/math/1540218853/87
読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。
すでに正解とPCでのカウントの照合が終わっているのに
読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。
とすると、部屋数が増えたり宝が増えても数式として算出可能なのでしょうか？
なんだ、このスレでデタラメ書きつづけてたのかよｗ
ご返答ありがとうございます。
もちろんそういう近づけかたを考えてもいいけど、その近づけかたを離散位相といってはいけない。
流石にこの程度の基本的な単語は正確に意味を確認するようにしないといかん。
>>152は無視してな
完全に寝ぼけてたわ
すみません。発言がよくわかりません。
「数式として算出」とは？ >>166は数式ではない？
ここでいう数式とは多項式などのΣのない形のものでしょうか？
「部屋数が増えたり」も、もともと部屋をn×(n+1)などとしていて大きさを変えられますよ？
宝の数が2以外でも（3なら3と）固定されていて部屋の形がn×(n+1)または(n+1)×nなら>>161を少し変えれば
Σのないnについての多項式が得られる、とは言えます。
念のため書いておくと
式>>166は部屋の縦横、宝の数が任意だが、Σがある。
部屋の形がn×(n+1)または(n+1)×nとして適当に場合分けすることにより数式処理ソフトで
Σの計算できるようにした、そのコードが>>161。
高校生にもわかるようにどなたか解説していただけませんか？(入試数学の解答のような形式であればありがたいです)
可能だけど立式するのは結構な手間
Σを用いた式として立式して sagemath というソフトで簡略化して n の多項式にしたのが>>161
どう考えても面倒なので161さん以外の誰もやっていなかった
宝箱の数をkとして立式することは可能だろうけれども、
更なる面倒さに付き合ってくれる人がいなければここには書かれない
こんなとこでどうでしょうか。
僕も今初めて使うので適当ですが例えば>>161の最初の
P1の式(Σを含むもの)の簡略化などはiPhoneアプリでも以下のようにして行えました
アプリ起動して「+」ボタンで新しい式を入力するモードにして
var l,a,n,k = var('l','a','n','k')
a=m/2
m=n+1
sum((m-1)*(m)-2*l-1, l,1,m-2)
+ sum(sum((m-1-k)*(m-k)+k-1-l, l,k,2*k) + sum((m-1-k)*(m-k)+3*k-2*l-1, l,2*k+1,m-2), k,1,a-2)
+ sum(sum((m-1-k)*(m-k)+k-1-l, l,k,m-2), k,a-1,m-2)
と入力して「evaluate」ボタンでこの式を評価(簡略化)
https://i.imgur.com/4bpYZLg.jpg
http://codepad.org/pbmWeZZ5
少し説明を加えておくと、マスに1から順に番号をあたえます。
配列P[c]には、c番目 のマスをPは何番目に調査するか
配列Q[c]には、c番目 のマスをQは何番目に調査するか を入れておきます。
i番目とj番目のマスに宝があるとき、P[i]とP[j]を比べて小さい方の値で、Pは宝を発見し、
Q[i]とQ[j]を比べて小さい方の値で、Qは宝を発見します。
この値を比べ、PとQどちらが早く発見したかを判定すると言うだけのものです。
縦、横のマスの数の変更や、宝の数の変更も難しくないと思うので、興味がある方はどうぞ。
（本当は、配列は一つで十分なんだけど、可読性や対称性を考えて書いておきました）
Pを通り方向ベクトル(1,2,0)に平行な直線lとBとの共有点を考えるとき、以下の問いに答えよ。
（1）lとBの共有点の個数を場合を分けて答えよ。
（2）共有点の個数が2個のときを考える。共有点の一方をS、他方をTとする。
P(x,y,z)とするとき、長さの積PS・PTをx,y,zで表せ。Pが表面上にあるときはP=Sとして考えよ。
（3)Bを平面x=u(-1≦u≦1)で切った切断面D_u上を点Pが動く。P(u,y,z)においてy^2+z^2の取りうる最大値Mをuで表せ。
さらにz=0のとき、積分 I_u = ∫[0→M] (PS・PT) dy をuで表せ。
（4）（3）で求めたI_uに対して定積分K = ∫[-1→1] I_u du を求め、さらに比の値K/(4π/3)を求めよ。
次の和を求めよ。
Σ[k=1,2,...,n-1] {(n,k)・(n+1,k-1)}
点Pが(a,b)にあるとき、偶数の目が出たら(a+1,b+1)に移動させ、奇数の目が出たら(a+1,b-1)に移動させる。
このとき、以下の事象が起こる確率を求めよ。
（1）Pが半直線y=x(x≧1)の上に乗る。
（2）Pが直線y=2x+1の上に乗る。
（3）m,nを整数の定数とし、Pが半直線y=mx+n(x≧1)の上に乗る。必要があればm,nの値に応じて場合分けして答えよ。
レスありがとうございました。
多項式で与えられたので他のソフトでも>163のように
簡単にグラフ化できました。
そういう意味で数式と書いたつもりでした。
いつもｃのコードありがとうございます。
このコードだと縦横マスを増やすのは容易でも、宝の数を増やすには for loopを
for(i=1,Pwin=Qwin=Draw=0;i<mn;i++)for(j=i+1;j<mn;j++) for(k=j+1;k<mn;k++) for(l=k+1;l<=mn;l++)
という具合に増やす必要がありますよね？
for(i=1,Pwin=Qwin=Draw=0;i<mn-2;i++)for(j=i+1;j<mn-1;j++) for(k=j+1;k<mn;k++) for(l=k+1;l<=mn;l++)
で、空回りを回避してます。
もし、このアルゴリズムで、宝の数を一般数化するなら、i,j,k,...の変数を配列にしてループにいれるか、
再帰関数化するか等の方がスマートですが、二つで固定なら、提示したような感じがシンプルですね。
しかし、宝の数可変を前提にプログラムを組むなら、別の方策を取ります。
Qは時刻 c に最初の宝を見つけるので、
・Pの宝の発見時刻が全てcより大きい → Qの勝ち
・Pの宝の発見時刻にcを含み、残りは全てcより大きい → 引き分け
・それ以外 → Pの勝ち
です。
Qは、時刻cに、マスcを調査するので、マスc+1、c+2、...の中に、P[x]>c を満たす
マスがいくつあるかをあらかじめカウントし、テーブル化すれば、あとは、
二項係数の積の和だけの、プログラムとなると思います。
m+1≧n≧1 のとき
Σ[k=1，n] C(n，k) C(m，k-1) = C(m+n，n-1)
∵ (1+x)^n (1+x)^m を展開したときの x^(n-1) の係数だから。
Σ[k=1,n-1] C(n,k) C(n+1,k-1) = C(2n-1,n-1) - C(n+1,2)
蝉「おまえさ、人としじみのどっちが偉いか知ってるか？」
伊坂幸太郎「グラスホッパー」角川文庫 (2007)
部屋の数=mn、宝の数trでmnＣtr個の組み合わせを返すサブルーチンが必要になって、ここが処理のボトルネックになるんじゃないかと思うのですが。
何が数学の本質なんでしょうか？
論理的な体系の構築？ 定理の創出？
>>60
まず、部屋を探る順番が一般の場合を考える。
部屋がNあり、その集合をＲとする。Ａ君、Ｂ君が探る順番を表わす全単射写像をそれぞれf,gとする：
f,g: Ｒ→{0,1,…,N-1} （順番は0から始まるとする。）
部屋自体の位置はなんら答えに影響しない。
σ=g・f^{-1} と置くと、σは{0,1,…,N-1}の置換。（・は写像の合成）
Ａ君がi番目に探る部屋はＢ君がσ(i)番目に探る部屋ということ。
以下、「Ａ君がi番目に探る部屋」のことを「部屋i」ということにする。
求めたいのは、「Ａ君がＢ君よりも早く宝を見つける宝の配置の数」であるが、宝の数をcとすると、それは
Σ[σ(i)>i] binomial(#{j| j>i, σ(j)>i}, c-1) （0≦i,j≦N-1、binomialは二項係数）
である。
なぜか？
「Ａ君が初めて宝を見つける部屋（部屋iとする）」で場合分けしよう。
（つまり部屋0～i-1には宝がなく、部屋iに宝がある場合）
部屋iはＢ君がσ(i)番目に探る部屋だからσ(i)>iでないと
少なくともＢ君はＡ君よりも前か同時に部屋iで宝を見つける
（Ｂ君はその前に別の部屋で宝を見つけることもある）ことになりＡ君は勝てない。
したがって、σ(i)>iが必要。
残りのc-1個の宝は部屋i+1～N-1にあるが、宝がある部屋を部屋jとすると、
やはりσ(j)>iでないといけない。逆に全部の宝でそうであればＡ君が勝つ。
よって、残りのc-1個の宝が置かれてもいい部屋の数は#{j| j>i, σ(j)>i}だけあり、
全部そこに置かれる場合はbinomial(#{j| j>i, σ(j)>i}, c-1)通り。
したがって、上記のようになる。
続く
続き
部屋が縦m、横nで、Ａ君は横1行を探し終えたらすぐ下の1行に移り、
Ｂ君は縦1列を探し終えたらすぐ右の1列に移るという場合を考える。
つまり、m=4,n=3の場合、Ａ君は
０１２３
４５６７
８９1011
Ｂ君は
０３６９
１４７10
２５８11
という順番で探す。
このとき、σ=0,3,6,9,1,4,7,10,2,5,8,11。
一般には、σ(nk+l)=ml+k (0≦k≦m-1, 0≦l≦n-1)。
ここまでをPythonで表すと：
#二項係数。SageMathでは定義ずみ
def binomial(n,r):
from math import factorial as f
return f(n)//f(r)//f(n-r) if r>=0 and n-r>=0 else 0
#置換p、宝c個で勝つ宝の配置の数
def nwinperm(p,c):
N = len(p)
return sum(binomial(len([j for j in range(i+1, N) if i<p[j]]),c-1)
for i in range(N) if i<p[i])
#部屋が縦m、横nのときの置換
def rectperm(m,n):
return [m*l+k for k in range(m) for l in range(n)]
#部屋が縦m、横n、宝がc個で横優先が勝つ宝の配置の数
def nwinrect0(m,n,c):
return nwinperm(rectperm(m,n),c)
続く
続き
部屋が縦m、横nの場合を考えているが、もう少し計算を進める。
#{j| j>i, σ(j)>i} をこの場合に具体的に表そう。
i,j (0≦i,j≦mn-1)をそれぞれ nk+l, nk'+l' (0≦k,k'≦m-1, 0≦l,l'≦n-1) とする。
σ(i)>i ⇔ lm+k>nk+l ⇔ (m-1)l>(n-1)k、
j>i ⇔ nk+l>nk'+l' ⇔ 「k=k' かつ l<l'」または「k<k'」、
σ(j)>i ⇔ l'm+k'>nk+l ⇔ l' + k'/m > (nk+l)/m [ここで nk+lをmで割った商をq、余りをrとすると]
⇔ l' + k'/m > q + r/m ⇔ 「q≦l'≦n-1 ただし l'=q, k'≦r を除く」
を使って
#{j| j>i, σ(j)>i} = #{(k',l')|『「k=k' かつ l<l'」または「k<k'」』かつ l'm+k'>nk+l}
に出てくる『「k=k' かつ l<l'」または「k<k'」』かつ l'm+k'>nk+lを満たす組(k',l')の数を求める。
k=k' かつ l<l'のとき σ(i)>iからlm+k>nk+lだからl'm+k'>nk+lは常に成り立つので、l<l'≦n-1でn-1-l個。
k<k' のとき l'm+k'>nk+l ⇔ 「q≦l'≦n-1, k<k'≦m-1 ただし l'=q, k<k'≦r を除く」だから
(n-q)(m-1-k) - (r-k)δ(r>k)個、ただしδ(P)はPが真なら1、偽なら0である関数。
よって、#{j| j>i, σ(j)>i} = (n-1-l) + (n-q)(m-1-k) - (r-k)δ(r>k)。
したがって、求める数は
Σ[0≦k≦m-1, 0≦l≦n-1, (m-1)l>(n-1)k] binomial((n-1-l) + (n-q)(m-1-k) - (r-k)δ(r>k), c-1)。
これを使ったPythonコード：
#nloc(m,n,k,l)は縦m、横nの部屋で横優先が部屋(k,l)で初めて宝を発見する場合で
#宝が置かれても縦優先に先を越されない部屋の数。
def nloc(m,n,k,l):
q,r = divmod(n*k+l,m)
return (n-1-l) + (n-q)*(m-1-k) - (r-k if r > k else 0)
#部屋が縦m、横n、宝がc個で横優先が勝つ宝の配置の数
def nwinrect1(m,n,c):
return sum(binomial(nloc(m,n,k,l),c-1) for k in range(m) for l in range(n) if (m-1)*l>(n-1)*k)
続く
続き
部屋がm×(m+1) (n=m+1) のとき。
(m-1)l>(n-1)k ⇔ 0≦k≦m-2 かつ k+1≦l≦m。
(nk+l)/m = k + (k+l)/m より k+l<mのときq=k,r=k+l、k+l≧mのときq=k+1,r=k+l-m。
r>k (k+l<m)とr≦k (k+l≧m)とに分けるように場合分けをする：
①0≦k≦[(m-1)/2], k+1≦l≦m-k-2 のとき r>k、
②[(m+1)/2]≦l≦m-1, m-1-l≦k≦l-1 または ③l=m, 0≦k≦m-2 のとき r≦k。
m=6のとき
×①①①①②③
××①①②②③
×××②②②③
××××②②③
×××××②③
×××××××
m=7のとき
×①①①①①②③
××①①①②②③
×××①②②②③
××××②②②③
×××××②②③
××××××②③
××××××××
後はΣの計算。>>60に合わせるとＱ君がＡ君の立場でmが>>60でのn。
#以下 SageMathコード
,var m,n,l,k,q,r,c
T2 = (n-1-l) + (n-q)*(m-1-k)
T1 = T2 - (r-k)
#mが奇数の場合：
Q1 = (sum(sum(binomial(T1.subs({n:m+1,q:k,r:k+l}),c-1), l,k+1,m-k-2), k,0,(m-1)/2-1)
+ sum(sum(binomial(T2.subs({n:m+1,q:k+1,r:k+l-m}),c-1), k,m-1-l,l-1), l,(m+1)/2,m-1)
+ sum(binomial(T2.subs({n:m+1,l:m,q:k+1,r:k}),c-1), k,0,m-2)
).subs({m:n,c:2}).simplify_full().factor()
#mが偶数の場合：
Q2 = (sum(sum(binomial(T1.subs({n:m+1,q:k,r:k+l}),c-1), l,k+1,m-k-2), k,0,m/2-2)
+ sum(sum(binomial(T2.subs({n:m+1,q:k+1,r:k+l-m}),c-1), k,m-1-l,l-1), l,m/2,m-1)
+ sum(binomial(T2.subs({n:m+1,l:m,q:k+1,r:k}),c-1), k,0,m-2)
).subs({m:n,c:2}).simplify_full().factor()
def Q1st(x):
return (Q1 if mod(x,2) == 1 else Q2).subs({n:x})
続く
続き
部屋がm×(m-1) (n=m-1) のとき。
(m-1)l>(n-1)k ⇔ 1≦l≦m-2 かつ 0≦k≦l。
(nk+l)/m = k + (l-k)/m より q=k,r=l-k。
r>k (l>2k)とr≦k (l≦2k)とに分けるように場合分けをする：
①0≦k≦[(m-3)/2], 2k+1≦l≦m-2 のとき r>k、
②1≦k≦[(m-3)/2], k≦l≦2k または ③[(m-1)/2]≦k≦m-2, k≦l≦m-2 のとき r≦k。
m=7のとき
×①①①①①
×②②①①①
××②②②①
×××③③③
××××③③
×××××③
××××××
m=8のとき
×①①①①①①
×②②①①①①
××②②②①①
×××③③③③
××××③③③
×××××③③
××××××③
×××××××
後はΣの計算。>>60に合わせるとＰ君がＡ君の立場でmが>>60でのn+1。
#mが偶数の場合：
P1 = (sum(sum(binomial(T1.subs({n:m-1,q:k,r:l-k}),c-1), l,2*k+1,m-2), k,0,m/2-2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,2*k), k,1,m/2-2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,m-2), k,m/2-1,m-2)
).subs({m:n+1,c:2}).simplify_full().factor()
#mが奇数の場合：
P2 = (sum(sum(binomial(T1.subs({n:m-1,q:k,r:l-k}),c-1), l,2*k+1,m-2), k,0,(m-3)/2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,2*k), k,1,(m-3)/2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,m-2), k,(m-1)/2,m-2)
).subs({m:n+1,c:2}).simplify_full().factor()
def P1st(x):
return (P1 if mod(x,2) == 1 else P2).subs(n=x)
以上、整理して少し異なったけど>>161の導出でした。
P1 == 1/24*(6*n^3 + 20*n^2 - n - 27)*(n - 1) # nが奇数のとき
P2 == 1/4*n^4 + 7/12*n^3 - 7/8*n^2 - 13/12*n + 1 # nが偶数のとき
Q1 == 1/24*(6*n^2 + 10*n - 3)*(n + 1)*(n - 1) # nが奇数のとき
Q2 == 1/24*(6*n^2 - 2*n - 5)*(n + 2)*n # nが偶数のとき
それだけ前置きやってｋを含めた式が作れないのですか？
「ｋを含めた式」って何？
あと>>178見て
宝の数が任意のものならΣが取れないでしょう。
正直、Pythonはわからないし、sagemathは数学そのものなのでまだ理解できていないので
読めたのはＣだけですが…
久しぶりにまともなスレになった気がします
P(A)をP(B)で割ることによって
P君の勝つ数とQ君の勝つ数が導ける
P(A)／P(B)＝（P君の勝つ数）／（Q君の勝つ数）
｛ｎ（ｎ＋２）－ｋ－１｝／｛ｎ^2（ｎ＋１）－ｋｎ｝
P(A)／P(B)＝――――――――――――――――――――
｛ｎ（ｎ＋２）－ｋ｝／｛ｎ（ｎ＋１）^2－ｋ（ｎ＋１）｝
＝（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝
∵［ｎ≧２，ｎ（ｎ＋１）－１＞ｋ≧１］
∵の範囲でｎとｋをいろいろと変えて見ることにより
様々な勝率が導ける
計算知能にそのまま入力するだけで通分と約分を
自動計算してくれるので試してごろうじろう
■Wolfram入力例
（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝，ｋ＝２，ｎ＝３
AからBCへ下した垂線をAD,
BCの中点をMとする。
BD > CDとすると
BD^2 - CD^2 = 2BC･MD を示せ。
中線定理やろ
「垂線の定理」と言うらしい
ナルシストやね
⇔
BD-MD=CD+MD
こんな感じか
BD^2 - CD^2 = (BD+CD) (BD-CD)
= BC ( (BM+MD) - (CM - MD) )
= 2BC･MD
(∵ BM=CM )
特に垂線である意味がないし「垂線の定理」って何かの間違いでは？
誰が面白い‥‥
AB^2 - AC^2 = (BD^2 + h^2) - (CD^2 + h^2) = BD^2 - CD^2 = ... = 2 BC・MD
なるほど
a+b^2=c^3
a^2-b(b+c)=a+b+c
（1）このような(a,b,c)を一組求めよ。
（2）（1）で求めたもの以外に(a,b,c)の組が存在するなら、全て決定せよ。
Prelude> [(a,b,c)|a<-[1..100],b<-[1..100],c<-[1..100],a+b^2==c^3,a^2-b*(b+c)==a+b+c]
[(4,2,2)]
[(2,2)]
最初にピタゴラスの定理を証明した人って、どういう発想で定理が正しいと考えたのでしょうか？直感でしょうか、経験的によく知られていたのでしょうか？
a=-1，b=-3，c=2
a=0，b=-1，c=1
a=0，b=0，c=0
a=-1，b=1，c=0
a=4，b=2，c=2
証明はどうだかわからないがピタゴラスが多いついたのはタイルを見て予想したと言われているらしい
直角二等辺三角形を敷き詰めると直角二等辺三角形の場合にはピタゴラスの定理が成り立つことがすぐにわかる
ピタゴラスはそこから直角三角形なら常に成り立つのではないかと考えたということのようだ
例えばAtiyahとかdeBrangeさんなど
http://www.wfg-bluebonnet.com/blog/garden-life/%E3%83%94%E3%82%BF%E3%82%B4%E3%83%A9%E3%82%B9%E3%81%AE%E3%82%BF%E3%82%A4%E3%83%AB.2jpg.jpg
そのためにはやはり、数学の全分野だけでなく、物理学とか哲学とか計算機科学の全分野も究めないと無理なレベルでしょうか？
証明は？
a=c^3-b^2
a^2-b(b+c)=a+b+c
a^2-b^2-bc=a+b+c
a^2-c^3+a-bc=a+b+c
a^2-c^3-bc-b-c=0
(c^3-b^2)^2-c^3-bc-b-c=0
の整数解を求める
宝の数を変化させるコードをHaskellに移植してみた。
import Data.List
import Data.List.Split
m = 5 -- 縦マス(短軸）
n = 6 -- 横マス(長軸)
k = 5 -- 宝の数
q = [0..m*n-1]
matQ = chunksOf m q
matP = transpose matQ --行列を転置して
p = concat matP -- 配列に変換
combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
treasure = combinations k q -- 宝の組み合わせ
ip y = minimum $ map(\x -> elemIndices x p!!0) y -- 宝の、配列ｐでのindex列を求めて最小値を返す
iq y = minimum $ map(\x -> elemIndices x q!!0) y
idxp = map ip treasure -- 宝の組み合せで実行して
idxq = map iq treasure
p_q = zipWith (-) idxp idxq -- 差をとって大小判別
p1st = length $ filter (<0) p_q -- 短軸方向探索ｐが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q
main = do
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw
Prelude> :main
p1st = 54036, q1st = 55469, draw = 33001
ネットで調べたけどそれらしい答えが無くて困ってる
無限ホテルが集合論のお話で、ホテルは可算無限集合、無限に居る宿泊客全員も可算無限集合で、どっちも無限としての大きさが合うから部屋は過不足なく用意されるって話だってところまではネットで読んだ
で、Wikipediaには順序数? の計算ルールが書いてあって、1+ωとω+1は違うってあったからこれが部屋移動の理由かと最初は思った
でも無限ホテルって無限人の来客があってもokってあるから、これってω+ωでどこに客をぶちこんでも意味変わらないなと
だからこの予想は違うと今は思ってる
この疑問のしっくり来る(理解できる)解説が見つからなくてずっとモヤモヤしてるので、誰か教えてくれるとありがたいです
import Data.List
import Data.List.Split
m = 5 -- 縦マス(短軸）
n = 6 -- 横マス(長軸)
k = 5 -- 宝の数
q = [0..m*n-1]
matQ = chunksOf n q
matP = transpose matQ --行列を転置して
p = concat matP -- 配列に変換
combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
treasure = combinations k q -- 宝の組み合わせ
ip y = minimum $ map(\x -> elemIndices x p!!0) y -- 宝の、配列ｐでのindex列を求めて最小値を返す
iq y = minimum $ map(\x -> elemIndices x q!!0) y
idxp = map ip treasure -- 宝の組み合せで実行して
idxq = map iq treasure
p_q = zipWith (-) idxp idxq -- 差をとって大小判別
p1st = length $ filter (<0) p_q -- 短軸方向探索ｐが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q
main = do
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw
>matrix.exe
p1st = 55469, q1st = 54036, draw = 33001
先頭を開けて1人追加するのは 1+ω = ω
倍の部屋番号へ移して ω 人追加するのは 2ω = ω
EFG
ｎ＝２の６マスでP君Q君のそれぞれのファーストの
組の総数をお願いします<(_ _)>
Q勝ち：BG CG BC CE CF
引き分け：AB AC AD AE AF AG BE
q win : BC, BF, BG, CF, CG
even : AB, AC, AE, AF, AG, BE
かと思った
DEF
P勝ち CD DE DF EF
Q勝ち BC BE BF CE CF
引分け AB AC AD AE AF BD
>>236に至ってはよく見るとABCDEFGの7種使ってるし
ありがとうございます
連続の質問になって申し訳ないんで付けど、2ωとω+ωってこれは違うものなんですか?
１個目だけじゃなく、２個目の宝を先に見つけることも考えたら
結局、PQで差はないという直感どおりの結果になるな
１２３
４５６
１２ ・Ｑ
１３ ・Ｑ
１４ ・Ｐ
１５ ・Ｐ
１６ ・・
２３ ＱＱ
２４ ・Ｐ
２５ ＱＰ
２６ Ｑ・
３４ ＰＱ
３５ Ｑ・
３６ Ｑ・
４５ ＰＰ
４６ Ｐ・
５６ Ｐ・
https://rio2016.5ch.net/test/read.cgi/math/1537516085/256
＞オイラー定数をγと置く。nの約数の総和をσ(n)と置く。RHは
＞
＞σ(n)＜(e^γ)*n*log(log n) (∀n＞5040)
＞
＞と同値であることが知られている。
とりあえず元論文はコレらしい
[24] G. Robin, Grandes valeurs de la fonction somme de diviseurs et hypoth`ese de Riemann,
J. Math. Pures Appl. 63 (1984), 187–213.
英語で読めるのないかなぁ？
D:\bin>for %i in (1,2,3,4,5,6,7,8,9,10) do treasure 5 6 %i
D:\bin>treasure 5 6 1
p1st = 14, q1st = 14, draw = 2
D:\bin>treasure 5 6 2
p1st = 203, q1st = 197, draw = 35
D:\bin>treasure 5 6 3
p1st = 1801, q1st = 1727, draw = 532
D:\bin>treasure 5 6 4
p1st = 11418, q1st = 11008, draw = 4979
D:\bin>treasure 5 6 5
p1st = 55469, q1st = 54036, draw = 33001
D:\bin>treasure 5 6 6
p1st = 215265, q1st = 211894, draw = 166616
D:\bin>treasure 5 6 7
p1st = 685784, q1st = 680768, draw = 669248
D:\bin>treasure 5 6 8
p1st = 1827737, q1st = 1825076, draw = 2200112
D:\bin>treasure 5 6 9
p1st = 4130886, q1st = 4139080, draw = 6037184
D:\bin>treasure 5 6 10
p1st = 7995426, q1st = 8023257, draw = 14026332
１：同等
１～８：短軸探索有利
９、１０：長軸探索有利
という結果になった。
Haskellのコードはここ
--exe Fileにコンパイルしてコマンドラインから実行できるように改変（但し、エラー処理皆無）
https://rio2016.5ch.net/test/read.cgi/math/1490734993/209
$Rscript main.r
P1st Q1st even
3 4 13
６マスで宝を３個にしてみた
$Rscript main.r
P1st Q1st even
3 4 13
P 1st
[,1] [,2] [,3]
[1,] C C D
[2,] D D E
[3,] E F F
Q 1st
[,1] [,2] [,3] [,4]
[1,] B B B C
[2,] C C E E
[3,] E F F F
even
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,] A A A A A A A A A A B B B
[2,] B B B B C C C D D E C D D
[3,] C D E F D E F E F F D E F
Rのスクリプトをここに置いたから数値を変更して実行可能
http://tpcg.io/X3bC0A
自由ついでに５×６マスで宝が５個、先に全部の宝を見つけた方が勝者とすると
これくらい差が出る
D:\bin>treasure2 5 6 5
p1st = 54036, q1st = 55469, draw = 33001
確率にすると
> treasure2(5,6,5)
P1st Q1st even
54036 55469 33001
> 54036/(54036+55469+33001)
[1] 0.379184
> 55469/(54036+55469+33001)
[1] 0.3892398
なので差があると直観するかどうかは個々人の感性だな。
使える人間にとっては電卓みたいなもんだよ。
log2の計算にいちいちマクローリン展開して手書き計算しないだろ。
juliaが流行しているのは自分の周りだけなのかな
（NGに巻き込まれて見えてないだけだったらゴメン）
どれかを移植して実力を示していただけたらうれしい。
５✕６マスで宝が１５個の時の計算とかまだ誰も出してない。
先に２個の宝をみつけた方なら
123
456
12 Q
13 Q
14 P
15 P
16 =
23 Q
24 P
25 P
26 =
34 Q
35 =
36 =
45 P
46 =
56 =
にならない？
２個を先にみつけるじゃなくて
これは１個めの発見はQの方が確率が高くて、２個めに発見はPの方が確率が高いというだけの話だったみたいね。
juliaが周りで流行ってるだけで自分自身はＣの人（Ｃソース書いてくれた人とは別人）
5x6ますで宝１５個とか、ID:TZqGbv4dにお願いしたらすぐやってくれるんじゃない？
完全に作り直してるし。
i 番目をどちらが先にみつけるかを計算してみた。
４×５マスに宝が５個あるとき
> treasures(4,5,5)
p1st q1st even
[1,] 1948 9680 3876
[2,] 5488 10016 0
[3,] 7752 7752 0
[4,] 10016 5488 0
[5,] 9680 1948 3876
１個め２個めは短軸方向探索のQが、４個め５個めは長軸方向探索のPが、先にみつける宝の配置の組み合わせが多い。３個めは同じ。
全体としてはイーブンだが、
勝者は１個めを先にみつけた方にするか、全部を先にみつけた方にするかで結果が変わる。
Rのコードはここに置いたので数値を変えて実行可能。
http://tpcg.io/Ph7TUQ
「プログラムで、ごり押し計算」
「マクローリン展開して手書き計算」
俺は後者の方が美しく感じるけどな。
実は前者で計算したのに後者を装ってほしいくらい。 (※ 私見です)
単なる確認なんだけれども、
「i番目をどちらが見つけるか」というのは
先にi個見つけた方を勝ちとするのではなくて
例えばi=2だと
Pが１つ発見、Qが１つ発見⇒2番目を見つけたQの勝ち、ということですか？
>254の計算は各人にとってi番めの計算。
例えばi=2だと
Pが１つ発見、Qが１つ発見だと勝敗は未決で
どちらが発見者にとって２個めを発見したらそれが勝者として数えた。
ここまで答が出せた
254 名前：１３２人目の素数さん[sage] 投稿日：2018/10/30(火) 13:03:57.24 ID:TZqGbv4d
>>252
i 番目をどちらが先にみつけるかを計算してみた。
４×５マスに宝が５個あるとき
> treasures(4,5,5)
p1st q1st even
[1,] 1948 9680 3876
[2,] 5488 10016 0
[3,] 7752 7752 0
[4,] 10016 5488 0
[5,] 9680 1948 3876
１個め２個めは短軸方向探索のQが、４個め５個めは長軸方向探索のPが、先にみつける宝の配置の組み合わせが多い。３個めは同じ。
全体としてはイーブンだが、
勝者は１個めを先にみつけた方にするか、全部を先にみつけた方にするかで結果が変わる。
Rのコードはここに置いたので数値を変えて実行可能。
http://tpcg.io/Ph7TUQ
全体を眺めると直感的通り互角。
局所でみると濃淡があるということと理解した。
「真理」の探究は意味があるのでしょうか？
ここまでは算出できたが、宝を１４にしたらエラー終了した。
D:\bin>treasure 5 6 11
p1st = 13346984, q1st = 13395944, draw = 27884372
D:\bin>treasure 5 6 12
p1st = 19312228, q1st = 19372871, draw = 47808126
D:\bin>treasure 5 6 13
p1st = 24301031, q1st = 24358063, draw = 71100756
■引き分けの組み合わせは勝敗と無関係なので除外
宝が２個以上の時、
スタート地点のAマスと対極にある最終マスのLには
P君もQ君もどちらも決してたどり着くことはできないので
このLマスと組みとなる宝の配置は重複情報で意味を持たない
ので除外する
Pが先に見つけるのは以下の21通り
CE,DE,DI,EF,EG,EH,EI,EJ,EK,FG,FH,FI,FJ,FK,GI,GJ,HI,HJ,IJ,IK,JK,
Qが先に見つけるのは以下の22通り
BC,BD,BF,BG,BH,BI,BJ,BK,CD,CF,CG,CH,CJ,CK,DF,DG,DH,DJ,DK,GH,GK,HK,
となる
何マスだろうが、宝が何個であろうが
出発点と終点が同じであれば
PQの宝を得られる個数の期待値は同じということだな
BE と CI の扱いは？
期待値は宝の数なわけで、元の問題は１個めをみつけるステップの数を比較しているんだと思う。
宝を先にみつけたら独り占め、同時にみつけたら折半 というルールなら手に入れる宝の数の期待値は同じになるだろうね。
BI, CIは引き分けで除外なのでは
>>265
メモリを食わないコードを書いてみた
今思ったけど再帰で書いた方が読みやすかったか
https://ideone.com/HaAqJO
>>270
n x n+1 の部屋を縦横に調べる2人の場合は先着する部屋数が等しくなるからそうなるね
p君 ABCDEFGHIJKL
q君 BCDEFGHIJKLA
とかなら殆どq君が独り占め
https://ideone.com/HaAqJO
ありがとうございます、これを待望しておりました。
ABC
DEF
P勝ち CD DE
Q勝ち BC BE CE
勝敗だけ知りたければデータ圧縮が可能
早い者勝ちなら先回りすることが勝つ秘訣
>244に数値を挙げたけど宝の数が増えると逆転しちゃう。
個人的にはどこが逆転する境なのか算出方法が知りたいところ。
質問してばっかりだったので反省して自分で調べてみたんですけど
ω+ω=ω×ω=ω2だってことでした
でも、これは「無限ホテルのω号室の次の部屋からω人の客を泊めた」って事ですよね?
だから無限ホテルの話にあるように1号室→2号室、2号室→4号室、3号室→6号室とずらして、間に入れ込めば2ω=ωになって万事解決って事で合ってますか?
ω×ωはちがうかった……
これじゃω^2になっちゃう
q5..q6..q7..q8
q9q10q11q12
p1..p4..p7..p10
p2..p5..p8..p11
p3..p6..p9..p12
同じ座標なら数字の小さいほうが勝ち
互いに数字の小さいほうを選んで勝負
q2 vs p4 で q2の勝ちとなる
この後にq10とp6の探査をしても
情報としての価値はゼロ
＞だから無限ホテルの話にあるように1号室→2号室、2号室→4号室、3号室→6号室とずらして、間に入れ込めば2ω=ωになって万事解決って事で合ってますか?
そうそう
プログラムを書いてみました。
多倍長を使える処理系を用いればいいのかもしれませんが、実数型で誤魔化しました。
故に大きな数字のところでは誤差があります。
http://codepad.org/VN03aiqT
同じ方針のものがPythonで>>194-198にある
読めよ。数学板なんだから。
いつもありがとうございます。
いやぁ、この出力は圧巻ですね。
Haskell先生もびっくり。
ありがとうございます
おかげさまですごくしっくりきました
失礼しました。
数列を無理矢理分数式化する人や、価値の無い長い文章を投下する人がいるので、
読み飛ばしていました。
宝箱が二つの場合は、多項式での表現が完成していたんですね。
あのようなσやδを含む式を整理する数式処理ツールがあったとは驚きです。
二個で可能だったのだから、もっと多くの場合でも、可能なんでしょうね。
P2 == 1/4*n^4 + 7/12*n^3 - 7/8*n^2 - 13/12*n + 1 # nが偶数のとき
Q1 == 1/24*(6*n^2 + 10*n - 3)*(n + 1)*(n - 1) # nが奇数のとき
Q2 == 1/24*(6*n^2 - 2*n - 5)*(n + 2)*n # nが偶数のとき
多項式ってこれだけ？
ｋは変えられないし出力は意味不明だしナニコレ？
>>204の式ならｋ＝５５４２２２，ｎ＝３２２３００９８８とかでも
数秒で出力してくれるよ
でn=3の場合66通り全部書きだして比較してみろよ。
実際書き出してみた正解とひとつも合わない式になんの意味がある？
Prelude Data.Ratio> print [(n+1)*(n^2+2*n-1-k)%(n^2*(n+2)-n*k)|let n = 3,k<-([0..14]++[16..30])]
[56 % 45,26 % 21,16 % 13,11 % 9,40 % 33,6 % 5,32 % 27,7 % 6,8 % 7,10 % 9,16 % 15,1 % 1,8 % 9,2 % 3,0 % 1,8 % 3,2 % 1,16 % 9,5 % 3,8 % 5,14 % 9,32 % 21,3 % 2,40 % 27,22 % 15,16 % 11,13 % 9,56 % 39,10 % 7,64 % 45]
Prelude Data.Ratio>
kに0～30何入れても正解なんぞ出てこんやろ？
ｋに５００～８００００だとどうですか？
でも出力できたよ
ためしてごろうじろう
k>15だとすべて4/3より大きい値しかでないからアウト。何入れてもだめ。
>>292
n = 3～100までいれて全滅の式にそんな値いれても糞の意味もない。
わかればいいと思う
ｋ＝５７２３４５７７５４２９９７４７２１２，ｎ＝３２１２３０１０９８８５５でも
出力できたぞ
>>204
> ＝（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝
宝箱の数 k=1のとき p の勝ち数 = q の勝ち数になるけど、
上記の式は
= (n+1)(n^2+2n-2) / {n^2(n+2)-n}
= (n^3+3n^2-2) / (n^3+2n^2-n)
だから間違ってるね
というか式の導出過程がどの1ステップも論理的じゃないから検算する必要もないんだけど
>あのようなσやδを含む式を整理する数式処理ツールがあったとは驚きです。
σを処理できないから>>195-196、δを処理できないから>>196-197を人手で行っている
SageMathにやらせているのはn乗の和の公式さえあれば高校生ができる計算
>二個で可能だったのだから、もっと多くの場合でも、可能なんでしょうね。
>>196-197のsubs({m:n+1,c:2})の2を3に変えれば宝が3個の場合の多項式が得られる
ｋ＝１７４５６６１９２５１，ｎ＝１３２１２３でちゃんと１が出力される
さすが
あほか？n=3～100で正しい数値出してない式になんの信憑性がある？
正しい答え出なきゃなんの意味もない。
>>195-196のPythonをHaskellにすればいい
Haskellにもリスト内包表記があるんだから
読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。
と書いたが犠牲者が出ているようだな。
そのまま消えてしまうだけで別にかまわないのに。
>>204
＞計算知能にそのまま入力するだけで通分と約分を
＞自動計算してくれるので試してごろうじろう
＞
＞■Wolfram入力例
＞
＞（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝，ｋ＝２，ｎ＝３
ってわざわざ全角で書いてコピペで入力できなくできないようにしてくれてる所。
wolfram 日本語版だけは全角でも入力できるけどその他のツールは全滅。
いちいち半角に打ち直さんといかん。
脳みそ1ccしかないんちゃうかと。
俺は>205の助言の意味が分かったので>204のidを速攻でNGidに登録したよ。
日本製のエディタには全角半角変換できるのがあるよ。
例えば、http://mana.ikuto.com/
x∈[-3,3]で動くとき最大値Mと最小値mを, aについて次の2つ場合分けすることによって与えよ。
(1)a≧??のとき,
x=3でM=??a-??,
x=-3でm=-??a-??
(2)a≦??のとき,
x=-3でM=-??a-??,
x=3でm=??a-??
となっているのですが、これで場合分けは足りているのですか？
あとで数式コピペしてソフトに貼り付けるなんて普通にするじゃん。
* はさすがに見苦しいから我慢するけど、全部大文字にするのは意味わからん。
しかも
＞計算知能にそのまま入力するだけで通分と約分を
＞自動計算してくれるので試してごろうじろう
といいながらだよ？
アホじゃね？
御助言にしたがってHaskelｌに移植しました。
import System.Environment
choose (n,r) = product[1..n] `div` product[1..n-r] `div` product[1..r]
nloc m n k l = do
let q = div (n*k+l) m
r = mod (n*k+l) m
in (n-q)*(m-k) + q-1-l + if r>k then k-r else 0
nwin m n c = sum[choose ((nloc m n k l), c-1) | k<-[0..m-1], l<-[0..n-1], k*(n-1) < l*(m-1)]
mwin m n c = sum[choose ((nloc n m k l), c-1) | k<-[0..n-1], l<-[0..m-1], k*(m-1) < l*(n-1)]
draw m n c = choose(m*n,c) - nwin m n c - mwin n m c
main = do
argList <- getArgs -- m : 縦マス(短軸） n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
putStrLn $ "p1st = " ++ show(mwin m n k) ++ ", q1st = " ++ show(nwin m n k) ++ ", draw = " ++ show(draw m n k)
おかげ様でこういうのも瞬時に計算してくれました。
１０×２０マスで宝が１００個
>takara 10 20 100
p1st = 15057759425309840160151925452579572328997602171271937639470, q1st = 15057796557877993527038542474310161591275806044157319150135, draw = 60432921540347294111327092128863840691952977587098698541050
数学板は例外かもしれないが、マクロウイルスが貼られるのの予防か半角で投稿すると拒絶されることがあるな。
httpを貼ろうとするとはねられるときには全角にすることもあるな。まあ、数文字大文字に留めるけど。
誰も答えていないしみんな困ってるんだと思うが、すべての場合を調べているわけではない、と考えればいいだけの話。
というかそうとしか捉えられないw
既約分数で表示してくれ
P(A)をP(B)で割ることによって
P君の勝つ数とQ君の勝つ数が導ける
P(A)／P(B)＝（P君の勝つ数）／（Q君の勝つ数）
｛ｎ（ｎ＋２）－ｋ－１｝／｛ｎ^2（ｎ＋１）－ｋｎ｝
P(A)／P(B)＝――――――――――――――――――――
｛ｎ（ｎ＋２）－ｋ｝／｛ｎ（ｎ＋１）^2－ｋ（ｎ＋１）｝
＝（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝
∵［ｎ≧２，ｎ（ｎ＋１）－１＞ｋ≧１］
∵の範囲でｎとｋの数値をいろいろと変えることにより
様々な勝率が導ける
計算知能にそのまま入力するだけで約分を
自動計算してくれるので試してごろうじろう
■Wolfram入力例
（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝，ｋ＝２，ｎ＝３
(n+1)(n^2+2n-1-k)/｛n^2(n+2)-nk｝,k=2,n=3
スタート地点のAマス以外のすべてのマスに
宝がある状態であるｋ＝ｎ（ｎ＋１）－１の時、
必ずP(A)／P(B)＝１になる
ｋ＝ｎ（ｎ＋１）－１の時にP(A)／P(B)≠１となるｎを
見つけることができれば反例になる
見つけてごろうじろう
このとき、以下の不等式を満たす(a,b,c)が存在するような自然数Nの最大値を求めよ。
N≦a^2+b^2+c^2≦2018
タイプミスで draw が間違ってますよ
ご指摘ありがとうございました。
× draw m n c = choose(m*n,c) - nwin m n c - mwin n m c
○ draw m n c = choose(m*n,c) - nwin m n c - mwin m n c
ご指摘を受けたのでデバッグしたのを投稿します。
import System.Environment
choose (n,r) = product[1..n] `div` product[1..n-r] `div` product[1..r]
nloc m n k l = do
let q = div (n*k+l) m
r = mod (n*k+l) m
in (n-q)*(m-k) + q-1-l + if r>k then k-r else 0
nwin m n c = sum[choose ((nloc m n k l), c-1) | k<-[0..m-1], l<-[0..n-1], k*(n-1) < l*(m-1)]
mwin m n c = sum[choose ((nloc n m k l), c-1) | k<-[0..n-1], l<-[0..m-1], k*(m-1) < l*(n-1)]
draw m n c = choose(m*n,c) - nwin m n c - mwin m n c
main = do
argList <- getArgs -- m : 縦マス(短軸） n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
putStrLn $ "p1st = " ++ show(mwin m n k) ++ ", q1st = " ++ show(nwin m n k) ++ ", draw = " ++ show(draw m n k)
変な問題だけど、次の2つの場合、すなわち
・x=3のとき最大、x=-3のとき最小 (a≧6のときか？)
・その逆 (a≦-6のときか？)
に分けて??を埋めよという問題なのだろうから、
その2つのときだけ考えて答えれば良いのではないだろうか
「分けて」ってのが変だよね
次の2つの場合について、ならわかるんだけど。
どうすればこれを実現できるのでしょうか？
自殺をしても無駄なのでしょうか？
行列の固有値と対角化
(4)が全然わかりません
よろしくお願いします
https://i.imgur.com/fxbCChT.jpg
縦５横６のとき宝の数を１から３０まで増やして長軸探索が先にみつける確率と短軸探索がさきにみつける確率の差を描いてみた。
http://i.imgur.com/7qGjOJX.png
縦５横６のときだと宝の数は９から２１のときが長軸探索が有利となった。
短軸有利→長軸有利→同等となるようで、再逆転はないもよう。
縦ｍ横ｍ＋１として長軸探索が有利になる宝の数の上限と下限を算出してみた。
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]
[1,] 0 2 2 6 9 13 17 23 29 36 43 52 61 71 82 93 105 118 132 147
[2,] 0 3 7 13 21 31 43 57 73 88 105 118 135 152 166 185 202 220 242 253
グラフにしてみた。
http://i.imgur.com/PiL9xyH.png
これをお願いします
直感の概算 (a,b,c)＝(40,20,4) N＝2016
微調整 (a,b,c)＝(44,9,1) N＝2018
なんか問題を勘違いしてるかな？
a=44,b=9,c=1のとき2018-a^2-b^2-c^2=0
2018-a^2-b^2-c^2，a=44,b=9,c=1
∴N=2018
(a,b,c)＝(44,9,1)、(43,12,5)
∴N=2018
(43,12,5)
(41,16,9)
(35,27,8)
(34,29,11)
(33,23,20)
a=35,b=28,c=3
a=35,b=27,c=8
∴N=2018
>>324
(34,29,11)は違う
17.27 正則行列Ａ = { [a,0,0] [0,b,c] [0,c,b] } について，次の問に答えよ。(九大*)
(1) 行列Ａの逆行列Ａ^(-1) の (2,3) 成分を求めよ。
(2) Ａの固有値を求めよ。
(3) Ａ^2 = { [4,0,0] [0,0,2i] [0,2i,0] } を満たす a,b,c の値を求めよ。iは虚数単位。
(4) nを自然数とし，Ａ^n (i,j)は行列Ａ^n の(i,j)成分を表わすものとする。
そのとき、Ａ^n (2,2) + Ａ^n (3,2) を n, b, c を用いて表わせ。
a[n] = a^n、b[n] = ((b+c)^n + (b-c)^n)/2、c[n] = ((b+c)^n - (b-c)^n)/2とおいて
A^n = [[a[n],0,0],[0,b[n],c[n]],[0,c[n],b[n]]]、
[1,0,0]A = a[1,0,0]A、[0,1,1]A = (b+c)[1,0,0]、[0,0,1]A = (b-c)[1,0,0]A。
(1) c[-1]。
(2) a,b+c,b-c。
(3) a^2=4 ⇔ a=±2、
(b+c)^2 = 2i、(b-c)^2 = -2i ⇔ (b,c) = (1, i)、(-1, -i)、(i, 1)、(-i, -1)。
(4) b[n] + c[n] = (b+c)^n。
(1)
det(Ａ) = a(bb-cc),
Ａ^(-1) = { [1/a，0，0] [0，b/(bb-cc)，-c/(bb-cc)] [0，-c/(bb-cc)，b/(bb-cc)] }
(2)
det(Ａ-λE) = det{ [a-λ，0，0] [0，b-λ，c] [0，c，b-λ] }
= (a-λ)(b-c-λ)(b+c-λ)
∴ λ = a，b±c,
(3)
Ａ^2 = { [a^2，0，0] [0，bb+cc，2bc] [0，2bc，bb+cc] }
∴ a = ±2，(b,c) = (0,±(1+i)) (±(1+i),0)
(4)
Ａ^n = { [a^n，0，0] [0，f_n，g_n] [ 0，g_n，f_n] }
ただし、f_n = {(b+c)^n + (b-c)^n}/2， g_n = {(b+c)^n - (b-c)^n}/2,
(f_n)^2 - (g_n)^2 = (bb-cc)^n,
あとは自分で考えて
Nの最大値は２０１８
顰蹙のプログラム解
Prelude> [(a,b,c)|a<-[1..45],b<-[a..45],c<-[b..45], a^2+b^2+c^2==2018]
[(1,9,44),(3,28,35),(5,12,43),(8,27,35),(9,16,41),(19,19,36),(20,23,33)]
そこから規則性が見いだせれば理論はあとからついてきたりすることもあるからね。
コラッツの問題みたいに未決のままのもあるけど。
プログラム書いて遊んでるだけ。
数学的な解出てもガン無視してるし。
>312のような指摘もとてもありがたい。
こんなんできた～ってひけらかしたいんだろ？
それを検証して解析解をPCでの計算に応用。
おかげで>142から>314に進化できた。
プログラミングのトレーニング課題を与えてくれた方に深謝。
引き分けのバグ指摘にも感謝。
数理展開が勉強になるようにコードの議論も俺には嬉しい。
このスレではじめてHaskellの存在を知った初心者なので>299のような適格なアドバイスは嬉しいね。
x^2-(4n-3)x+a/(n^2+n+1) ≦ 0
を満たす整数xが存在するためにn,aが満たすべき条件を述べよ。
ａｂｃｄに当てはまる数字は？
※答は１通りしかないようです。
aとcで割れば？
細かい条件は自分でやって
１を許すと沢山ある(1,1,1,1),(1,1,2,1),(1,1,2,2),(1,1,3,1),(1,1,4,1),(1,1,5,1),(1,1,6,1).....けど
(2,2,2,2)が答？
ごめん、ａｂｃｄ は４桁の整数
[2592]
× a^b＋c^d＝1000a+100b+10c+d
○ (a^b)*(c^d)＝1000a+100b+10c+d
これまた顰蹙のダンプリストｗ
Prelude> [(a,b,c,d)|a<-[1..10],b<-[1..10],c<-[1..10],d<-[1..10],a^b*c^d==1000*a+100*b+10*c+d]
[(2,5,9,2)]
失礼しました
>342のコードが正しい
＞(4) nを自然数とし，Ａ^n (i,j)は行列Ａ^n の(i,j)成分を表わすものとする。
＞ そのとき、Ａ^n (2,2) + Ａ^n (3,2) を n, b, c を用いて表わせ。
元の質問者の方向きに解き方の解説
行列のn乗の計算は
A を A’ = P ^-1 A P (A’ は対角行列) と対角化して
A’^n = (P ^-1 A P)^n
⇔ A’^n = P ^-1 A^n P
⇔ P A’^n P^-1 = A^n
ここでA’ は対角行列なので
A’^n は各要素をn乗するだけという流れ
問い (1)～(3) は対角化の仕方を調べているうちにわかると思うので略
x>0
のときの証明方法を教えて下さい
４マス３行（３ターン）と３マス４列（４ターン）で一つの宝と出くわす
確率は同じにならない
■３マス４ターンで少なくとも一つの宝と出くわす確率は
#A＝３^4－２^4＝６５なので
P(A)＝６５／８１
■４マス３ターンで少なくとも一つの宝と出くわす確率は
#B＝４^3－３^3＝３７なので
P(B)＝３７／６４
∴P(A)＞(B)
∵P(A)＝６５／８１＝０．８０２
∵P(B)＝３７／６４＝０．５７８
トランプカードをよくシャッフルする
この山札から１ターン３枚を４回ですべて引くのと
１ターン４枚を３回ですべて引く場合も同じ
しょっちゅう確率の問題に手を出してる。
しかし一度たりとも正解の数値と合ってる式出した事ない。
まぁ本人自分の出した答えが間違ってる事すら理解出来てないのである意味で幸せなのかもしれない。
苦労して立式して合うはずの答えが何故か合わないあの苦々しさに耐えないで済むんだから。
１ターン３枚を４－１回で引く時に
ハートのエース１枚が出る確率は
P(A)＝１－（３／４）（２／３）（１／２）＝３／４
１ターン４枚を３－１回で引く時に
ハートのエース１枚が出る確率は
P(B)＝１－（２／３）（１／２）＝２／３
∴P(A)＞(B)
３列ｘ４ターンと４行ｘ３ターンの探査で同じ確立になるという
計算式をお願いします<(_ _)>
■■■■
■■□■
■■■■
だれかこれをお願いします。
nが自然数なので2次不等式を解いてもあまり上手くいきそうにありません
x = 2n-1で成立
⇔ 4*n^4－2*n^3－4*n+2 ≧ a
ありがとうございます
1つの例ではなくて必要十分な形で占めしていただけませんか
x^2-(4n-3)x+a/(n^2+n+1) ≦ 0 ... (1)
左辺は x=2n-3/2 のとき最小、
整数の範囲では x=2n-1 または x=2n-2 で最小値
(a - 4n^4+2n^3+4n^2-2n) / (n^2+n+1) となる
この式の分母は正なので分子が0以下なら(1)を満たす
立山秀利「入門者のPython」講談社BlueBacks (2018/Sep)
398p．1404円
http://bookclub.kodansha.co.jp/product?item=0000275970
【執筆時に使用した環境】
・Microsoft Windows 8.1 および 10
・Python version 3.6
・Anaconda 5.2 for Windows
・Spyder 3.2.8
上記以外の環境でご利用の場合、本書の解説どおりに操作を行えない可能性があります。予めご了承ください。
本書に掲載されている情報は、2018年8月時点のものです。実際にご利用になる際には変更されている場合があります。
【サポートページ】
http://tatehide.com/bbpython.html
左辺 - 右辺のグラフを書いてみた。
https://i.imgur.com/oCCfFGR.png
4^k+1を4×4^kと見なすことで
成り立つと仮定された不等式を援用して新たな不等式を考えているらしいことはわかりますが
どう計算したら24k-5>0になるのかがわかりません
4^k-（8k+1）
https://i.imgur.com/JzIyFJd.jpg
24k≧72
24k-5≧72-5=67>0
0より大きい24k-5よりもさらに両辺の差は大きいのでもちろんそれは0より大きく、よって不等号の正しさが証明されていたのですね
二回両辺の差を考えようとしてどうしても24k-5が作り出せず混乱していました
ありがとうございました
http://datecocco.hatenablog.com/entry/2015/11/26/000000
はじめて作ったボードゲームを売った話
http://nrmgoraku.hateblo.jp/entry/2017/05/17/210000
ボードゲームイベント「ゲームマーケット」から業界が見えた！
https://entertainmentstation.jp/61107
ゲームマーケットに挑む人向けガイド
http://spa-game.com/?p=4830
ボードゲームはどう作るのか、自分なりに考えた
http://roy.hatenablog.com/entry/2014/07/09/124824
オトナも遊べるボードゲーム！自作するといくらになるのか
http://www.d-laboweb.jp/special/sp312/
自作ボードゲーム販売への道・販売場所編
https://kdsn.xyz/create_game_selling_area/
はじめての同人ボードゲーム作り
https://ameblo.jp/subuta96/entry-11932093993.html
アナログゲーム市場が「クラウドファンディング」で盛り上がるワケ
https://www.sbbit.jp/article/cont1/34394
y'+(x^2)y=1
http://m.wolframalpha.com/input/?i=y%27%2B%28x%5E2%29*y%3D1+for+y
こっちだな
http://m.wolframalpha.com/input/?i=y%27%2B%28x%5E2%29*y%3D1+
□■■■ ■■■■
■■■■ ■■■■ PとQが同時に見つける
■■■■ ■■■□
■■■■ ■■■■ ■■■■ ■■■■ ■■■■
□■■■ ■■■■ ■□■■ ■■■■ ■■■■ Pが先に見つける
■■■■ □■■■ ■■■■ ■□■■ ■■□■
■□■■ ■■□■ ■■■□ ■■■■ ■■■■
■■■■ ■■■■ ■■■■ ■■□■ ■■■□ Qが先に見つける
■■■■ ■■■■ ■■■■ ■■■■ ■■■■
EFGH
IJKL
のように命名すると
縦:m 横: 宝:k での配置を列挙するコードは既出。数値を変えて実行可。
http://tpcg.io/X3bC0A
配置の列挙は確率ではありませんよ
宝が一つの時、縦探査のP君が決して取れない宝は２マス
□□□■
□□□■
□□□□
宝が一つの時、横探査のQ君が決して取れない宝は３マス
□□□□
□□□□
■■■□
決して宝を取れないマスが一マス多いQ君が
P君と同じ確立になるのはなぜ？
P(A)＝１－（３／４）（２／３）（１／２）＝３／４
P(B)＝１－（２／３）（１／２）＝２／３
∴P(A)＞(B)
>>354と同じ
>>370
P（短軸探索）が先、Q（長軸探索）、同時 の配置を表示するスクリプトを書いてみた。
数値を変えて実行できる。
ｍ：短軸 ｎ：長軸 ｋ：宝の数
http://tpcg.io/1woIN1
P1st Q1st even
26 27 13
P (= long axis searcher) finds first.
■ ■ □ ■
□ ■ ■ ■
■ ■ ■ ■
■ ■ ■ □
□ ■ ■ ■
■ ■ ■ ■
■ ■ ■ □
■ ■ ■ ■
□ ■ ■ ■
以下略、
５ｘ６で宝が２個のとき
P1st Q1st even
203 197 35
引き分けになる配置は３５通り、３例ほど挙げるとこんな感じ
とても手作業で列挙する気にはならん。
■ ■ ■ ■ □ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
□ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ □ ■ ■ ■
■ ■ ■ ■ ■ ■
■ □ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ □ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ □ ■ ■ ■ ■
y(x) = u(x)e^(-xxx/3)
を与式に入れると
du/dx = e^(xxx/3),
u(x) = u(0) + ∫[0,x] e^(ttt/3) dt,
y(x) = e^(-xxx/3) {y(0) + ∫[0,x] e^(ttt/3) dt },
かな
この方程式が解けません
教えて下さい
であってますか？
初期条件y(0)=0
当たりをひく確率は同じですか？
複数回ひく場合でも前にひいたくじがなくなる訳ではなく
毎回同じ確率で抽選されるという仮定の場合です
24/9^3で2回当たり
192/9^3で1回当たり
8^3/9^3で全て外れ
「一回でも当たる確率」は、(1+24+192)/9^3 = 217/9^3 = 1-(8/9)^3 < 243/9^3 = 1/3
なので、当選確率1/3のクジを一回引くのより小さい
しかし、「当たる回数の期待値」は
(3*1+2*24+1*192)/9^3=(3+48+192)/9^3=243/9^3=1/3
なので、当選確率1/3のクジを一回引くのと同じ
離散量の確率は場合の数をいかに効率的にカウントするかによるね。
手作業だと漏れがでるからプログラムの利用は必須
>377参照。
投稿前に自分でシミュレーション検証して投稿すれば、
>302のように こいつ 呼ばわりされなくて済むんだけどね。
自分で算出した値が別の言語の算出結果と一致したと投稿されるとシミュレーションの正しさが確認できていいね。
俺が鈍足のRコードのをだすと高速のｃが投稿されたり、解析解が投稿されて数理とプログラム論理の勉強になって嬉しいね。
sed使えよｗ
顰蹙のシミュレーション検証
１００万回シミュレーションして頻度をだしてみた
確率1/3のくじを１回ひく
> mean(replicate(k,sample(x,1))
[1] 0.333435
確率1/9のくじを３回ひくのでは、)
> mean(replicate(k,sum(sample(y,3))))
[1] 0.333176
普通の話は出来れば
普通の話とは出来ればわけておいてほしいな、と思う
この7種類の硬貨を使って1円～70円の70通りの支払いができます
ただし一度に使用できる硬貨は3枚以下(同じ硬貨2度使いは可)です
7種類の硬貨はそれぞれ何円だったのでしょうか？
最低額は1円、最高額は24円以上
不可だと無理なんじゃね？ 知らんけど
３ｘ４の１２マスで宝が一つだけの時、
P君とQ君は互いに最終列と最終行の宝は
取ることができない
□□□■
□□□■
□□□□
□□□□
□□□□
■■■□
つまり、P君の探査範囲は縦３マスｘ３列
Q君の探査範囲は横４マスｘ２行になる
それぞれの探査範囲内でP君とQ君が
少なくとも一つの宝を見つけるという
事象Aと事象Bを考える
P(A)／P(B)＝（P君の勝つ数）／（Q君の勝つ数）
∴P(A)／P(B)＝｛１－｛ｎ／（ｎ＋１）｝^ｎ｝／｛１－｛（ｎ－１）／ｎ｝^（ｎ－１）｝
ｎ＝３のとき、P(A)／P(B)＝３３３／３２０
互いの最終列と最終行にある宝の取れないマスが一つ多い
Q君よりもP君のほうが僅かに確率が上がることが
如実に示される
■Wolfram入力例
｛１－｛ｎ／（ｎ＋１）｝^ｎ｝／｛１－｛（ｎ－１）／ｎ｝^（ｎ－１）｝，ｎ＝３
｛1-｛n/(n+1)｝^n｝/｛1-｛(n-1)/n｝^(n-1)｝,n=3
もちろんいいけど、解が2組以上もあるとは思えないというのが最初の直感
こちとらヒントを見ても、さっぱり分からん
1, 4, 5, 15, 18, 27, 34 かな？
素晴らしい❗
Prelude Data.List> filter (<=70) $ let x = [0,1,4,5,15,18,27,34] in nub$sort$[a+b+c | a<-x,b<-x,c<-x]
[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70]
Prelude Data.List> length $ filter (<=70) $ let x = [0,1,4,5,15,18,27,34] in nub$sort$[a+b+c | a<-x,b<-x,c<-x]
71
(1,４,５)は全くの盲点だったわ
まぁ確率なんか計算できなくても社会生活困るわけでもないし。
なんとなく割り算して、それで楽しいならそれでいいのかなと思う。
先にジョーカーを引いたほうが勝ちとする
(引いたカードは山に戻さない)
①トランプ52枚 + ジョーカー1枚
②トランプ52枚 + ジョーカー2枚
先攻勝率は①②で同じ 27/53
後攻勝率は①②で同じ 26/53
不思議だと思わない？(計算めんどい)
これはクリプテックスの確率だった
>>374よりも精度を上げることができた
-p(a+b)+p^2ab+1=0
pa+b=pb
を満たす整数の組(a,b)が存在するかどうか判定せよ。
y(0)=0
この方程式はどのように解けばいいですか？
答えもお願いします
こうです
y = ∫0→x ∫0→x y dxdx
y(0)=0
y(0)=0
dy(0)/dx=0
初期条件抜けてました
1枚だけで可能 … 1,4,5,15,18,27,34
2枚で可能 … 2,6,8～10,16,19,20,22,23,28,30～33,35,36,38,39,42,45,49,52,54,61,68
残り … 3,7,11～14,17,21,24～26,29,37,40,41,43,44,46～48,50,51,53,55～60,62～67,69,70
与式 → y(0) = 0,
xで微分すると
y '(x) = ∫[0,x] y(t) dt → y '(0) = 0,
もう一度xで微分すると
y "(x) = y(x),
これより
y(x) = a･e^x + b･e^(-x)
y(0) = a+b = 0,
y '(0) = a-b = 0,
∴ a=b=0
y(x) = 0
>>408
初期条件を付記する必要はありません。(式から出ます。)
|x| が小さいところでは、マクローリン展開より
y(x) = e^(-xxx/3){y(0) + ∫[0,x] e^(ttt/3) dt}
= e^(-xxx/3) y(0) + x + Σ[k=1,∞] (-1)^k /{4･7････(3k+1)} x^(3k+1)
これはどうやって出されたのですか？
組み合わせは1,198,774,720通りなので総当たりは断念しました。
表示可能な数字は
[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,72,73,76,79,81,83,86,88,95,102]
で71は無理だな。
私は 398投稿者ではありませんが、次のプログラムで答が出せます。
単純に、C[70,7]通りとはせず、ある程度変数の範囲を押さえ込むことがこつかも。
void next(int depth,int *x){
if(depth<7){
int i,is=x[depth-1]+1,ib=x[depth-1]*3+1,ie=ib<70?ib:70;
for(i=is;i<=ie;i++){x[depth]=i;next(depth+1,x);}
}else{
int i,j,k,mem[211],flag;
for(i=1;i<211;i++)mem[i]=0;
for(i=0;i<7;i++){mem[x[i]]++;mem[2*x[i]]++;mem[3*x[i]]++;}
for(i=0;i<6;i++)for(j=i+1;j<7;j++){mem[x[i]+x[j]]++;mem[2*x[i]+x[j]]++;mem[x[i]+2*x[j]]++;}
for(i=0;i<5;i++)for(j=i+1;j<6;j++)for(k=j+1;k<7;k++){mem[x[i]+x[j]+x[k]]++;}
for(i=1,flag=1;i<71;i++)if(mem[i]==0)flag=0;
if(flag==1)printf("%d,%d,%d,%d,%d,%d,%d\n",x[0],x[1],x[2],x[3],x[4],x[5],x[6]);
}
}
int main(){int x[7];x[0]=1;next(1,x);return 0;}
http://codepad.org/ZhDv5l7X
いつものｃでの高速解をありがとうございます。
他の組合せがないことがわかって助かりました。
ありがとうございました
何度もすみません
これはどんなふうに解けばいいですか？
エネルギー積分 y'^2+2/y=C を求め、変数分離形の公式通りに
∫ dy/√(C-2/y)=∫ dt=t とする
を一つの数式で表すとどうなりますか？
（1）この4つの複素数全てについて、その実部と虚部が共に整数となる例を1組挙げよ。
（2）（1）において、4つの複素数のいずれも0でないことはあるか。
（3）少なくとも1つが実数でなく、その絶対値が1でないようなa,b,c,dはあるか。
間違えた
1 3 7 13 22 34 50 70 95 125 161 203
を一つの数式で表すとどうなりますか？
a_n = Σ(k=1，n+1) [kk/4]
生成関数は x/{(1+x)(1-x)^4}
1, 3, 7, 13, 22, 34, 50, 70, 95, 125, 161, 203, 252, 308, 372, 444, 525, 615, 715, …
http://oeis.org/A173196
Sn=1+3･2+5･2^2・・・+(2n-1)･2^n-1の問題で
Snを掛ける２してSn-2Snで引くまではわかったんですが、
Sn-2Snで引いた-S=1+2･2+2･2^2+・・・ 2･2^2n-1-(2n-1)･2^nから
2+2^2+2^3・・・+2^n-(2n-1)･2^-1←この一番最後にある「-1」がなぜ付いてるのかがわからないんですがなぜこうなってるんでしょうか？
お願いします
おお、これは助かりまする
1 3 7 13 22 34 50 70 95 125 161 203
252 308 372 444 525 615
事前準備と見事に一致
2つの数式で書けば
a_{2m} = m(m+1)(4m+5)/6,
a_{2m+1} = (m+1)(m+2)(4m+3)/6
蛇足だけど。
http://oeis.org/A002623
1つの数式で書けば
a_n = Σ(k=1，n+1) [kk/4] = [ (n+1)(n+3)(2n+1)/24 ],
http://oeis.org/A002623
括弧の位置おかしくね？
r≠1 とする。
S_n = Σ[k=1，n] (2k-1)･r^(k-1),
r･S_n = Σ[k=1，n] (2k-1)･r^k = 1 + Σ[k=1，n+1] (2k-3)･r^(k-1),
辺々引くと
S_n - r･S_n = -1 + 2Σ[k=1,n] r^(k-1) - (2n-1)･r^n
= -1 + 2(r^n - 1)/(r-1) - (2n-1)･r^n,
以下略
2+2^2+2^3・・・+2^n-(2n-1)･2^n-1←でしたすいません
両辺を二乗して
(y "(x))^2 - 1 = 0とし
y "(x) = ±1
を特性方程式を作って後は解くだけでよいでしょうか？
101たす102は3
8たす1は1
では2たす2は？
いろんなサイト見てるんですけど>>417のような問題が無くて解けないんです
よかったらヒント下さい
完全追尾型多項式が完成しました
宝の個数を２で固定します
P1st =｛12n^4+28n^3-42n^2-52n-3(-1)^n+51｝/48
Q1st =｛12n^4+20n^3-18n^2-20n-3(-1)^n+3｝/48
■Wolframに入力すると既約分数表示になるので御注意
P1st／Q1st
=｛8(n-1)｛(n-2)n-6｝/｛2n(n+2)(6n^2-2n-5)-3(-1)^n+3｝｝+1
>>161に正しい答えがあるよ
わかりにくいなら最後の辺りだけ見て
P1st は nが奇数の時P1偶数のときP2
Qも同様
実際に多項式になっていない
C[70,7]通りのRのスクリプトを書いてみた。
正確にはC[70,8]*C[8,3]=528659651520通リｗｗｗ
他の言語に移植する人いるかなぁ？
is.1_70 <- function(x){
total=NULL
for(i in x){
for(j in x){
for(k in x){
ijk=i+j+k
if(!(ijk %in% total)) total=append(total,ijk)
}
}
}
all(1:70 %in% total)
}
M=69
for(a in 0:M){
for(b in a:M){
for(c in b:M){
for(d in c:M){
for(e in d:M){
for(f in e:M){
for(g in f:M){
for(h in g:M){
y=c(a,b,c,d,e,f,g,h)
if(is.1_70(y)) print(y)
}
}
}
}
}
}
}
}
Haskellに移植。
とりあえずコンパイルエラーは出なかった。
朝までに計算が終わるかどうかは不明。
import Data.List
m = 69
sub x = do
let ijk = filter (<=70).nub $ sort [i+j+k| i<-x,j<-x,k<-x]
all (\y -> elem y ijk ) [0..70]
main = do
print $ [(b,c,d,e,f,g,h)| b<-[0..m],c<-[b..m],d<-[c..m],e<-[d..m],f<-[e..m],g<-[f..m],h<-[g..m],sub [0,b,c,d,e,f,g]]
import Data.List
firstUnavailable x = let y = 0:x in head $([1..71] ¥¥)$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
next x = [n:x|n<-[head x+1..firstUnavailable x]]
xss = iterate (¥xs->concat [next x|x<-xs]) [[1]]
isGood x = let y = 0:x in (==70)$length $intersect [1..70]$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
main = do
print [x|x<-(xss !! 6),isGood x]
-- 最終行にhが抜けてたので修正。
import Data.List
m = 69
sub x = do
let ijk = filter (<=70).nub $ sort [i+j+k| i<-x,j<-x,k<-x]
all (\y -> elem y ijk ) [0..70]
main = do
print $ [(b,c,d,e,f,g,h)| b<-[0..m],c<-[b..m],d<-[c..m],e<-[d..m],f<-[e..m],g<-[f..m],h<-[g..m],sub [0,b,c,d,e,f,g,h]
いつもありがとうございます。
お見事に算出されました。
*Main Data.List> :main
[[34,27,18,15,5,4,1]]
import Data.List
firstUnavailable x = let y = 0:x in head $([1..] ¥¥)$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
next x = [n:x|n<-[head x+1..firstUnavailable x]]
xss = iterate (¥xs->concat [next x|x<-xs]) [[1]]
isGood x = (>70) $ firstUnavailable x
main = print [x|x<-(xss !! 6),isGood x]
自宅のパソコンだとghc -O2 で22秒でおわったけどcodepadだとTimeoutした。
この中から1枚だけ選んで10回投げたところ、10回連続で表が出た。
このとき、この選ばれた1円玉が両面表である確率は
普通の1円玉である確率より高い？低い？
ghciでやったんだ。流石にその勇気はなかったｗｗ
6種のコイン2枚以下で34円まで表せないといけない
逆に6種を2枚以下で34円まで表せるなら1枚35円の
コインを追加した7種が3枚までで70円まで表せる。
よってまず6種2枚までで1～34円が全て表せるかを調べて、
それが無理ならコインの価値は最大34円までと限定できる
この先も上から攻めていけば多少探索範囲を限定できると思うが定かではない
＞最大価値のコインの価値が1枚35円以上だと、残りの
＞6種のコイン2枚以下で34円まで表せないといけない
この部分はアプリオリではないか
正しいような気はするが少なくとも数行では証明できなさそう
P(本物|10連続表)
=P(本物＆10連続表)/P(10連続表)
=999/1000・1/1024/P(10連続表)
=999/1024/P(10連続表)/1000
P(偽物|10連続表)
=P(偽物＆10連続表)/P(10連続表)
=1/1000・1/1/P(10連続表)
=1/P(10連続表)/1000
∴ P(本物|10連続表)＜P(偽物|10連続表)
y"(x) = y(x)^(-2)
以下より(x)を省略
(y'^2)' = 2y'y"
(y^(-1))' = -y(x)^(-2) * y'
ここで
y" = y^(-2)の両辺にy'をかけて
y"y' = y(x)^(-2) * y' となり
1/2 * (y'^2)' = - (y^(-1))'
(y'^2)' = (-2y^(-1))'
故に y'^2 = -2y^(-1)
ここまではあってますか？
ここから先が解けるかどうかわかりませんがもう少し考えてみます
なるほど、スッキリした
① P(本物＆10連続表)＝(999/1000)*(1/1024)
② P(偽物＆10連続表)＝(1/1000)*(1)
①：②＝999：1024
P(本物|10連続表)＝①/(①+②)＝999/2023
P(偽物|10連続表)＝②/(①+②)＝1024/2023
1025枚以上から結論が逆になるんだな
故に (y ')^2 = c - 2/y,
ですね。cは積分定数です。
そこから先は xをyの関数と見て
x = ∫(1/y ') dy
= ∫√{y/(cy-2)} dy
= (1/c)√{y(cy-2)} + {2/c^(3/2)}･log[c√y + √{c(cy-2)} ] + c'
のような式になり、逆関数を求めるのは難しい。
・用途
クーロン散乱・ラザフォード散乱で正面衝突する場合(θ=0)とかに使えるかなぁ。
y'' * y' = 1/ y' とか y'' * y'^2 = 1
と入力したら答えでるね
それによると
y'(x) = v(x) とおけば
v' * v^2 = 1
積分して
v^3 / 3 = x + c
以下略
コインを1000回投げて10回以上連続して表がでる確率は？
1000回投げて10回以上連続して表
[1] 0.3854498
1000回投げてちょうど10回連続して表
[1] 0.1700181
1000回投げて20回以上連続して表
[1] 0.000468149
1000回投げてちょうど20回連続して表
[1] 0.0002342865
http://tpcg.io/TXVPdJ
これを参考にプログラム
https://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q1149349046
確率誤答の達人が全角文字で組んでくれるかもｗｗｗ
A[n]：○と×併せてn個を一列に並べた文字列で、○の１０連以上が無いものの数。
B[n]：○と×併せてn個を一列に並べた文字列で、○の１０連以上を含むものの数。A[n]+B[n]=2^n
A[n;k]：A[n]の中で、最後に○がk個連続しているもの。A[n]=A[n;0]+...+A[n;9]
P[n]：A[n;k](k=0～9)とB[n]を並べた、11成分の縦ベクトル
P[1]={1,1,0,0,0,0,0,0,0,0,0}^t
X={{1,1,1,1,1,1,1,1,1,1,0},
{1,0,0,0,0,0,0,0,0,0,0},...,
{0,0,0,0,0,0,0,0,1,0,0},
{0,0,0,0,0,0,0,0,0,1,2}}
P[1000]=X^(999) P[1] の第11成分を2^1000で割ったものが、求める確率
41301272734778977984946818232089531229879543376756574850136155867680807079676964
05909423852137579237591446526939613263507523948827986893531646240157193872907615
64116695521478307244714549348159061083607249922721310512099499789154886902065157
8128373092635280064104398841562373328900104830268510093352961/2^1000
≒0.38544975241248163591...
y = - (9/2)^(1/3)･|x-c '|^(2/3)
と解ける。
合計11枚が机の上に裏向きに置いてある。
ランダムに1枚ずつ引いていった場合の、得られた数字の総和の期待値を求めよ。
ただし、ジョーカーを引いた時点で終了するものとし、
Aは数字扱いではなく、最終的に得られた数字の総和が２倍になるものとする。
△ABK=△BCLとなる点Kの位置を求めよ。
シミュレーションしたら４５くらいになった。
> summary(re)
Min. 1st Qu. Median Mean 3rd Qu. Max.
41.80 44.29 45.01 45.03 45.81 48.46
１００万回のシミュレーションで総和の分布はこんな感じになった。
http://i.imgur.com/twslLtg.png
33になる？
2~10のうちジョーカーより左にくるものの期待値は44/2=22。
確率1/2で2倍になるから33。
シミュ合ってる？
54/2の3/2倍で41だ。
ジョーカーが早めに出る場合はAで2倍になる確率が低い
Aが出るまでの値だけを２倍？
袋の中から球を1個取り出すことをn回繰り返すとき、同じ色がk回連続して取り出され、かつk+1回以上は連続して取り出されない確率をP(n,k)とおく。
以下の問に答えよ。
（1）P(n,k)/P(n+1,k+1)をnとkの式で表せ。
（2）lim[n→∞] Σ[k=1 to n] kP(n,k) を求めよ。
自分は媒介変数表示のまま、√（x'^2 +y'^2）を積分して回答しました。
https://i.imgur.com/LU5BKnS.jpg
そう、Aを引いたら全部2倍
f <- function(x){
i=1
y=numeric()
while(x[i]!=11){ #11:ジョーカーでなければ
y[i]=x[i] # yに保存
i=i+1
}
if(1 %in% y) return(2*(sum(y)-1)) # 1があれば総和から1引いて2倍
else return(sum(y))
}
# simulation ,sample関数で1から11をランダムに並べ替え変え
re=replicate(10^6,f(sample(11))) #100万回fを繰り返す
summary(re)
hist(re,col='lightblue',xlab='sum',main='')
詳しい解説ありがとうございます
微分方程式のサイトを調べてる感じだと後は何とかなりそうな気がします
全部ってのがJoker引くまでの全部なのかAをひくまでの全部なのかを聞かれてるんだよ。
23A47J
の場合
(2+3)×2 + 4 + 7
なのか
(2+3+4+7)×2
なのかどっちにも取れるんだよ。
自分が書いてる文章でホントにちゃんと意味が伝わるのか考えながら書かないと。
もとの問題文に明確に書かれているし、
>>476の質問者も単に「全部」と「Aが出るまでの値だけ」の2択で質問しているから紛らわしさはない
図形的な解釈はハイレベル理系数学などに出てるから本屋で見ろ
数式で示すなら
x = rcosθ
y = rsinθ
にパラメータ表示の公式を適用して整理すればいい（ r はθの関数であることに注意）
23A47Jの場合、A=1 J=11
>480のfで
> f(c(2,3,1,4,7,11))
[1] 32
でいいんだよね？
それでよければシミュレーションでいいと思うんだが。
言いたかったことを代弁してくれて、ありがとう
ありがとうございます
媒介変数表示の時の図形的意味はなんとなく分かっているので
これで一応理解できたと思います
∫√(r^2+r'^2)から直接図形的意味を説明できませんでしょうか？
とても知りたいです
そういう変形を使わず余弦定理とかを使って証明する方法があれば知りたいということです
> トランプのA～10の10枚とジョーカー1枚の
> 合計11枚が机の上に裏向きに置いてある。
> ランダムに1枚ずつ引いていった場合の、得られた数字の総和の期待値を求めよ。
> ただし、ジョーカーを引いた時点で終了するものとし、
> Aは数字扱いではなく、最終的に得られた数字の総和が２倍になるものとする。
この文章のどこにエースを引くまでの全部と読み取れる要素があるねん？
むしろ最終的に得られた数の二倍ってJ引くまでのトータルを二倍としか読めない気がするけど。
コレがエースを引くまでのみを二倍の意味ならどう考えても言葉の選択間違ってるやろ？
１１個の順列を列挙して>480のｆを適用して平均値をだそうと思ったのだが、
11! = 39916800 なので PC処理が終わらないので断念。
なかなかやるじゃないですか(´・ω・`)
fの逆関数をgとすると
g(x) = x^(1/3)
となって、導関数は
g'(x) = (1/3) x^(-2/3)
となります。するとx=0でg'(x)は無限大に発散して連続でないように思えます。
g(x)は1階連続微分可能ではないということでよいのでしょうか？
もし1階連続微分可能でないとすると、f(x)とg(x)のグラフは回転・反転
させただけで滑らかさは全く変わらないことと不整合なように思えるのですが
どのように考えればよいのでしょうか？
∫√(r^2+(dr/dθ)^2) dθ で、dθが正になるように積分範囲を決定すると
= ∫√((rdθ)^2 + (dr)^2)
だから、三平方の定理で、直角を挟む辺の長さが(rdθ)、(dr)の直角三角形を
考えると斜辺の長さは √((rdθ)^2 + (dr)^2)
この斜辺の長さを足しあわせたものが曲線の長さになる。
非常に物理的な大雑把で直観的な説明。
ハイレベル理系数学持ってないけど、多分似たような説明だと思うので絵だけでも本屋で見てください。
袋の中から球を1個取り出すことをn回繰り返すとき、同じ色がk回連続して取り出され、かつk+1回以上は連続して取り出されない確率をP(n,k)とおく。
以下の問に答えよ。
（1）P(n,k)/P(n+1,k+1)をnとkの式で表せ。
（2）lim[n→∞] Σ[k=1 to n] kP(n,k) を求めよ。
×確率1/2で2倍になるから、54/2の3/2倍で40.5 (54の3/4倍で40.5)
○確率2/3で2倍になるから、54/2の5/3倍で45 (54の5/6倍で45)
任意の数字をＮとすると、Ｎ,Ａ,Ｊを引く順番と期待値は
Ｎ→Ａ→Ｊ (1/6)*(2Ｎ)
Ｎ→Ｊ (1/6)*Ｎ
Ａ→Ｎ→Ｊ (1/6)*(2Ｎ)
Ａ→Ｊ (1/6)*0
Ｊ (2/6)*0
合計すると (5/6)*Ｎ
求める期待値は (5/6)*(2+3+4+5+6+7+8+9+10)＝(5/6)*54＝45
なる！
(2)
Q(n,k) を少なくとも一回同じ色がk回連続引く確率とする。
m = [log √n]、l = [n/m] とおく。1～n回のコイントスのなかから連続する m 回のコイントスを重複しないように l 回に分けることが出来る。
T[1]～T[l] をそのような m 回のコイントスとしT[i]がすべて表になる事象をX[i]とするとき
P(X[i])=2^(-[log√n]) > 2^(-log√n) > 1/√n
である。
すべての i でX[i]が起こらない事象をYとするとき
P(Y) < (1-1/√n)^l < (1-1/√n)^(n/m) < exp (n/m) log(1-1/√n) < exp(-(√n/m))。
よって
P(n 回中 [log √n] 回連続表がでる)>(1-exp(-(√n/m)))。
∴Σ[k=1 to n] kP(n,k)
＞Σ[k=[log √n] to n] [log √n]P(n,k)
＞[log √n](1-exp(-(√n/m)))→∞。
a_1 = 1
a_(3n+1) = a_(2n+1)
a_(3n-1) = a_(2n-1)
a_(3n) = -a_n,
を満たす。この時、 lim(n→∞) (1/n)Σ(k=1，n) a_k を求めよ。
(面白スレ２８ より)
このスレも残り半分になりました。
考え方を教えてください
https://i.imgur.com/wq2ieeN.png
Prelude> [(a4,a3,a2,a1,x2,x1)|a4<-[1..9],a3<-[0..9],a2<-[0..9],a1<-[0..9],x2<-[1..9],x1<-[0..9],
(a4*1000+a3*100+a2*10+a1)*(x2*10+x1)==119868,(a4*1000+a3*100)*(x2*10+x1)==117600]
[(1,4,2,7,8,4),(2,8,5,4,4,2),(4,2,8,1,2,8)]
a=4281 x=28
aの下二桁とxを掛けると2268
aの上二桁とxを掛けると1176
それぞれ素因数分解する
それらを見比べるとaの下二桁、上二桁にそれぞれ必ず含まれる因数がわかる
そこに残りの因数をどれだけ移せるかを考える
4281と28で合ってる？
差を取って素因数分解すれば2*2*3*3*3*3*7になるので
あとは自分で考える
https://i.imgur.com/nQfh0oK.jpg
x = ∫(1/y ') dy
= ∫√{y/(cy-2)} dy
= (1/c)√{y(cy-2)} + {2/c^(3/2)}･log[c√y + √{c(cy-2)} ] + c'
ここの2行目から3行目の変形わかりませんでした
出来るもんだと思い込んでただけで全くわかりません
どうやって導き出したんでしょうか？
リーマンの写像定理から自明でした
んなわけない。
開集合じゃなきゃどうすんの？
{x^2+y^2≦1} ∪ [1,2] ∪ {(x-4)^2 + y^2≦4}
とか。
いくらでも複雑な例作れるよ？
んな簡単なわけない。
(A)a^bは自然数
(B)任意の有理数pに対して、a^pは無理数
申し訳ありません
途中から開集合で考えていたので整合性がとれていませんでした
開集合でなくてもいいならトポロジストの正弦曲線からも作れたりしますよね
a=e、b=log2
2π*(3+1)/2π*(1)＝4回転
限界などない
無限にいける
https://i.imgur.com/koStKCg.png
y = (2/c) cosh(t)^2,
cy - 2 = 2 sinh(t)^2,
dy = (4/c) sinh(t) cosh(t) dt,
とおいてみる。
∧__∧
( ´･ω･)∧∧l||l
/⌒ ,つ⌒ヽ ) <>>353
(___ ( __)
"''"" "'ﾞ''` 'ﾞ ﾞﾟ' ''' '' ''' ﾟ`
ありがとうございます
coshとか使うの初めて
でも、固定円の円周をちょんぎって直線にしたら1回転になる？
不思議と言えば不思議だな
とりあえずxi≧0に限定してti = √xiとおいて
f(x) = Σti^(p/2)、Σti = 1
なので凸不等式使えばいいのでは？
訂正 ti = xi^2です。
ぐるぐるした線なら円以上にたくさん回転するよ
こういうの
https://i.imgur.com/tyuKM7x.jpg
各成分の正負が一致してる場合はその方法やラグランジュなりで解けるのですが
一致しない場合、各成分で正負が違う場合の求め方が何とも…
1/6+2/6=1/2の確率でNを引かずにおわるかなぁ？
成分が負の場合を扱うなら 負^p を扱う事になるのでそこをどうするか決めないと答え出ない希ガス。
pが整数の時しか考えないとか。
そのＮは9枚じゃなくて1枚ね。
最後に9枚分の期待値を単純合計してもいいという理屈は
自分でもあまり良く分かってない。(問題>>467)
いきなり終了のＪが1/3もあるから、そんなに直感に反しないと思うんだけど
i:2～10に対し確率変数X[i]を
X[i] = 2i (i A J)
i (i J A)
2i (A i J)
0 (A J i)
0 (J i A)
0 (J A i)
とおく。
E(X[i]) = 5/6iである。
よって
E(得点) = Σ E(Xi) = Σ 5/6i = 5/6(2+3+…+10) = 5/6×54 = 45。
x = ∫(1/y ') dy
= ∫√{y/(cy-2)} dy
y = (2/c) cosh(t)^2,
cy - 2 = 2 sinh(t)^2,
dy = (4/c) sinh(t) cosh(t) dt
とおくと
x = ∫{4/c^(3/2)}･cosh(t)^2 dt
= ∫{2/c^(3/2)}･(cosh(2t) + 1) dt
= {1/c^(3/2)}･sinh(2t) + {2/c^(3/2)}t + c'
ここで
y = (2/c) cosh(t)^2
= (2/c)･{(e^(2t) + e^(-2t) + 2)}/4
e^(2t) = s とすると
y = (2/c)･{(s + 1/s + 2)}/4
ｓでそろえると
s^2 - 2(cy-1)s + 1 = 0
{s - ((cy-1))}^2 = 0
s = cy-1
e^(2t) = cy-1
両辺にlogすると
t = (1/2)log(cy-1)
ゆえに
x = (1/c)√{y(cy-2)} + {1/c^(3/2)}･log(cy-1) + c'
となりましたが答えが合いませんでした
どこで間違えたのでしょうか？
(B)はどうやって証明する？
[,1] [,2] [,3]
[1,] A J N
[2,] A N J
[3,] J A N
[4,] J N A
[5,] N A J
[6,] N J A
この各行が同様に確からしく起こるってことでいいんだな。
OM ≧ OS が成り立つのはなぜでしょうか？
https://ja.wikipedia.org/wiki/%E3%83%AA%E3%83%B3%E3%83%87%E3%83%9E%E3%83%B3%E3%81%AE%E5%AE%9A%E7%90%86
s^2 - 2(cy-1)s + 1 = 0
まで正しいが、次から違っている。
{s - (cy-1)}^2 = cy(cy-2),
∴ s = (cy-1) ± √{cy(cy-2)},
e^t = √s = [ √(cy) ± √(cy-2) ] / √2,
t = log[ √(cy) ± √(cy-2)] - (1/2)log(2),
100円玉２個でやってみた
助かった ありがとう
OM ≧ OS が成り立つのはなぜでしょうか？
https://imgur.com/OX9wyZE.jpg
Mのx座標≧Sのx座標、Mのy座標≧Sのy座標だから。
そんなもん有難がってもいい事ないぞ
・a^2≦0かつa^2>0
虚実数は、実数kと虚数単位iでは表せないことを示せ。
ありがとうございました。
動点Pについて、f(P)=PA^2+PB^2+PC^2とする。
（1）kをk≧√3である実数とする。動点Pが平面z=0上を動くとき、f(P)=kとなる点Pの軌跡で囲まれる面積を求めよ。
（2）動点Pが空間を動くとき、f(P)=kとなる点Pの軌跡で囲まれる面積を求めよ。f(P)=kとなる点Pの軌跡で囲まれる体積を求めよ。
（2）は鮮やかに解いてください
積分してもいいですが
タイヤが道路と設置した長さだけ自動車は移動していると言われれば納得できるな。
設置→接地
|p-a|^2+ |p-b|^2+ |p-c|^2
=3|p-g|^2-2(|g-a|^2+ |g-b|^2+ |g-c|^2)
考え方を教えてください
https://i.imgur.com/wq2ieeN.png
上で回答ついてるじゃんか
差をとると2268、これはaの下二桁とｘを掛けた値。
素因数分解すれば2*2*3*3*3*3*7になる。
aの下二桁をｙとすると一方が22以下なら他方が2桁にならないから23以上。
2268の平方根は47.62352なのでx,yのいずれかは47以下である。
2 2 3 3 3 3 7の積がこの範囲にあるのは27 28 36 42の4つ
そのときの他方の数は84 81 63 54
あとは自分で考える。
組み合わせを考えるのが面倒だから
2268の約数で23以上、47以下はプログラムにやらせた。
> x=2268/(23:47)
> y=(23:47)[x-floor(x)==0]
> z=2268/y
> rbind(y,z)
[,1] [,2] [,3] [,4]
y 27 28 36 42
z 84 81 63 54
ついでだから続きも書いておく。
xの候補が27 28 36 42 84 81 63 54に絞られたので
119868を割り切るのは
28 42 84の３つ
そのときの商は 4281 2854 1427
でこれがaの候補。
最大は4281でそのときのｘは28
> x=2268/(23:47)
> y=(23:47)[x-floor(x)==0]
> z=2268/y
> b=c(y,z)
> c=119868/b
> d=b[c-floor(c)==0]
> (a=max(119868/d))
[1] 4281
> 119868/a
[1] 28
>
Haskellだと1行ですんだ。
Prelude> [(a,x)|a4<-[1..9],a3<-[0..9],a2<-[0..9],a1<-[0..9],x2<-[1..9],x1<-[0..9], let a=(a4*1000+a3*100+a2*10+a1),let x= x2*10+x1,a*x==119868,(a4*1000+a3*100)*(x2*10+x1)==117600]
[(1427,84),(2854,42),(4281,28)]
Prelude> [(a,x)|b<-[10..99],c<-[0..99],x<-[10..99],let a=100*b+c, a*x==119868,100*b*x==117600]
[(1427,84),(2854,42),(4281,28)]
ありがとうございます
式を変形して綺麗にしてあるのもようやく理解出来ました
わからないところは>>460の x = (1/c)√{y(cy-2)} + {2/c^(3/2)}･log[c√y + √{c(cy-2)} ] + c'
この式の log[c√y + √{c(cy-2)} ] の部分で+の部分が±ではない理由は何でしょうか？
>>504にあるけど1176のほうも使った方が簡単じゃね？
の組み合わせのうち（3,4,5）から数えて7番目は何になるかという問題がわかりません
なるほどね
1176=2*2*2*3*7*7
2268=2*2*3*3*3*3*7
でxは公約数か
Prelude> ps = [(a,b,c)|a<-[1..100],b<-[a..100],c<-[b..100],a^2+b^2==c^2]
Prelude> ps !! 7
(10,24,26)
顰蹙のダンプリスト
Prelude> [(a,b,c)|a<-[1..100],b<-[a..100],c<-[b..100],a^2+b^2==c^2]
[(3,4,5),(5,12,13),(6,8,10),(7,24,25),(8,15,17),(9,12,15),(9,40,41),(10,24,26),(11,60,61),(12,16,20),(12,35,37),(13,84,85),(14,48,50),(15,20,25),
(15,36,39),(16,30,34),(16,63,65),(18,24,30),(18,80,82),(20,21,29),(20,48,52),(21,28,35),(21,72,75),(24,32,40),(24,45,51),(24,70,74),(25,60,65),
(27,36,45),(28,45,53),(28,96,100),(30,40,50),(30,72,78),(32,60,68),(33,44,55),(33,56,65),(35,84,91),(36,48,60),(36,77,85),(39,52,65),(39,80,89),
(40,42,58),(40,75,85),(42,56,70),(45,60,75),(48,55,73),(48,64,80),(51,68,85),(54,72,90),(57,76,95),(60,63,87),(60,80,100),(65,72,97)]
Prelude> ps !! (7-1)
(9,40,41)
７番目だと変わらないみたいだけど
整数の候補の上限を変えたら答が変わるな。
pitNth n m = do
let ps = [(a,b,c)|a<-[1..m],b<-[a..m],c<-[b..m],a^2+b^2==c^2]
map (\x -> ps !! x) [0..(n-1)]
２桁の９９までだと２０番目は,(20,21,29)
Prelude> pitNth 20 99
[(3,4,5),(5,12,13),(6,8,10),(7,24,25),(8,15,17),(9,12,15),(9,40,41),(10,24,26),(11,60,61),(12,16,20),(12,35,37),(13,84,85),(14,48,50),(15,20,25),(15,36,39)
,(16,30,34),(16,63,65),(18,24,30),(18,80,82),(20,21,29)]
３桁の９９９までだと２０番目は,(18,24,30)
Prelude> pitNth 20 999
[(3,4,5),(5,12,13),(6,8,10),(7,24,25),(8,15,17),(9,12,15),(9,40,41),(10,24,26),(11,60,61),(12,16,20),(12,35,37),(13,84,85),(14,48,50),(15,20,25),(15,36,39)
,(15,112,113),(16,30,34),(16,63,65),(17,144,145),(18,24,30)]
> 整数の候補の上限を変えたら答が変わるな。
二秒で分かりそうなもんだけどww
すると７番目が(9,40,41)というのはどうやって確信できるんだろう？
AはBに10秒で10mの差をつける
BはCに10秒で10mの差をつける
CはDに10秒で10mの差をつける
DはEに10秒で10mの差をつける
AがEに10ｍの差をつけるのは何秒後？
ABCDEの5人が一緒に直線道を同時にスタート
10秒後にAはBに10mの差をつけた
20秒後にBはCに10mの差をつけた
30秒後にCはDに10mの差をつけた
40秒後にDはEに10mの差をつけた
AがEに10ｍの差をつけたのは何秒後？
計算力とかプログラムの力ではなくてエレガントな解法を知りたい
[A]と[f(A)]はそれぞれAとf(A)の閉包を表しています。
ＮＧに紛れ込んで見えてなかったらすまんとしか言いようないけど
だって計算機つかえば一瞬で答え出るような問題頭使う気しない。
a^2 + b^2 = c^2 より
c=b+n とすれば a^2 = 2bn+n^2
nが自然数なら b の最大値は n=1 のとき (a^2-1)/2
a が√201 ≒ 14.17 を超えない限り2桁のbまで調べれば十分
既約ピタゴラス数にも通し番号付ければいいのにね。
ケッヒェル番号、ドイチュ番号、ホーボーケン番号、などなど。
xyz空間に、原点Ｏを中心とする半径１の球体Ｇがある。
また、放物線 y = 1-x^2, z=0 をz軸方向に平行移動して得られる曲面によってxyz空間を2つに分解したとき、原点Ｏを含まない方をＴとする。
ＧとＴの共通部分Ｇ∩Ｔの体積を求めよ。
「GUT は善だ」（ドイツ人）
「GUT は腸だ」（英米人）
y=b で切ったときの断面を考える。(0≦b≦1)
-√(1-bb) ≦ x ≦ -√(1-b) および √(1-b) ≦ x ≦ √(1-bb),
-√(1-bb-xx) ≦ z ≦ √(1-bb-xx),
∫ dz = 2√(1-bb-xx),
S(b) = 2∫[√(1-b), √(1-bb)] 2√(1-bb-xx) dx
= 2 [ x√(1-bb-xx) - (1-bb) arctan{√(1-bb-xx) /x} ]
= 2{ -(1-b)√b + (1-bb) arctan(√b) },
V = ∫[0,1] S(y) dy
= 2{-4/15 + (π/3 -32/45)}
= 2(π/3 -44/45),
V = ∫[0,1] S(y) dy
= 2∫[0,1] { (1-yy)arctan(√y) - (1-y)√y } dy
= 2 [ (1/3)(2-y)(1+y)^2 arctan(√y) - (1/45)(30+35y-21yy)√y ](0,1)
= 2 (π/3 - 44/45)
= 0.13883954683764
ずっとZで切った断面積考えててわからなかったわ、わざわざすみません。
解けるでしょ
5人とも一定の速度で走るとすれば
v(A) - v(B) = 1 (m/s)
v(B) - v(C) = 1/2 (m/s)
v(C) - v(D) = 1/3 (m/s)
v(D) - v(E) = 1/4 (m/s)
辺々たすと
v(A) - v(E) = 1 + 1/2 + 1/3 + 1/4 = 25/12 (m/s),
よって
10 / (25/12) = 4.8 (s)
放物線と円の共通してる面積がうまく表せないんや。教えてください。
ちょっと解いてみてw
全く同じとまでは言わないけど
かなり類似だろ。
鳩ノ巣原理が一対一の全単射関係の濃度なら
ハッシュテーブルは箱と中身で同値類と代表元なんだから。
ボロノイ図やゲージ固定も類似だね。
ありがとうございます。
お礼にaの上限を３０にして算出してみました。
Prelude> m=30
Prelude> [(a,b,c)|a<-[1..m],b<-[a..floor(m^2/2-1/2)],c<-[b..2*b],a^2+b^2==c^2]
[(3,4,5),(5,12,13),(6,8,10),(7,24,25),(8,15,17),(9,12,15),(9,40,41),(10,24,26),(11,60,61),(12,16,20),(12,35,37),(13,84,85),(14,48,50),
(15,20,25),(15,36,39),(15,112,113),(16,30,34),(16,63,65),(17,144,145),(18,24,30),(18,80,82),(19,180,181),(20,21,29),(20,48,52),
(20,99,101),(21,28,35),(21,72,75),(21,220,221),(22,120,122),(23,264,265),(24,32,40),(24,45,51),(24,70,74),(24,143,145),(25,60,65)
,(25,312,313),(26,168,170),(27,36,45),(27,120,123),(27,364,365),(28,45,53),(28,96,100),(28,195,197),(29,420,421),(30,40,50),(30,72,78),(30,224,226)]
[(a,b,c)|a<-[1..m],b<-[a..floor(m^2/2-1/2)],c<-[b..2*b],a^2+b^2==c^2]
が遅いので速度を上げようとしたけど下記ではエラーが返ってきた。達人にデバックを期待(._.)
[(a,b,c)|a<-[1..m],b<-[a..floor(a^2/2-1/2)],c<-[b..floor(sqrt(a^2+b^2))],a^2+b^2==c^2]
[(a,b,c)|a<-[1..m],b<-[a..floor(m^2/2-1/2)],let c = sqrt(a^2+b^2), fromIntegral(floor(c))==c]
今回はダンプリストではなくてRのスクリプト（HaskellやPythonは独学中ｗ）
a^2+b^2の平方根が整数の組み合わせを考えればいいんだから、簡単にプログラムが組めた。
A=100
pita=NULL
for(a in 1:A){
B=floor(a^2/2-1/2)
for(b in a:B){
c=a^2+b^2
if(floor(sqrt(c)) == sqrt(c) ){
pita=rbind(re,c(a,b,sqrt(c)))
}
}
}
> pita[7,]
[1] 9 40 41
> pita[77,]
[1] 42 56 70
> pita[100,]
[1] 50 120 130
> pita[777,]
a b c
216 288 360
というか、計算機に答を出す命令を組むのが楽しいんだよね。
このあたりは価値観の問題だよね？
２の平方根の１００桁めの数字を出すのは不毛に思えるけど
１００個目のピタゴラス数を計算するのは不毛に思えない人がいるのがこのスレだと思っている。
原始ピタゴラス数を表示するプログラムと解説は、そのスレッドの119と120がお勧めです。
https://rio2016.5ch.net/test/read.cgi/math/1478040803/
レスありがとうございます。
おかげて次のステップのプログラムができるようになりました。
これが分からないんですが
たとえば仮に軸をy=x、θを45°とした場合
このような薄い立体の体積がなぜ、側面積*⊿xで求まるのかが分かりません
側面積*⊿√2xとならないのはなぜですか？
お願いします
ごめん風疹で頭いかれてる
スレ汚し失礼
M⊂RがコンパクトであることとMの最小値の存在が同値であることってどう示すんですか？
Mに最小値がないとする。
単調減少列x[i]∈Mをlim X[i] = -∞ or lim x[i] = inf M ととれる。
このとき M ⊂ ∪ (x[i],∞) であるが有限個ではM全体を被覆しない。
Mが最小限mをもつとする。
被覆 M ⊂ ∪U[i] に対し x∈U[i0] である i0 をとれば M ⊂ U[i0] である。
https://i.imgur.com/DHcSQPj.jpg
[命題: Mはコンパクトである ←→ Mに最小値が存在する]
(←) Mに最小値 α が存在する時
任意の無限開被覆 {(x_λ, +∞) ; λ ∈ Λ } に対して α ∈ (x_ξ, +∞) となる ξ ∈ Λ が存在する.
この時、 (x_ξ, +∞) ただ１つで 有限開被覆となる. よってコンパクトである.
(→) 対偶で示す. Mに最小値が存在しない時
M の下限 β をとる. β= -∞ なら、有限開被覆は常に不可能.
βが有限なら、Mの無限開被覆 {(β + 1/n, +∞) ; n=1,2, ... } から有限開被覆は取り出せない.
よってコンパクトではない.
x=a で切ったときの断面を考える。(-1≦a≦1)
-√(1-aa-yy) ≦ z ≦ √(1-aa-yy), … 円の内部
y ≧ 1-aa,
なので弓型である。
S~(a) = ∬ dz dy
= ∫[1-aa, √(1-aa)] 2√(1-aa-yy) dy
= [ (1-aa)arcsin{y/√(1-aa)} + y√(1-aa-yy) ](y=1-aa,√(1-aa))
= (1-aa)arcsin(a) - a(1-aa)^(3/2),
あるいは
S~(a) = ∬ dy dz
= 2∫[0, a√(1-aa)] {sqrt(1-aa-zz) -(1-aa)} dz
= [ (1-aa)arcsin{z/√(1-aa)} + z√(1-aa-zz) -2(1-aa)z ](z=…)
= (1-aa)arcsin(a) - a(1-aa)^(3/2),
V = 2∫[0,1] S~(x)dx
= 2∫[0,1] { (1-xx)arcsin(x) - x(1-xx)^(3/2) }dx
= 2[ (1/3)x(3-xx)arcsin(x) + (1/45)(9x^4 -23xx +44)√(1-xx) ](x=0,1)
= 2(π/3 - 44/45),
>>599
それは解けぬ...
M = (Jf)^{-1}|x=0 と置くと、
F[i] = M[i,k] { f[k] - .. } より
JF[i,j] = ∂F[i]/∂x[j] = M[i,k] ∂f[k]/∂x[j] = M[i,k] Jf[k,j] = (M. Jf)[i,j] = δ[i,j] (x=0)
F(0) = 0, C^∞ は明らか.
>>620
分かりやすい解答ありがとうございます
なんで解けない方ばかり行くんだろうねぇ
である。
（1）切り分けられた立体のうち、原点O(0,0,0)を含む方の体積Vをaで表せ。
（2）（1）の結果を用いて、次の定積分の値を求めよ。ここでg(x)はf(x)=sinxの逆関数であり、定義域は0<x<90°とする。
∫[0→sina°] g(x) dx
梅沢富美男ぢゃないんだから…（ローソンのCM)
z=c で切ったときの断面を考える。(-1≦c≦1)
三日月形(？)になる。
{1-√(1-4cc)}/2 ≦ x^2 ≦ {1+√(1-4cc)}/2,
x1 = √{[1-√(1-4cc)]/2},
x2 = √{[1+√(1-4cc)]/2},
とおくと
Sz (c) = 2∫[x1, x2] {√(1-cc-xx) - (1-xx)} dx
= [ (1-cc)arcsin(x/√(1-cc)) + x√(1-cc-xx) - 2{x - (1/3)x^3} ](x=x1,x2)
V = ∫[-1/2, 1/2] Sz(z) dz = …
かなり面倒だ…
(2)がわかりません。ちなみに私立の推薦なので答えは不明です。どなたかよろしくお願いします。
V = (2/3) tan a°
p = 1/2 のとき 0 でない値に収束。このとき a = √5
OR = t √( t^2 + a^2 t^{2p} ) / (√( t^2 + a^2 t^{2p} ) - a t^p )
p=1 の場合
OR = t √( 1 + a^2 ) / (√( 1 + a^2 ) - a ) → +0 ≠ 10 よって不可.
p＞1 の場合
OR = t √( 1 + a^2 t^{2p-2} ) / ( √( 1 + a^2 t^{2p-2} ) - a t^{p-1} ) → +0 ≠ 10 不可.
0＜p＜1 の場合
OR = a t^{p} √( a^{-2} t^{2-2p} + 1 ) / { a t^{p-1} (√( a^{-2} t^{2-2p} + 1) - 1) )
= ( t + (1/2) a^{-2} t^{3-2p} +... ) / { (1/2) a^{-2} t^{2-2p} + ... }
OR → 10 が可能となるのは、2a^2 = 10, 2-2p = 1 の時
すなわち、 a=√5 , p = 1/2
ABCの3つの箱から当たりのある箱を選ぶ
最初に選んだ箱をAとする
当たりが●、ハズレが○、？は●と○が不確定な状態
一つも開示されない状態の箱は
A○○？ B○？ C？
Cを開けると○
AとBの○が一個減るので
A○？ B？
となるから、Bが当たりになる確率が上がるって話？
君がモンティホール問題と呼んでいる問題の問題文を端折らずに書いてみてくれないか
期間限定公開の
数字のいたずらっていう動画でやってたやつ
https://youtu.be/7BNcQpDhV94?t=78
モンティホール問題は
「司会者は答えを知っていて、自分が開けるようなヘマはしない」
という構造の問題にすぎないので、何か深淵な数理的秘密があるのだろうと
思っていると訳が分らなくなります。
> A○○？ B○？ C？
> A○？ B？
の意味がわからない
変えた方が確率が高いのは、変えると箱を二つ選ぶのと同じことになるからだよ
二つ選んだ上でその中に当たりがあったら教えてもらえるというのと同じことだから
確率が上がるという説明がどうにも納得できない
最初に選んだ方が当たりであることっていうのが1/3しかないから
2/3に乗り換えられるなら乗り換えた方が確率は高くなる
上で書いたけど要するに
最初A選んだときに「ではAならその1個だけですけどそれをやめてBとCの2個選んでもいいですけどどうしますか？」って言われてるのと同じ
初め当たりを引いていた場合 (1/3)*0
初めハズレを引いていた場合 (2/3)*(1/2)
加えると 1/3
初めから変えなければもちろん 1/3
同じではないのか。
初めハズレを引いていた場合 (2/3)* 1
の間違い。理解した。
Aを選んでハズレのDが開けられて、ABCの3択にされて
Bに変えた後にハズレのCを開けられて、ABの2択にされた。
ABの当たり確率は？ (司会者は当たりを知っている)
もちろん1つだけ当たりの場合だけど
ヒントは条件付確率
(ハズレのCが開けられる確率は？)
(a) 当たりになる場合＝最初に選んだ箱が当たりである場合。
その確率は、
1/3
(b) はずれになる場合＝最初に選んだ箱がはずれである場合：
その確率は、
2/3
(2) 箱を必ず変更する方針の場合
(a) 当たりになる場合＝最初に選んだ箱がはずれである場合。
その確率は、
2/3
(b) はずれになる場合＝最初に選んだ箱が当たりである場合：
その確率は、
1/3
正解
つまりCDがハズレでプレイヤーがCDを選ばないことが予め決まっているなら
ABを選んだ場合の当たる確率はそれぞれ当初ABに当たりが入れられた確率
等確率ならどちらも1/2
そうではなくABCDに当たりが入れられた確率は等確率で、
司会者はプレイヤーの選択に応じて開けられるハズレを
選んで開けているだけ(今回たまたまDとC)なら
A 1/4 B 3/4 かな
最初に選んだ箱がハズレなら変えれば必ず当たるんじゃない？
>>630
遅くなりましたがありがとうございます！
モンティが確定情報をもとにハズレのドアを開けるから
プレイヤーが最初に選択したドアの確率は最後まで１／４のまま
こんなんネットで引いたら丁寧に解説してるサイト死ぬほど出てくるやん。
https://mathtrain.jp/monty
とか。
Aが当たりの確率は1/4です
BCDが当たりの確率はどれもそれぞれ1/4ですが、CDはハズレです
ABのどちらを選びますか？
↓ (選A、開E)
② 当(A,B,C,D)＝(1/5、4/15、4/15、4/15)
↓ (選B、開D)
③ 当(A,B,C)＝(9/29、8/29、12/29)
当(A,B,C)＝(1/4、3/8、3/8)
(選B、開C)
① P(当A ∧ 開C)＝(1/4)*(1)
② P(当B ∧ 開C)＝(3/8)*(1/2)
①：②＝4：3
P(当A｜開C)＝①/(①＋②)＝4/7
P(当B｜開C)＝②/(①＋②)＝3/7
ω=1/z=x+yiと置くとxyの一次式になってこういうゴリ押しでも解けたので、
zのままゴリ押しても行けるのではないかと思ったのですがどうしても分かりません。助けて下さい。
https://i.imgur.com/zDMY5iK.jpg
k=±1の場合は成立しないので除外しても構いません。
数3範囲外だったっけ・・・？
1/z=wとおいて|w|≧1/2が成立する範囲。
wの軌跡は
|w-(-i)| = |w-k|
だから(-i)とkの垂直２等分線。
条件満たすのは|k|≧1/√3のとき。
ありがとうございます。
1/zを置く方法では一応解けたので、zのままで虚部実部を分割して代数的に解く方法をお願いします。
z は原点を通る円になるので、直径≦2 になればよい
1枚目は最後の行でlog4-log2の直前に掛かる2はどこから来たのですか？
2枚目はマルで囲んだ引かれる積分の部分に(1+x)が掛かってないのはなぜですか？
https://i.imgur.com/J8uNVmX.jpg
https://i.imgur.com/MR5r7Uf.jpg
-{log(3-cos0)-log(3+cos0)}
= (log4-log2)
-(log2-log4)
cos0を0にしてました……
2枚目もどうしてもわからないのでお願いします……
2枚目は途中で切れてる
途中できれてると言いますと？
1+xは全体にかかるのでは？と思ったのですが
難問
(1+x)がどこかにいっちゃってるなあ
その先もなかったことで話が進んでるの？
やっぱどっかいっちゃってますかね。ありがとうございます。
すぐ下で件の項のn→無限での極限を求めてますがそのままです。
旺文社の大学入試全レベル問題集というやつでした。
1+xを掛けたままの状態で極限がゼロなことはどうやったらシンプルに証明できますか？
皆様のおかげで解決しました。ありがとうございます。
それぞれの扉が当たりの確率は10%です。
これを3グループ、5枚,3枚,2枚に分けます。
それぞれのグループをA,B,Cグループとします。
グループごとの当たりがある確率は50%,30%,20%です。
ここで、プレイヤーはAグループを選びます。
すると、モンティが残った2グループのうちの
Cグループの扉を開けてハズレだと教えました。
このとき、Aグループに当たりがある確率と
Bグループに当たりがある確率は同じでしょうか？
自然数nに対して有理数n/Nの循環節の長さをf(n)とおくとき、以下の各Nに対してf(n)を最大にするnを1つ求めよ。
（1）N=7
（2）N=17
（3）N=37
こいついい加減死なないかな
「(1)N=7」
痴呆症かな？
精度が大幅にアップグレード
P(A)／P(B)＝（P君の勝つ数）／（Q君の勝つ数）
｛｛n/(n+1)｝^n-1｝｛k(n-1)-n^3+n((n-1)/n)^n+n｝
P(A)／P(B)＝――――――――――――――――――――――
｛｛｛｛(n-1)/n｝^n-1｝n+1｝｛k-n^2+(n/(n+1))^n-n｝｝
∵［n≧2，n(n+1)-1＞k≧1］
Nとnを求めよ。
f(x)=f'(x)となるようなa（すなわちe）が存在することを中間値の定理を用いて示せ。
これら13個の点から相異なる2点を無作為に選んで結ぶとき、その線分の長さの期待値Eと1/2の大小を比較せよ。
aは負でもいいのか、そのときxの範囲として1/(奇数)とかを含めるかどうか
流石に実数値関数だろうからR全体ではないだろう
あとa=0のときの0^0はどうするのか、xの定義域に0を含めなければa=0でf(x)=f'(x)となるが
俺レベルでもひと目で自明とわかる問題をバンバン出してるし
中学生とかならいいけど、大学生以上でこれやってたら悲惨だなー
多分後者っぽい気がすんだよね
淡々と問題貼り続けるキチっぽさが
ご指導ありがとうございます
aは正の実数です
lim[t to 0] t^t = 1 と定義させていただきます
流石に(e^x)’=e^xを使うと自明になってしまうからダメだろなとまでは思うけど、じゃ(a^x)’=a^x log aはいいのかという話になる。
でも高校の教科書の定義は底がeの時の対数関数だからやっぱりダメっぽい。
するとそもそも論としてa^xの微分可能性は使っていいのかもかなり怪しくなる。
この問題何は仮定してよくて何は証明しないといけないのかがそもそもサッパリ。
此れは解答できるでせう
わざわざ素因数という言葉があるくらいだから高校数学の用語としては1は因数ということになると思う。
すると条件は
Nとnは互いに素、N-n=16と言ってるのと同じでこんなもん死ぬほど解ある。
・a[1] = c (1/3 < c ≦ 2/3)
・0 < a[3n] ≦ 1/3、1/3 < a[3n+1] ≦ 2/3、2/3 ≦ a[3n+2] < 1
・lim[n→∞] a[n] は収束する。
このとき、L = lim[n→∞] a[n] の取りうる値、もしくはその範囲を求めよ。
傑作だと思う。
カオス理論から帰着した
失礼。
lim[n→∞] a[3n]
だわ
失礼。
lim[n→∞] a[3n]
だわ
ln(x) < x/e より
(3/2)ln(2) = ln(2√2) < (2√2)/e,
ln(2) < (4√2)/(3e) < 1/√2 = 0.7071… (← e > 8/3)
7^6 = 117649 > 10^5,
7 > 10^(5/6),
log(7) > 5/6 = 0.8333…
∴ ln(2) < 1/√2 < 5/6 < log(7),
正n角形のとき
線分（辺または対角線）の長さは
2sin(kπ/n) (k=1,2,…,n-1)
確率はいずれも等しく 1/(n-1),
E_n = {1/(n-1)}Σ[k=1,n-1] 2sin(kπ/n)
= {2/(n-1)}cot(π/2n)
～4n/{(n-1)π}
→ 4/π (n→∞)
√3 (n=3) からnとともに減少する。
積和公式
2sin(kθ) = {cos((k-1/2)θ) - cos((k+1/2)θ)}/sin(θ/2),
から
Σ[k=1,n-1] 2sin(kθ) = {cos(θ/2) - cos(nθ -θ/2)}/sin(θ/2),
ここで θ=π/n とおけば、nθ=π より
= 2cot(θ/2)
実数a>0 と xに対して
f(x) = a^x
が定義されているとする。
f'(x) = lim(h→0) {f(x+h) - f(x)}/h
= lim(h→0) {a^(x+h) - a^x}/h
= (a^x) lim(h→0) (a^h - 1)/h
= (a^x) g(a),
とおく。
g(a) は連続函数で
g(1) = 0
a>1 のとき
g(a^m) = lim(h→0) {a^(mh) - 1}/h = m･lim(H→0) (a^H -1)/H = m g(a),
a がm乗になると、g(a) はm倍になる。
アルキメデスの原理により、これはいくらでも大きくなる。
中間値の定理より
g(e) = 1
を満足する e>1 が存在する。
N/d = x, n/d = y とおくと
x^d - y^d = 16/(d^d),
∴ d^d は 16 を割り切る。
∴ d=1，2
d=1 のときは >>690
d=2 のとき
x^2 - y^2 = 4,
(x，y) = (±2，0)
(N，n) = (±4，0) となる。(不適)
この関数の逆関数を求めるにはどうすればいいですか？
y = log(x) + x^2 = log(x) + log(e^{x^2}) = log(x *e^{x^2} )
e^y = x * e^{x^2}
2e^{2y} = 2x^2 *e^{2x^2}
W( 2e^{2y} ) = 2 x^2
∴ x = √( W( 2e^{2y} )/2 )
※ W(x)は ランベルトのW関数. f(x) = x e^x の逆関数として定義される.
lim[t→∞] L(t)/{f(t+1)-f(t)} を求めよ。
2a1+1=an+3
お願いします
f(x) = x^3 -3AAx -B,
とおくと
f '(-A) = f '(A) = 0
さらに
A = 1 +3t +3tt,
B = ±(-2A+1)(A+3t+1)(A-3t-2)
とおけば
f(-2A+1) = f(A+3t+1) = f(A-3t-2) = 0, or
f(2A-1) = f(-A-3t-1) = f(-A+3t+2) = 0,
正弦定理より a + b + c = 2 (sinA + sinB + sinC)
拘束条件は A + B + C = π
ラグランジュ未定乗数を μ として
F(A,B,C) = 2 (sinA + sinB + sinC) - μ*( A + B + C )
∂F/∂A = 2cosA - μ = 0, ... , ...
A = B = C = arccos(μ/2) = π/3 以下略
Mp を G = (V, E) の完全マッチングとする。
明らかに、
|Mp| = |V| / 2
が成り立つ。
よって、すべての完全マッチングの辺の数は等しく、その数は、 |V| / 2 である。
Mmax を G = (V, E) の最大マッチングとし、 |Mmax| > |Mp| と仮定して矛盾を導く。
明らかに、
2 * |Mmax| ≦ |V|
が成り立つ。
∴ |Mmax| ≦ |V| / 2 = |Mp|
これは矛盾である。
よって、完全マッチングは最大マッチングである。
訂正します：
G = (V, E) を完全マッチングをもつグラフとする。
Mp を G = (V, E) の完全マッチングとする。
明らかに、
|Mp| = |V| / 2
が成り立つ。
よって、すべての完全マッチングの辺の数は等しく、その数は、 |V| / 2 である。
Mmax を G = (V, E) の最大マッチングとし、 |Mmax| > |Mp| と仮定して矛盾を導く。
明らかに、
2 * |Mmax| ≦ |V|
が成り立つ。
∴ |Mmax| ≦ |V| / 2 = |Mp|
これは矛盾である。
よって、完全マッチングは最大マッチングである。
指定の値域で導関数がゼロになるものが一つしかない場合
論述でこれ書いたら雑な方法認定されて減点されますかね？
導関数が値域のどこかで正か負の無限大にならない場合、
+0+、-0-、+0-、-0+ の4パターンしかないですから端点と0の点だけ調べればいけますよね？
やっぱ増減表書かないとまずいでしょうか？
https://i.imgur.com/2tESt3F.jpg
答えは10%なんですが、過程式がわからないです…
よろしくお願いします
水と食塩を別々に考える
誰も解けないかんじですか？難問ですが解ける人いたらお願いします。
e^{2πia/b} 以下略
このとき、B[a]がA加群として有限生成なら、生成系はあるnが存在し{1,a,......,a^n}と取れることの証明を教えて下さい。
これはおかしいですね
AをB代数、b∈Bとして、A[b]がA加群として有限生成なら、生成系として{1,b,....b^n}がとれるでお願いします
BはA代数です
全部の解答解説お願いします。
(1) 整数 k, k’ が exp(i2π k/b) = exp(i2π k’ /b) となる必要十分条件は
2π k/b = 2π k’ /b + 2π n (nは適当な整数) の関係にある事である.
すなわち k ≡ k’ (mod b) であり、k=0, 1, ..., b-1 が相異なる exp(...) を与える. よって #P = b.
(2) a k ≡ 1 (mod b) を与える k が存在する. (∵ a, b は互いに素). (1)よりそれが求めたかった k である.
(3) 明らかに Q ⊂ P である. また(2)よりQは P の生成元を含む、よって P ⊂ Q.
(4) (3)より a1 = a2 = 1 としても同じ事である.
適当な k,k’ を選べば k/b1 + k’/b2 = k(b2 k + b1 k’)/(b1 b2) = 1/(b1 b2) とできる. (∵例えばユークリッド互除法)
よって (1)～(3)より #(Q1Q2) = b1 b2
A, Ab, Ab+Ab^2, ...... はそれぞれ有限生成でA[b]も有限生成だから、あるnが存在してA+Ab+...Ab^n = A[b]となる
よって1,b, ...... , b^nがとれる
とネーター加群の真似をしてみたのですが、これは正しいでしょうか？
15ｇだ。
残り500ｇの中に食塩は何ｇある？
75ｇだ。
250ｇ足したら食塩水は何ｇになった？
750ｇだ。
750ｇの食塩水の中に75ｇの食塩がある。何％だ？
10％だ。
式か？
100×0.15=15
600-100=500
15×(500/100)=75
500+250=750
(75/750)×100=10
この五式で満点だろう。
一般にMが有限生成、(m[i])がMの元の集合でM = Σ[i∈I]m[i]AとするとIの有限部分集合FがとれてM = Σm[i∈F]A。
(∵) M = Σ[j=1～n]n[j]Aとする。
各 j に対し有限集合 F[j] と a[i,j]∈Aで
n[j] = Σ[i∈F[j]]m[i]a[ij]
となるものがとれる。
F = ∪ F[j] とすれば n[j] ∈ Σ[i∈F] m[i]Aであるから
M ⊂ Σ[j=1～n]n[j]A ⊂ Σ[i∈F] m[i]A である。
問題 15%の食塩水600mLから100mLを使い、その後、水を250ｇ入れると[ ]%の食塩水になる。
と改変すると比重を考える必要が出てきて難問化するね。
有限生成であることの同値な言いかえとしてそのようなことが成り立つのは知りませんでした
ありがとうございます
必要であれば以下の事実を用いて良い。
「a^b-c^d=1を満たす2以上の自然数a,b,c,dはただ一組しか存在しない」
m≧2、n≧2
ありがとうございます
問題①
ある仕事を仕上げるのにAは10日、Bは15日、Cは20日かかる。この仕事を3人で協力して行ったが、途中Aが休んだので、仕上げるのに6日かかった。Aは何日仕事を休んだか？
答え 3日ですが、 過去問に過程式がないので理解できないです。
問題②
2人が自転車に乗って走っている。その速さの比は11:8である。この2人が周囲480mの円形の池の同じ地点から同時に同方向にスタートしてまわるとき、速い人は遅い人を4分ごとに追い越す。2人の毎分の速さを求めよ
答え 440m. 320m. 同じく過程式が知りたいです。
よろしくお願いします
問題①
ある仕事を仕上げるのにAは10日、Bは15日、Cは20日かかる。この仕事を3人で協力して行ったが、途中Aが休んだので、仕上げるのに6日かかった。Aは何日仕事を休んだか？
答え 3日ですが、 過去問に過程式がないので理解できないです。
計算しやすいように仕事量を60u(unitの略)とすると
A,B,C が1日にこなす仕事量は6u,4u,3uとなる。
休んだ日数をxとすると。
6u*(6-x)+4u*6+3u*6=60u
2人が自転車に乗って走っている。その速さの比は11:8である。この2人が周囲480mの円形の池の同じ地点から同時に同方向にスタートしてまわるとき、速い人は遅い人を4分ごとに追い越す。2人の毎分の速さを求めよ
答え 440m. 320m. 同じく過程式が知りたいです。
よろしくお願いします
a=11u m/min
b=8u m/min
とおいて
480/4=(11-8)u
u=40
Aが一日で片づける仕事を３０とすると１０日で３００
Bが一日で片づける仕事は２０
Cが一日で片づける仕事は１５となる
三人がフルで６日間働くと
（３０＋２０＋１５）ｘ６＝６５ｘ６＝３９０の仕事量
『この仕事』の仕事量はAの１０日分で３００
本来３９０できたはずの仕事が３００しかできなかったので
差分は９０
Aが休んだ日数は
９０／３０＝３で三日間となる
任意の整数nに対してf(n)は整数か。
素早い回答ありがとうございます。すごく助かります！
（2）√2+√3とπの大小を比較せよ。
Ax^2 + B^x + C = 0 ･･･②
R = (aC - Ac)^2 - (aB - Ab)(bC - Bc)
①、②が実数係数の2次方程式でいずれも2実数解をもつとする。①の2解はα、β；②の2解はγ、δ。
このとき、
R = a^2(Aα^2 + Bα + C)(Aβ^2 + Bβ + C) = A^2(aγ^2 + bγ + c)(aδ^2 + bδ + c)
を示せ。
2実数解→異なる実数解
ax^2 + bx + c = 0 ･･･①
Ax^2 + Bx + C = 0 ･･･②
R = (aC - Ac)^2 - (aB - Ab)(bC - Bc)
①、②が実数係数の2次方程式でいずれも異なる2実数解をもつとする。①の2解はα、β。②の2解はγ、δ。
このとき、
R = a^2(Aα^2 + Bα + C)(Aβ^2 + Bβ + C) = A^2(aγ^2 + bγ + c)(aδ^2 + bδ + c)
を示せ。
α～δは無次元、
(aC - Ac)^2 - (aB - Ab)(bC - Bc)は[U^2]、
a^2(Aα^2 + Bα + C)(Aβ^2 + Bβ + C)、A^2(aγ^2 + bγ + c)(aδ^2 + bδ + c)は[U^4]。
(1) 90ﾟ，120ﾟ
(2) √2 + √3 > π,
あ、こっちは答だけぢゃねぇのか。
θ = 15ﾟ = 60ﾟ - 45ﾟ = 45ﾟ - 30ﾟ,
から加法公式により
sin(15ﾟ) = (√6 - √2)/4,
tan(15ﾟ) = 2 - √3,
が求まる。これらを Snellius-Huygens の式
2sinθ + tanθ > 3θ,
に入れると
(√6 - √2)/2 + (2 - √3) > π/4,
よって
√2 + √3 > 2(√6 - √2) + 4(2 - √3) > π,
(*) (√2 + √3) - 2(√6 - √2) - 4(2 - √3)
= (1/4)(√2 - 1)^2･(√3 - 1)^4･(√3 - √2)
> 0,
不等式スレ９ - 761 (3), 762
くだらない問題
計算がちょっと長いだけだった
ＡＣＢ
ＢＡＣ
ＢＣＡ
＋ＣＡＢ
--------
３１２３
Ａ，Ｂ，Ｃは？
M^* を G の最大マッチングとする。
M を G のマッチングとする。
このとき、 G には、 M に関する、点を共有しない |M^*| - |M| 個の単純な増加パスが存在することを示せ
ざっとググった感じarcsinがでてくるらしいのですができればarcsin使う方向でお願いしたいです
https://i.imgur.com/PRH9wUj.jpg
ご協力ありがとうございます！
= ∫ (1 - x) * (arcsin(x))' dx
=
(1 - x) * arcsin(x) + ∫ arcsinx dx
=
(1 - x) * arcsin(x) + x * arcsin(x) + sqrt(1 - x^2)
= arcsin(x) + sqrt(1 - x^2)
Prelude> [(a,b,c)|a<-[1..9],b<-[1..9],c<-[1..9], (a+a+b+b+c)*100+(b+c+a+c+a)*10+(c+b+c+a+b)==3123]
[(3,7,8)]
（A,B,C)＝（3,7,8）
Ａ，Ｂ，Ｃが１桁の整数とは一言も書いてないけどな
俺には配列の演算が配列の要素通しの演算になるRが使い勝手がいいな。
（）で目がチカチカするがｗ
for(A in 1:9){
for(B in 1:9){
for(C in 1:9){
if (sum((c(A,B,C)+c(A,C,B)+c(B,A,C)+c(B,C,A)+c(C,A,B))*c(100,10,1)) == 3123)
print(c(A,B,C))
}
}
}
[1] 3 7 8
わかる方詳しい解説お願いします😭✨
M^* を G の最大マッチングとする。
M を G のマッチングとする。
このとき、 G には、 M に関する、点を共有しない |M^*| - |M| 個の単純な増加パスが存在することを示せ。
x=1-ay
(2a+2)(1-ay)-y=2a+6
(-2a^2-2a-1)y=4
-2(a-1/2)^2-1/2<0
(-2a^2-2a-1, y)=(-1,-4),(-2,-2),(-4,1)
(a,x,y)=(-1,5,4),(0,1,4)
ありがとうございます！
それだと(a,x,y)=(-1/2,-3,-8)のようなaが分数の場合が考慮されてないです。
問題①
2、7、15、26、40、( )
問題②
1、2、5、10、( )、26
答えは①57 ②17 です。よろしくおねがいします。
すいません、自己解決しました。
x = √( W( 2e^{2y} )/2 )
これはyに後は数値を代入して計算ソフトなどで計算するだけでしょうか？
ランベルトのW関数についていろいろ調べたのですが数値を出す例がほとんどでした
W( 2e^{2y} )をランベルトのW関数使わずにyの関数で表す方法はないのでしょうか？
例えばW( ye^y ) = y のように
∫[0→a] f(x) dx
を求めよ。aは正の実数である。
赤玉がn個と青玉がn-k個あり、これらをでたらめに左から右に横一列に並べる。
このとき
「ある連続する4つの玉からなる部分で、左から『赤赤赤青』となっている部分が存在する」
ような確率をn,kで表せ。
Zは整数全体
50∈Z が単位元となるZ上の群構造はあるか調べよ
∫[0→a] f(x) dx = a W(a) - ∫[0, W(a)] x e^x dx
= aW - [ x e^x - e^x ]{0,W}
= a W - W e^W + e^W - 1
= a( W(a) + 1/W(a) - 1) - 1
袋から玉を無作為に取り出す操作を繰り返す。取り出した玉は袋に戻さない。
袋の中の玉で、一番はじめに赤玉がなくなった場合「勝利」とし、同様に青玉がなくなった場合「敗北」とする。
また袋の中に赤玉も青玉も残っている状態で白玉を取り出した場合、操作を終了し「引き分け」とする。
（1）c=0のとき、勝利する確率を求めよ。
（2）c=1のとき、勝利する確率を求めよ。また（1）で求めた確率との大小を比較せよ。
次項から自項を引く
① ５、８、１１、１４、（１７） だから３づつ増えている
② １、３、５、（７）、９ だから奇数の列が隠れている
y=-4/(2a^2+2a+1)
よりaが有理数であることに注意してx,yが共に整数となるようなaを探せばいい
それ以上、条件が絞れないんですか？
その場合どういうふうに探せばいいんですか？
aが分数もありえるので
解決しました。
すぐ総当たり法に頼って「解けたぞ！」は、さすがに違うだろ...と思う。
答え40%です。ちと問題の意味がわかりません。過程式をよろしくお願いします
x = (5/4) * 1.12
〔問題文〕
AB＝AC=ADである四面体ABCDにおいて、辺BCの中点をEとし、点Aから底面BCDに垂線AHを引く。
このとき、点Hは△BCDの外心であることを「三垂線の定理」を用いて証明せよ。 〔以上〕
だそうです。よろしくお願いします。
△ABH ≡ △ACH ≡ △ADH。
三垂線の定理の使いどころがわからない
どこで使うのこれ
わがんね
E はどう使うの？
全くわからんね
>>787
多分三垂線の定理を適応させるために用意したものかな？？
https://i.imgur.com/LGcBsow.jpg
手書きですまんが
＜問題＞
50人の生徒からなるクラスがある。
そのうち30人は男子、20人は女子であり、男子のうち18人、女子のうち2人は眼鏡をかけている仮定とする。
（１）男女の別、眼鏡の有無のそれぞれが持つ平均自己相互量を求めよ。
（２）男女の性別が判っているという条件のもとで、眼鏡の有無が持つ条件付き自己情報量を求めよ。
答えは、
（１）H(X) = 0.97ビット, H(Y) = 0.97ビット
（２）H(Y|X) = 0.77ビット
となっております。
得意な方がいましたら、（２）の答えを出すまでの計算過程を教えていただきたいです。
よろしくお願い致します。
単純に -(18/30)log2(18/30)-(2/20)log2-(2/20) ちゃう？
題意より AB = AC
∴ ΔABCは2等辺Δ
∴ Aから底辺BCに下した水栓は中点Eで交わる。
散水栓の定理より、Hから辺BCに下した水栓も中点Eで交わる。
∴ ΔHBCも2等辺Δ
∴ HB = HC
同様にして
HB = HC = HD
3点B,C,Dは点Hを中心とする円周上にある。
点Hは△BCDの外心である。
https://mathtrain.jp/hinetsuhi
高校生なのですが、これで分からないところがあるのですが（純粋に数学的操作なのでここで質問させていただきます）
https://i.imgur.com/yWbiiCK.png
これの「両辺で積分」とありますが、何を変数として積分しているのでしょうか？
P,V,Tの微小変化量を⊿P、⊿V、⊿Tとする、というところからのみ話を勧めてて謎なのですが
まさか何で積分してもよいということはないですよね？時間とかですか？
⊿じゃなくて、dで考えると
dP/P + γdV/V = 0
両辺に積分記号をくっつけて（積分して）
∫1/P dP + γ∫1/V dV = 0
以下略
気になるなら右辺はＣでも。
簡単な微分方程式の本（昔の高校教科書レベル）を読むとわかりやすいかも。
気持ちが悪ければΔVで割り算して、Ｖに関して積分すれば
ええんでない？
あ、それぞれ別の変数で積分してよいのですか。
難しい……
>>795
これは試してみて納得しました。難しいですね……
ありがとうございました。
なんでもいいんだけど例えば V=V(T) と置いて置換積分
∫ 1/V(T) dt
= ∫ 1/V(T) V'(T) dT
= ∫ 1/V (dV/dT) dT
= ∫ 1/V dV
Pも同様
変数の間に関係が成り立つから、実は別の変数ではないんだよなあ
でもどんなパラメータで媒介変数表示しても、結局置換積分でパラメータは見えなくなるから
別々の変数で積分したような見た目になる
変数で積分してるんじゃないよ
細かい議論はすっ飛ばして言えば
辺々を順番に足し合わせていくことで
Σ(⊿P/P + γ⊿/V) = Σ0
で、⊿→d になるように極限をとれば、（細かい議論を吹っ飛ばして）
∫記号に変わるってこと。
∫f(x)dxはf(x)をxで積分してるんじゃなくて、f(x)dx を範囲の分だけ足し合わせてる感覚。
お見事
勉強してきます😭
A^5 -5A+E=0となるときAは対角化可能であることを示せ
初めて聞く言葉なので興味が湧いて
https://logics-of-blue.com/information-theory-basic/
を読んでみた。
# https://logics-of-blue.com/information-theory-basic/
"予想がつかない→不確実性（情報エントロピー）が大きい→平均情報量も大きい"
ent <- function(x){ # 情報エントロピー（平均情報量）
x=x/sum(x)
entropy=0
for(i in x) entropy=entropy+i*(-log2(i))
return(entropy)
}
ent(c(30/50,20/50)) # gender
ent(c((18+2)/50,(50-18-2)/50)) # glass
"各々の確率分布の情報量の差分の期待値をとります
確率分布が異なっていれば、情報量があるとみなすのが、
カルバック・ライブラーの情報量です。"
rel_ent <- function(P,Q){ # 相対エントロピー
n=length(P)
if(n!=length(Q)) return(NULL)
P=P/sum(P)
Q=Q/sum(Q)
re=numeric(n)
for(i in 1:n) re[i] = Q[i]*(-log2(P[i])-(-log2(Q[i])))
return(sum(re))
}
#
"相互情報量は不確実性（情報エントロピー）の減少量とみなすことができます"
"
＜問題＞
50人の生徒からなるクラスがある。
そのうち30人は男子、 20人は女子であり、
男子のうち18人、女子のうち2人は眼鏡をかけている仮定とする。
"
30/50*ent(c(18/30,12/30)) + 20/50*ent(c(2/20,18/20))
> 30/50*ent(c(18/30,12/30)) + 20/50*ent(c(2/20,18/20))
[1] 0.7701686
Rなしで計算式を書くと
30/50 * ( 18/30*(-log2(18/30))+ 12/30*(-log2(12/30))) + 20/50 * ( 2/20*(-log2( 2/20))+ 18/20*(-log2(18/20)))
括弧を見やすくすると
30/50 * [ 18/30*{-log2(18/30)}+ 12/30*{-log2(12/30)} ] + 20/50 * [ 2/20*{-log2( 2/20)}+ 18/20*{-log2(18/20)} ]
■P1stを求める
宝一つの時の自陣当たり数
(n(n+1)/2)-1 ……①
その中での宝二個の組み合わせ数
((n(n+1)/2)-1)(((n(n+1)/2)-1)-1)/2 ……②
最終マスと①との組み合わせ数
(n(n+1)/2)-1 ……③
自陣の当たりと相手の当たりで自分が勝つ
組み合わせは②と差分の和
差分は1 3 7 13 22 34 50 70 95 125 161 203
252 308 372 444 525 615……
それを表す関数
(4n^3+6n^2-4n-3+3(-1)^n)/48
nが一つずれているのでn-1に補正
｛4(n-1)^3+6(n-1)^2-4(n-1)-3+3(-1)^(n-1)｝/48 ……④
計算知能で②x2+③+④を入力すると
P1st =｛12n^4+28n^3-42n^2-52n-3(-1)^n+51｝/48 ……⑤
全n(n+1)マスで宝二個の組合わせ数
n(n+1)｛n(n+1)-1｝/2 ……⑥
引き分け数は、n(n+1)-1と同着数の和
同着数は1 2 4 6 9 12 16 20 25……
これを表す関数は ｛2n^2-1+(-1)^(n)｝/8 ……⑦
n(n+1)-1 ……⑧
計算知能で⑦+⑧を入力すると
even =(10n^2+8n+(-1)^n-9)/8 ……⑨
計算知能で⑥-⑤-⑨を入力すると
Q1st =｛12n^4+20n^3-18n^2-20n-3(-1)^n+3｝/48
Prelude> let entropy x = sum $ map (\i -> -i*(logBase 2 i)) ( map(/sum(x)) x )
Prelude> 30/50 * entropy [18, 12] + 20/50 * entropy [2, 18]
0.7701685941085136
バカすぎてぜんぜんわからんのでお願いします。
まじすか
固有値が-1の5次のJordan cellをJとすると標数5では
J^5=-E=5J-E。
答えは5πであっていますか？
違っていたら解説お願いします。
>>810
x = 2 e^{t i} + e^{-2t i} (周長比 1:3 から 2項の向きが揃うタイミングが分かる)
r^2 = |x|^2 = 5 + 2e^{3*t i} + 2e^{-3*t i} = 5 + 4 cos(3t)
tanθ := Im{x}/Re{x} = (2s-s2)/(2c+c2)
(dθ/dt) /cosθ^2 = { 2(c-c2)(2c+c2) + 2(s+s2)(2s-s2) }/(2c+c2)^2
(dθ/dt) r^2 = 2 - 2 cos(3t)
( >>410 は θ ≠ t である事を見落としたと思われる)
S = (1/2) ∫ [0→2π]dθ r^2 =(1/2) ∫ [0→2π]dt (dθ/dt) r^2
= (1/2) ∫ [0→2π]dt (2 - 2cos(3t)) = 2π
念のためプロットしてみた
https://i.imgur.com/uCqzGmh.png
まーこんなもんじゃないでしょうか。小円の半径は√2 (面積 2π)
固有多項式が重根を持たないので最小多項式も重根を持たない。
すまん。一般のn 次だった。
x^5 - 5x +1 は最小多項式で割り切れる。
最小多項式が重根を持たないのは明らかなので対角化可能。
最小多項式重解持ち得るよん。
最小多項式で割り切れるのはわかりますが重解を持たないのは言い切れますかね？
最小多項式が重根を持てば
f(x) = x^5 - 5x +1 も重根を持つ
⇔ f(x) = 0 , f’(x) = 0 が共通解を持つ
なら大丈夫ですね。
a+c=-4/3, b+4ac+d=-2, ad+bc=4, bd=1のとき、
(a^2-b)(c^2-d)<0を示せ
ありがとうございます。まだまだ勉強が足りてませんでした。
xの多項式{(1+x)^a}{(1-x)^b}を展開したとき、係数の絶対値が最大となる項の次数をa,bで表せ。
高校数学の内容だけで解く場合はどうなりますか？
ベクトルで解くとr^2=4cos2θ+5となってしまいます。ご教示ください
(x,y) = 2* (cos(t), sin(t)) + 1* (cos(-2t), sin(-2t))
第１項を公転成分、第２項を自転成分と思ってください.
そして t は "公転角" と同時に "接触点の偏角" であり, "点 P の偏角 θ" ではない事に注意.
【自転角速度が -2 の理由】
周長比 1:3 なので 小円は計3回大円の周をナメるわけです.
つまり 1ナメ目の 公転角 t=2π/3 でPは大円と2度目の接触をします(t=0 が1度目), このとき自転角は -2*2π/3 の逆回りでと公転角の "方向" と一致するわけです.
【θとt の関係】
tanθ = y/x = (2s - s2)/(2c + c2). この両辺を t で微分 (s,s2 等の略記は省スペースのため)
[左]=(dθ/dt) ( 1 + (tanθ)^2 ) =(dθ/dt)( x^2 + y^2 )/x^2 = (dθ/dt) r^2 / x^2
[右]={ 2(c - c2)(2c + c2) +2(s + s2)(2s - s2) }/x^2 = ( 2 - cos(3t) )/ x^2
∴ (dθ/dt) r^2 = 2 - cos(3t)
【面積S】
微小三角形(面積: (1/2)*r*rΔθ) の極限和を求めればよいので,
S = (1/2) ∫ [θ:0→2π] dθ r^2 = (1/2) ∫ [t:0→2π] dt (dθ/dt) r^2 = (1/2) ∫ [t:0→2π] dt ( 2 - cos(3t) ) = 2π
面積だけ求めたいのなら (x, y) や r^2 を偏角 θ で表す必要は無いのです. (簡単な形にはならない気がする)
(積分の変数変換の辺りが高校数学範囲内なのかは知らない)
面白スレの解答は…
f(x) = (x^2 +2ax+b)(x^2 +2cx+d) = x^4 -(8/3)x^3 -2x^2 +8x+1,
(a^2 -b)(c^2 -d) = (7/6 -f/2)^2 -(20/27)^2/(7/6 +f/2)
= (1/3)[10 - (13+2√11)^{2/3} - (13-2√11)^{2/3}]
= -0.2376189664261441
< 0,
面白スレ２８-319,321
Ａ・Ｂ・Ｃそれぞれ1/3の確率で引けるとする。
これをできるだけ少ない計算で楽に解く方法ないですか？
400/3*6じゃだめ？
６回やればA２回B２回が期待できる。
300×2+100×2=800(点)
ありがとうございました。
1回だけ引いた場合の期待値×繰り返す回数って計算でいいんですか？
これってカードが4枚や5枚になったり、点数が変わっても同じですか？
なるほどそのf(x)の係数になっているのか
だとするとf(x)=0の4つの解が異なる2実解と互いに共役な複素数解であることを使えば
もっと簡単に導けるな
毎回同じ条件（箱から1枚引いては戻す）場合はそう。「反復試行」と呼び、「二項分布」に従う。高校数Ｂでやるはず。教科書にものってんじゃないかな
ので高校数学の範囲では期待値求める問題でないし、期待値に関する公式も原則使えない。
どうでもいいですが～♬
平均値って期待値じゃないの？
統計でどう教えるんだろ？
Ａ・Ｂ・Ｃそれぞれ1/3の確率で引けるとする。
この方が面白いね。
これも期待値８００でいいかな？
a<bのとき、以下の不等式が成り立つことを示せ。
(b/a)^2 < (1-cosβ)/(1-cosα) < (β/α)^2
これを平面図形で示せといわれたのですが分かりません。
テスト勉強しているのですが2.(3)が分からないのでどなたかご教示下さい
原点で不連続である。
ε ∈ (0, 1) とする。
δ を任意の正の実数とする。
P = ((δ/2) * cos(δ/2), (δ/2) * sin(δ/2))
とすると、
原点と点 P の距離は、 δ/2 であり、 δ よりも小さい。
|f(P) - f(0)| = |1 - 0| = 1 > ε
James R. Munkres 著『Analysis on Manifolds』のp.9 Exercise 3に例があります。
(a^2 -b)(c^2 -d) = (7/6 -f/2)^2 - (20/27)^2/(7/6 +f/2)
= - (1/3) [(13+2√11)^{1/3} - (13-2√11)^{1/3}]^2
= - 0.2376189664261441
< 0,
conv(A) = A なら A 自体がその凸集合である。
どこに悩む要素があるのか...
「11から20までの整数のうち、連続する自然数の和では表せない
ものをすべてあげなさい。」
うまい解き方あるのかな？奇数が連続する2つの自然数の和に
なることはほとんど自明だから、偶数だけチェックすればいいけど。
自作問題：素数が3つ以上の連続する自然数の和では表せないことを示せ。
有名どこでは
nがa～bの和なら2n=(b+a)(b-a+1)より2nは2べきでなくnも2べきでない。
逆にnが2べきでないなら2nも2べきでなく2n=xy、x>y、x、yの奇遇がことなるを満たすものがとれてnはa=(x-y+1)/2～
b=(x+y-1)/2の和になるってのがあるね。
20,18,14も候補から外れるな
14は残るか
2+3+4+5=14
3+4+5=12
(8+5-1)/2=6
(8-5+1)/2=2
20=2+3+4+5+6
1～:3 6 10 15 21 28 35 45 55
2～: 5 9 14 20 27 34 44 54
3～: 7 12 18 25 32 42 52
4～: 9 15 22 29 39 49
5～: 11 18 25 35 45
6～: 13 20 30 40
7～: 14 24 34
8～: 17 27
9～: 19
出てこないのは1,2,4,8,16。
その個数が奇数の場合、個数をa、真ん中の数をbとしてS=abと表され、
個数が偶数の場合、個数を2b、真ん中の2つの数の和をaとしてS=abと表される。
いずれの場合もaは3以上の奇数。よって、Sは必ず3以上の奇数を約数として持つ。
（すなわち、2以外の素因数を持つ）
逆に、Sが3以上の奇数の約数aを持っていれば、S=abと分解した上で、
そのa,bを用いて上記2通りのアプローチで少なくとも連続した2個以上の
「整数」の和で表すことができる。
そして、それが2個以上の「自然数」の和となる条件を調べると、
２つのアプローチの片方が必ず実現可能であることがわかる。
よって、Sが2個以上の連続した自然数の和で表されるための必要十分条件は
Sが2以外の素因数を持つこと。
１６だけが表せないでいいのかな？
opをこの空間の開核作用素
clをこの空間の閉包作用素
とし、op,clをP(X)からP(X)への写像とみなす。(P(X)はXの巾集合)
この写像の合成についてなりたつ式って何でしたっけ？
op･ｃｌ･op = op だったっけ？
op･cl･op･cl = op･cl だったっけ？
上はダメ
反例R＼{0}わ。
下は言える。
閉集合 F に対し ici F= i F が言えれば良い。
ci F ⊂ F ゆえ ici F ⊂ i F。
i F ⊂ ci F ゆえ i F ⊂ ici F。
じゃあ同様の議論で
cl･op･cl･op = cl･op も言えそうだな
詳しい解答解説お願いします。
絶対値外して部分積分
そんなこともできないのかゴミ野郎w
x≦-1 のとき
F(x) = exp(x){(e-1)x +1} < 0,
x≧0 のとき
F(x) = exp(-x-1){(e-1)x +(e-2)} > 0,
-1≦x≦0 のとき
F(x) = ∫[x,0] t･exp(t) dt + ∫[0,x+1] t･exp(-t) dt
= {(1-x)exp(x) - 1} + {1 - (2+x)exp(-x-1)}
= (1-x)exp(x) - (2+x)exp(-x-1),
F(-1/2) = 0,
点(-1/2，0) について対称
F '(x) = (x+1)exp(-|x+1|) - x･exp(-|x|) = 0,
より
x = -e/(e-1) で最小 { F(x) = -(e-1)exp(-e/(e-1)) }
x = 1/(e-1) で最大 { F(x) = (e-1)exp(-e/(e-1)) }
ありがとうございます。
ゴミですいません。
A, B⊂Xについて
δ(A∪B)≦d(A, B)+δ(A)+δ(B)
であることの証明が考えつきません。
どなたかお教えいただければ幸いです。
なお、A, B⊂Xに対して
δ(A)=sup{d(x, y)|x, y∈A}
d(A, B)=inf{d(x, y)|x∈A, y∈B}
A, B⊂Xについて
δ(A∪B)≦d(A, B)+δ(A)+δ(B)
であることの証明が考えつきません。
どなたかお教えいただければ幸いです。
なお、A, B⊂Xに対して
δ(A)=sup{d(x, y)|x, y∈A}
d(A, B)=inf{d(x, y)|x∈A, y∈B}
ひぇー、即答ですね。確かにその通りですね。恐れ入りました。
a_1∈A, b_1∈B s.t. δ(A∪B) = d(a_1, b_1)
a_0∈A, b_0∈B s.t. d(A, B) = d(a_0, b_0)
となっているとする。
三角不等式より
d(A∪B) = d(a_1, b_1) ≦ d(a_1, a_0) + d(a0, b_1)
≦ d(a_1, a_0) + d(a_0, b_0) + d(b_0, b1)
≦ δ(A) + d(A, B) + δ(B)
（最後の不等号はδの定義より）
これまたお見事ですね。
初等的に導かれて面白い問題ですが、知る人ぞ知る問題なのかな。
a_1∈A, b_1∈B s.t. δ(A∪B) = d(a_1, b_1)
a_0∈A, b_0∈B s.t. d(A, B) = d(a_0, b_0)
となっているとする。
三角不等式より
d(A∪B) = d(a_1, b_1) ≦ d(a_1, a_0) + d(a_0, b_1)
≦ d(a_1, a_0) + d(a_0, b_0) + d(b_0, b_1)
≦ δ(A) + d(A, B) + δ(B)
（最後の不等号はδの定義より）
A, Bが開集合ならa_i∈A, b_j∈Bとするとまずい。
sup, infなので。
δ(A∪B)=d(a_1, b_1)
d(A, B)=d(a_0, b_0)
となるa_0, a_1∈A, b_0, b_1∈Bの存在があやふやですよね。
これでは証明になっていないと思います。
存在があやふやだと？どこまで自分の頭で考えたんだ？
A, Bが開のときは、それらは、 Xに対するA, Bの補集合の元。
統計学でt検定ってデータでいうと1変数じゃないですか？
n変数(次元)のデータに対して各クラスに有意差があるないってどういうふうに検定したらいいですか？
つまりA, Bの元でないこともあるってことですよね笑
d(A,B)の定義より
∀(ε＞0) ∃(x'∈A, y'∈B) d(x', y') < d(A,B) + ε
(1)∀(x∈A, y∈B) d(x,y)≦ d(x,x')+d(x',y')+d(y',y) < δ(A)+δ(B)+d(A, B) + ε
δ(A),δ(B)の定義より
(2)∀(x,y ∈A) d(x,y)≦ δ(A) ≦ δ(A)+δ(B)+d(A, B)
(3)∀(x,y ∈B) d(x,y)≦ δ(B) ≦ δ(A)+δ(B)+d(A, B)
(1)-(3)より
(4)∀(x,y ∈A∪B) d(x,y) < δ(A)+δ(B)+d(A, B) + ε
δ(A∪B) の定義より
(5) ∀(ε'＞0) ∃(x'',y'' ∈A∪B) δ(A∪B) - ε' < d(x'',y'')
(4),(5)より
∀(ε, ε'＞0) δ(A∪B) - ε' < δ(A)+δ(B)+d(A, B) + ε
∴ δ(A∪B) ≦ δ(A)+δ(B)+d(A, B)
(左)
AB = c とおく。
第二余弦定理より
1-cosα = {aa-(b-c)^2} /2bc = (a-b+c)(a+b-c)/2bc,
1-cosβ = {bb-(c-a)^2} /2ac = (b+c-a)(a+b-c)/2ac,
より
aa(1-cosβ) - bb(1-cosα)
= aa{1 - (cc+aa-bb)/2ac} - bb{1 - (bb+cc-aa)/2bc}
= {a(b+c-a) - b(a-b+c)}(a+b-c) /2c
= (b-a)(a+b-c)^2 /2c
> 0, (← b-a>0)
また、
1-cosβ > 1-cosα,
cosβ < cosα,
β > α,
(右)
sin(x) は 0<x<π で上に凸だから
sin(α/2) > (1-α/β)sin(0) + (α/β)sin(β/2) = (α/β)sin(β/2),
(1-cosβ)/(1-cosα) = {sin(β/2)/sin(α/2)}^2 < (β/α)^2,
緊急です。
ほんとお願いします。
https://i.imgur.com/X8VmzMu.jpg
昨日この質問をした者です。
納得できました。
本当にありがとうございます。
お願いします
0 < (b-a)/2R (← 題意)
= sinβ - sinα (← 正弦定理)
= 2sin((β-α)/2)cos((α+β)/2) (← 和積公式)
ここで 0 < (α+β)/2 < π/2 だから
sin((β-α)/2) > 0,
β > α,
非自明な結び目(e.g. trefoil knot)の管状近傍を用意してS^3からくり抜いたものをMとおく。
∂M上のループlでそのMでのホモロジー類が0であるものをとる。
lとちょうど一個共有点を持つループをmとおく。
整数i(i≠0)を選びホモロジー類がm+|iであるループxを選ぶ。
xの帯状近傍に沿ってD^2×Iの側面∂D^2×Iを貼り付けたものをNとおく。
∂NはS^2なのでここにS^3を貼り付けたものをXとおく。
Xはホモロジー3球面になる。
証明はXの基本群をVan Kampen's theoremで計算する。
a^3+b^3=m
の時aとbを求めなさい
の 4次方程式と三次方程式の連立方程式になり、最終的に二次方程式として解けるようなのですが、うまくいきませんどなたかご教授お願いします
そもそも実数だか整数だかの条件が無い
∪_n ∩_m A_(n,m) は開集合となるか？
よろしくお願いします。
(個人的には上手な足し合わせによって開集合になりそうな気がするんだけどな)
は成り立たないね
(Un),(Vm)は自然数n,mを添え字とする開集合の列
とする時、
∪_n Un ＼ ∪_m cl(Vm)
∪_m Vm ＼ ∪_n cl(Un)
は適切な和の取り方によって同時に開集合と出来ますか？
(1) a^2+b^2= (a+b)^2 - 2ab = n
(2) a^3+b^3 = (a+b)^3 - 3(a+b)ab = m
より ab を消去
(a+b)^3 - 3(a+b)((a+b)^2 - n)/2 = m
∴ (a+b)^3 - 3n(a+b) + 2m = 0
これを (a+b)についての３次方程式として解き, (1) から ab の値を得る.
後は ２次方程式 x^2 - (a+b)x + ab = 0 を解けば (a, b) が求まる.
常套的な方法で、対称式 x=a+b、y=ab とおいてx、y の連立方程式を導き、それから2次方程式を解く、かな？
具体的には
a^2+b^2=n から x^2-2y=n
a^3+b^3=m から x(n-y)=m
この2式から y を消せば x の3次方程式 x^3-3nx+2m=0 が得られるので、それを解けばよい。
忘れてくれ
Solve[{a^2+b^2==n,a^3+b^3==m},{a,b}]
(a+b)^3 -3n(a+b) +2m = 0,
の根は
a+b = 2(√n)cosθ,
ここに cos(3θ) = -m/{n^(3/2)}
(*) 本問では n^3 - mm = (3aa-2ab+3bb)(ab)^2 > 0,
> ここに cos(3θ) = -m/{n^(3/2)}
>
a bは複素数らしいけどね。
まぁだからcosθの値域についてますます気にする必要ないんだけど。
https://i.imgur.com/VUu8EPe.jpg
答えに出てくる
Root[f, i] は方程式fのi次の根、f(x)=0のi番めの解を意味している。
で、& は無名の関数を作る記号で #1 ってのは関数の1番めの引数
(この文中の「#」は#じゃないけど#みたいな記号の意)
つまり Root[2#1^6-3#1^4n-2#1^3m+3#1^2n^2+m^2-n^3 &,1] は
方程式 2x^6 -3x^4n-2x^3m+3x^2n^2+m^2-n^3=0 の1番めの解
解説ありがとうございました。
mm-n^3≧0 のときは
x^3 -3nx +2m = 0,
の根は
a+b = - [m +√(mm-n^3)]^{1/3} - [m - √(mm-n^3)]^{1/3},
ある物理量Pに関してp1,p2,p3,・・・が与えられた時、p1でp2,p3,・・・を無次元化することを考える。
単純にp1で除せばいいのかと思いましたが、p1とp2,p3,・・・ではp1だけ符号が異なっているとします。
この場合、無次元化して正の数で表したい場合はどうすればよいのでしょうか？
説明が下手で申し訳ないです。数学なのかも怪しいですが、どなたか教えていただける方がおりましたらどうかよろしくお願いします。
連続ってどういうこと？ときいてきたのでつながっている事と答えたら
それはれん○○？？の事だと言っていたのですが
何と言っていたか思い出せないのですが何かそういう言葉はありますか？
連続は近づいていけることだと言っていました
この円周上を相異なる2つの点P,Qが、PQ=1となるように動く。
（1）PQ⊥ABのとき、PAの長さを求めなさい。ただしPA≦PBとします。
（2）A,P,B,Qをこの順に結んで出来る図形が凸四角形であるとき、その面積の最大値を求めなさい。凸四角形とは、へこんでいない四角形を指します。
どこがおかしいですか？
https://i.imgur.com/r0tQ20D.jpg
(2) Aは円孤AQ上にある. ∠ABP = ∠AOP / 2 , ∠ABQ = ∠AOQ / 2 , ∠AOP + ∠AOQ = π/3
面積[APBQ] = (1/2)*2cos(∠ABP)*2sin(∠ABP) + (1/2)*2cos(∠ABQ)*2sin(∠ABQ)
= sin(∠AOP) + sin(π/3 - ∠AOP) = 2 sin(π/6) cos( ∠AOP - π/6 ) = cos( ∠AOP - π/6 )
(1)と同じ配置 ∠AOP=π/6 にて 最大値 1 となる
o Aは円孤PQ上にある
不定積分だとしたら定数(例えば C) を追加したほうがいいでしょう。
∫ sin(x)cos(x) dx = -cos(x)^2 /2 + C
∫ sin(x)cos(x) dx = sin(x)^2 /2 + C’
どちらでもよいのです。sin(x)^2 = - cos(x)^2 + 1 ですので、定数の差が 1 だけズレているだけですね。
定積分ならその差は結果に影響しません。
あーーなるほど！ありがとうございました
Sが一意に定まらないからこういう置き方はしちゃダメってことなのね。
これで 0=1 の証明をされたら間違い箇所を訂正するのに苦労しそう。
完全に忘れてましたｗポンコツですね……
∫sinxcosx dx = S(x) とおくと
S(x) = -cosx cosx - ∫(-cosx)(-sinx)dx
= -(cosx)^2 - S(x)
よって
2S(x) = -(cosx)^2 …(1)
また
S(x) = sinx sinx - ∫(sinx)(cosx)dx
= (sin x)^2 - S(x)
よって
2S(x) = (sinx)^2 …(2)
(2)-(1)より
0 = (sinx)^2 + (cosx)^2 = 1
よって
0=1 が証明された。
さて、どこが間違っているのでしょうか。
logも関数の形によっては原始関数の表現を結構いじれますよね。同じ話だと思います。
無限次元部分空間 S ⊂ X で, 線形空間 X は有限次元 (n次元) と仮定します.
Xの基底を {e_1, e_2, ..., e_n} と置きます.
Sの基底集合から n + 1 個選択して { g_1, ...., g_[n+1] }
それぞれをX基底で展開します. g_j = Σ[i=1,n] a[i , j ] e_i
g_j [k=1...n+1] が一次独立なので、「係数行列 a[i , j] の列は一次独立」です.
しかし a[i,j] は (n,n+1)型の行列なので、それは不可能です。(∵ 例えば左基本変形による掃き出し法)
矛盾が示せたので、Xの次元は無限です.
ありがとうございます
連結で調べてみたのですが集合で使う言葉なのか難しそうですね
高校の微分とか積分で連結って重要ですか？
連結じゃないと中間値の定理が使えないけど高校数学では(多分)区間上の関数しか考えてないから知らなくても問題ないと思います
a/b:c/d=ad/bc:1ですね
a/b ÷ c/dは1あたりいくつですから、比のc/dを1にした時のもう片方が答えですね
商品の価値は、その商品を生産する為に消費される energy ( 単位は erg )
即ち、熱量 ( 単位は cal.) によって決定される。
人：Ａ が、自分が所有する商品ａよりも、人：Ｂ が所有する商品ｂの方が
価値が高いと感じ、一方、人：Ｂ は自分が所有する商品ｂよりも、人：Ａ が
所有する商品ａのほうが価値が高いと感じる時、交換が起きる。
ａとｂとが等価値であると Ａ，Ｂ の双方が感じたならば、交換など起きない。
故に、Karl H. Marx ( 1818 – 1883, 65 ↟ ) が完成したと言われている、
労働価値説は完全なる誤謬なり。■
もっとわかりやすい理由はありませんか？
(a)Cの弦ABをとり、弦ABと劣弧ABで囲まれる部分の面積がS/3となるようにする。
(b)ABを1:2に内分する点Pをとる。優弧AB上に点Qをとり、弦ABと優弧ABで囲まれる部分の面積をPQで2等分する。
割算と掛算が互いに他方の逆演算になっていると考えてみたらどうかな。
0は、ひとまず考えに入れないで。
0を入れて考えたらどうなりますか？
まずx=1/bとは「bを掛けると1になる数」のことです
つまりbx=1となる数xのことですね
当然0に何を掛けても1にはならないのでこの場合は普通考えません(だからこの段階で0を入れて考えるとどうなるか、というのはナンセンスです)
で、a/bは「a掛ける1/b」です
さて、(a/b)/(c/d)は「a/b掛ける1/(c/d)」です
c/dに何を掛けたら1になるか？もちろんd/cですよね
つまり1/(c/d)=d/cです
これにa/bを掛ければ(a/b)/(c/d)=(ad)/(bc)となります
失せろゴミクズ
鈍い愚図は一昨日来やがれ
5のべき乗って下の方の位って同じ数字の並びがよく出るよなぁと思って見てたら
5^(2^n)≡5^(2^(n-1)) (mod 10^n)
が成り立ちそうな気がしてきたんですが、成り立ちますか？証明できなくてモヤモヤしてます！
5^(2^n) - 5^(2^(n-1)) = ( 5^(2^(n-1) )*( 5^(2^(n-1) - 1 )
右辺第1乗数には 5の因子が n 個以上含まれる. (∵ 2^(n-1) ≧ n 等号は n=1,2 の時のみ)
右辺第2乗数には...
5^(2^(n-1) - 1 = (5^(2^(n-2) - 1) ( 5^(2^(n-2) + 1 )
= (5^(2^(n-3) - 1) (5^(2^(n-3) + 1)( 5^(2^(n-2) + 1 ) = ...
= (5^(2^0) - 1)(5^(2^0) + 1) (5^(2^1) + 1)....( 5^(2^(n-2) + 1 )
= 4 (5^(2^0) + 1) (5^(2^1) + 1)....( 5^(2^(n-2) + 1 )
2の因子が n+1 個含まれる. (∵ 各項のmod 4 )
よって
5^(2^n) - 5^(2^(n-1)) ≡ 0 (mod 10^n) (n=1,2,3, ... )
5^(2^n) - 5^(2^(n-1)) ≡ 0 (mod 10^(n+1)) (n=3,4,5,...)
早い！ありがとうございます！
2の因子の括り出し方に感動しました！
キチガイ、失せろ。
どなたかこの問題をお願いします
積分しようにもできませんでした
A : X上の滑らかな関数全体のなす環
M : X上のベクトル場
T_x : x∈Xでの接空間
R : 実数全体のなす加法群でa∈R, f∈Aに対しfa=f(x)aとしA加群とみる
とするとき、
A加群としてT_xとMテンソルRが同型になることのイメージを教えて下さい
小学校以来の商を求める筆算の手続きが割算だと思っているぽい
https://i.imgur.com/XAszsNu.png
オレンジ領域の面積が等しくなるようにすればよい.
α - sinα*cosα = π/3
PQ*cosγ = cosα + cosβ
PQ*sinγ = (1/3)*sinα + sinβ
RO = cosβ - sinβ*cotγ = cosβ - sinβ *(cosα + cosβ) / ((1/3)*sinα + sinβ)
β - RO*sinβ = (RO + cosα) * sinα *(1/3)
綺麗な形にはならないのでプログラム組んだり, 適当なsolver使えって問題な気がする.
PQ = 1.36907... を得た.
1.26…じゃないかな
∫0～x √(1-(ξ-1)^2) dξ = π/3
積分することは困難ではないが、そうやってできた
x-1=cos(π/3+(x-1)√(1-(x-1)^2))
のような方程式を手で解くのはちょっとやりたくない
行列
(2)がわかりませんお願いします。
https://i.imgur.com/8cBJDmR.jpg
ページがまたがったため編集したら見にくくなってしまいました…ごめんなさい
高校のときはどう計算してたっけ……
A-rEが逆行列を持たないからってやってたはずで
行列式の言葉は習わなくても、逆行列を持たない条件として|ad-bc|=0は習ってるはずだし
当時の入試でも条件を求める問題は出題されてたみたいよ
自分自身は、予備校というか塾で tr とか det とか習って行列式の性質も少し習った記憶がある。
「A加群」の定義を誤解してないか？
このとき、任意のyについてf(x-y)/f(x)→1となるだろうか（x→∞）
W.ハイゼンベルグ「部分と全体」～私の生涯の偉大な出会いと対話～
424p．4860円 山崎和夫：訳 みすず書房（1999/Nov)
http://www.msz.co.jp/book/detail/04971.html
・カスタマーレヴュー
この出版社の価格では、若い人に読んでもらえるだろうか。
たとえば、文庫本として、近づきやすい価格で発行されることを希望する。 (2006/07/04)
答え時速16kmですが、なんでそうなるかわからないです。よろしくお願いします
その差は上ったときの水流の速さと下ったときの水流の速さを足したもの
水流の速さは下ったときは上ったときの2倍ということなので（以下略
すいません、理解できまへん。
答えが間違ってるのかな？過程式よろしくおねがいします
すみません
どういうことでしょうか
x:船の速度
r:川の流れ
(x-r)*1.5=18
(x+2r)*45/60=18
方程式を解くのが面倒くさいのでWolframに計算してもらいました。
https://www.wolframalpha.com/input/?i=solve+(x-y)*1.5%3D18,+(x%2B2*y)*45%2F60%3D18+for+x,y
船の速度16km/h
川の流れ4km/h
ありがとうございます。助かりました
c=2a,d=2b,10e+f=2(10c+d)+1
循環小数pで、p=0.abcdefabcdef...と表せるものを全て求めよ。
これ方程式なしで解くなら
上りの速度18/1.5＝12
下りの速度18/(45/60)=24
これは上りで川の流速分減速、下りで川の流速の２倍の加速だから
流速の３倍差がついている。
故に川の流速は(24-12)/(1+2)=4
船の静水速度は12+4もしくは24-4*2で16と出せる。
単調増加数列 a[n] (n=1,2,...) を次の漸化式で定義する.
a[1] := 1
a[2k] := a[2k-1] exp(10) - 10
a[2k+1] := a[2k] + 10 (k = 1,2,3,...)
a[n] を用いて 単調増加関数 f を以下のように定義する.
f(x∈(-∞,a[2])) := min(1, x)
f(x∈[a[2k], a[2k+1])) :=a[2k-1]* exp(x - a[2k])
f(x∈[a[2k+1], a[2k+2])) := a[2k+1] (k=1,2,...)
長いチャージ区間と 長さ10のexpダッシュ区間が交互に現れるように作った(連続)関数です.
f(x) ≦ x なのは明らか. ( f(x) < x にしたければ 0.9 f(x) で再定義)
例えば
a[2k-1]+1 < x ≦ a[2k] の時, f(x-1) / f(x) = 1 である.
a[2k]+1 < x ≦ a[2k]+10 の時, f(x-1) / f(x) = exp(-1). である.
x はいくらでも大きく取れるので f(x- 1) / f(x) は 収束しない.
(1)固有値と固有ベクトル
わざわざ固有方程式解いてもいいんだけど、基底ベクトルが線形写像でどう動くか見ればすぐわかる
λ1 = +1, λ2 = 0, λ3 = -1
v1 = +cos(θ/2) e_x + sin(θ/2) e_z
v2 = e_y
v3 = -sin(θ/2)e_x + cos(θ/2) e_z
(2) v^t S v を新たな直交基底 v1 , v2, v3 を用いて書き直す
v = x e_x + y e_y + z e_z = X v1 + Y v2 + Z v3 と置くと
v^t S v = v^t ( 1* X v1 + 0* Y v2 + (-1)*Z v3 ) = X^2 - Z^2 = c
元の基底からの回転を考慮すると, これは双曲面(x^2 - z^2 = c) を角度 θ/2 だけ回転させた図形である
https://ja.wolframalpha.com/input/?i=x%5E2+Cos%5B%CE%B8%5D+-+z%5E2+Cos%5B%CE%B8%5D+%2B+2+x+z+Sin%5B%CE%B8%5D%3Dc++plot
双曲線が出てきた。
サンガツ、もう少し条件ないとダメそうですね
Prelude> [(a,b,c,d,e,f)|a<-[0..9],b<-[0..9],c<-[0..9],d<-[0..9],e<-[0..9],f<-[0..9],c==2*a,d==2*b,10*e+f==2*(10*c+d)+1]
[(0,0,0,0,0,1),(0,1,0,2,0,5),(0,2,0,4,0,9),(0,3,0,6,1,3),(0,4,0,8,1,7),(1,0,2,0,4,1),(1,1,2,2,4,5),(1,2,2,4,4,9),(1,3,2,6,5,3),(1,4,2,8,5,7),(2,0,4,0,8,1),(2,1,4,2,8,5),(2,2,4,4,8,9),(2,3,4,6,9,3),(2,4,4,8,9,7)]
■①^2+④でもP1stは求められる
((n(n+1)/2)-1)^2+｛4(n-1)^3+6(n-1)^2-4(n-1)-3+3(-1)^(n-1)｝/48
計算知能で①^2+④を入力すると
P1st =｛12n^4+28n^3-42n^2-52n-3(-1)^n+51｝/48
http://d.hatena.ne.jp/yoshitake-h/20111210
L^≠φとするとL∪{f(L^)}がｆ鎖になり
の証明が分かりません
よろしくお願いします
https://rio2016.5ch.net/test/read.cgi/math/1543158054/
(2)
(左辺) = cosθ･(xx-zz) ＋ 2sinθ･xz
= sin(2φ)(xx-zz) ＋ 2cos(2φ)･xz (2φ = π/2-θ)
= 2sinφ･cosφ(xx-zz) + 2{(cosφ)^2 - (sinφ)^2}xz
= 2(cosφ･x - sinφ･z)(cosφ･z + sinφ･x)
= 2 X･Z,
∴ 直角双曲線柱 (?)
pを奇素数とすると
p^{2^n} - p^{2^{n-1}} = p^{2^{n-1}} * (p^{2^{n-1}} -1),
p^{2^{n-1}} -1 = (p-1)(p+1)(p^2 +1)(p^4 +1)……(p^{2^{n-2}} +1),
ところで、p^2 ≡ 1 (mod 8) ゆえ
p^2 +1, p^4 +1, ……, p^{2^{n-2}} +1 には 2の因子が1個づつ含まれる。
(p-1)(p+1)に2の因子がm個含まれるとすると、 m≧3 ゆえ
全体では (n-2)+m 個となる。
ありがとうございました
C上を2点A(α)、B(β)が独立に動く。C上の定点P(γ)の接線をlとするとき、2αβ/(α+β)がl上にあるための条件をα、β、γで表せ。
γ= 1/2 + 1/2 * e^{iθ} とすると
L(t) = γ + t e^{iθ+ iπ/2} = γ + it (2γ - 1) (t は実数パラメータ)
2αβ/(α+β) = L(t) = γ + it (2γ - 1)
(2αβ/(α+β) - γ)/(2γ - 1) = i t
∴ Re{ (2αβ/(α+β) - γ)/(2γ - 1) } = 0
( 2αβ/(α+β) を提示された意図がわからない...もっとスッキリした形になるのだろうか)
＞γ= 1/2 + 1/2 * e^{iθ} とすると
これ結局使われてないね
いちどα,β,γぜんぶこの形式に変換してみてはどうか
「e^{iθ} に直交する e^{iθ+ iπ/2} が接線の方向ベクトルを与える」
という事を明示したかった。
いきなり L(t) = γ + it (2γ - 1) とか書いても、は？ってなりそうだし。
Re{ (2αβ/(α+β) - γ)/(2γ - 1) } = 0
|2α - 1| = |2β - 1| = |2γ - 1| = 2
二番目の式を忘れるとこだった。 これで必要十分条件となる.
これ以上シンプルになるのか知らんけど.
理解しました
円の中心をμとして、接点γにおける接線を {γ+iτ(γ-μ)|τ:実数} とするのはむしろ公式かと思っていましたが、その理由を掘り下げていたのですね
2αβ/(α+β) は、(α+β=0 の場合を除き)つねに元の円の上にある
よってシンプルに書けば 2αβ/(α+β) = γ
答え2580mなんですが、問題の意味が理解できません。過程式よろしくおねがいします。
> 2αβ/(α+β) は、(α+β=0 の場合を除き)つねに元の円の上にある
これどうやって示すのか迷ったがやっと分かった.
円反転の性質を知っていれば瞬殺だ.
q := 2αβ/(α+β) の逆数( 単位円に関する反転 & 複素共役写像 ) をとると、1/q = 1/2* (1/α + 1/β) .
この点は 直線: Re(z)=1 上にある. そして直線: Re(z)=1 上の点の逆数をとれば元の円周上に移る.
つまり q (= 1/(1/q) ) は元の円周上にある.
複素数で図形を表す問題だけれど、計算のみで解くのが困難か、難しい問題の作り方教えてください。
過程式もなにも図示できたら自ずから答えてでる。
A 1 2 3 4 5 6 7 8 9 ...,71,72,73,..,86,87,88,...,155,156
156 155 154 153 152 151 150 149 148 ...,86,85,84,..,71,70,69,..., 2, 1 B
電柱の数：72+85-1=156
(87-1)*30=2580
敢えて式を書けば (72+85-70-1)*30
空白がずれで見づらいから画像にしてみた。
http://i.imgur.com/2kMUaUs.png
ありがとうございます。助かりました
分散分析
郡内分散と群間分散の比をF分布で検定する。
この比をF-ratioと呼ぶといやらしいｗ
胡散臭いけどやばい論文です https://ieeexplore.ieee.org/document/8350369 …
Learning From Pseudo-Randomness With an Artificial Neural Network?Does God Play Pseudo-Dice?
- IEEE Journals & Magazine ニューラルネットに π や eといった数学定数とか Mersenne Twister の次の桁を予測させると統計的に優位に当たってる
3:08 - 2018年11月25日
各サイコロの個性を反映して、当たる確率が 1/6 より高くなるかどうかです。
もっと性能の良いハードウェア乱数発生器でもテストして欲しかったですね。
「あらゆる状況において、人や組織がとる決断は、
（過去の状況と現在の状況は現段階では全く無関係であったとしても）
過去のにその人や組織が選択した決断によって制約を受ける」
という理論です
新しいスレッドを立ててください。
life time: 36日 0時間 33分 3秒
運営にご協力お願いいたします。
───────────────────
《プレミアム会員の主な特典》
★ 5ちゃんねる専用ブラウザからの広告除去
★ 5ちゃんねるの過去ログを取得
★ 書き込み規制の緩和
───────────────────
会員登録には個人情報は一切必要ありません。
月300円から匿名でご購入いただけます。
▼ プレミアム会員登録はこちら ▼
https://premium.5ch.net/
▼ 浪人ログインはこちら ▼
https://login.5ch.net/login.php
レス数が1000を超えています。これ以上書き込みはできません。
ニュース
- 【滋賀】いなり寿司を万引きしたとして誤認逮捕された女性が苦しい胸の内を明かす [おっさん友の会★]
- コロナワクチン接種後死亡で国を集団提訴 訴状では河野太郎大臣の「2億回打って死んでる人は1人もいない」発言に焦点 [Hitzeschleier★]
- 日経平均株価1200円安 イラン領内の爆発報道で急落… [BFU★]
- テレビ朝日「モーニングショー」 迷子の親子「顔出し放送」事情説明 「顔を映さない約束は確認できませんでした」 [ひかり★]
- 早稲田大まで行ったのに 34歳で月収12万円の後悔 貯金もゼロで「自分は何をやってたんだ」 [PARADISE★]
- 【野球】「プロ野球チップス」で誤字 日本ハム・伊藤大海投手を「176m」と記載してカルビー謝罪 [Ailuropoda melanoleuca★]
- 【フジテレビ】2024 FORMULA 1【NEXT】Lap268
- 【DAZN】フォーミュラ 【F1 F2 F3 P】lap1429
- 【WOWOW】UEFAチャンピオンズリーグ・ヨーロッパリーグ 決勝トーナメント★22
- 巨専】ワッチョイ
- とらせん
- 【D専】
- 【悲報】いなり寿司誤認逮捕、ガチで闇が深い… [207516352]
- 狂犬イスラエル、イラン・イラク・シリアを同時攻撃 [668970678]
- イスラエルのミサイル攻撃、誤報…イラン「爆発音は防空システムが作動しただけ」 [434776867]
- 


急かされて立てたお🏡



- おい！ あの美少女モデルがAVデビューしてないか？！？！ ロリコン界隈が騒然😱🤯😍😭 [485983549]
- ヒカルって実際ファン多そうだよな