素数

10兆までの素数のリストを作ってみませんか？
http://itpro.nikkeibp.co.jp/article/Watcher/20100519/348242/

参考:2桁素数
以下ソース
dim isPrime as boolean
dim s as string = ""

for i as integer = 2 to 100
　isPrime = true
　for j as integer = 2 to i-1
　　if i mod j = 0 then
　　　isPrime = false
　　　exit for
　　end if
　next j
next i
msgbox s

そうっすね

まあ素数！

落ち着かない

1 2 3 5 7 9 11 13 17 19 22

cubeスレ

素敵に見えた

cabuスレ

2 3 5 7 11 13 17 19 23
29 31 37 41 43 47 53 59 61 67
71 73 79 83 89 97 101 103 107 109
113 127 131 137 139 149 151 157 163 167
173 179 181 191 193 197 199 211 223 227
229 233 239 241 251 257 263 269 271 277
281 283 293 307 311 313 317 331 337 347
349 353 359 367 373 379 383 389 397 401
409 419 421 431 433 439 443 449 457 461
463 467 479 487 491 499 503 509 521 523
541 547 557 563 569 571 577 587 593 599
601 607 613 617 619 631 641 643 647 653
659 661 673 677 683 691 701 709 719 727
733 739 743 751 757 761 769 773 787 797
809 811 821 823 827 829 839 853 857 859
863 877 881 883 887 907 911 919 929 937
941 947 953 967 971 977 983 991 997 1009

す・・素数

落ち着いた

相手をジャガイモだと思い込む
手に人と書いてから、飲み込む
素数を数える

dim s as string = ""
これってVB6で使えたっけ？

「１」を素数に含めないのは、
ぜったい素数を考えた人の気まぐれで決めたんだと思う。

なんでVBなんだよ

.Net系っすか

マセマティカでやれ

↑素数
↓素数

67

素敵に見えた

ジョジョスレ

VBかよ

>>14
問題があるから含めないだけだよ

醤油

落ち着く〜

素数が無限にあることを証明しなさい
５分で初段

>>23
どんな問題があるんだ？

6部は黒歴史

素数だ素数を数えるんだ…

プッチ神父

>>27
数�Tの整数問題作る時に

リーマン予想とはなんだったのか

素数と素麺

>>32
NHKを見るまでサラリーマン川柳の一種かと思ってたw

プッチ神父

落ちつくなあ・・

>>1
何も表示されねえぞそれ

もとかず

素数は素敵

>>32
このスレとなんの関係があるんだよ

>>13
VB 2010だってよ

も…素数

>>41
よく見たら.netだった
最後でmsgbox使ってたからVB6かと早とちりしたわ

にしても>>1のプログラム
文字列Sをどこでも使ってないんだが

素数と来たら神父、みたいなお決まりの流れって飽きてこないか

>>43
無言で素数を数え続けるプログラムか
素敵じゃないか

P≧NPだっけ

素数やん・・・
どうせ素数の並びに規則性なんて無いんだろ
リーマンw

はよ、リーマン解けや
超絶エリートのお前らなら余裕だろ

>>43
元記事の画像
http://itpro.nikkeibp.co.jp/article/Watcher/20100519/348242/zu02.jpg

>>1がテキストに起こすときミスったんだろ

>>40
超簡単に言うと、リーマン予想は
「素数の並び方って何か規則性あんじゃね？」
ってことだ

>>49
なるほど
書き間違いか
それにしてもこんなシンプルなプログラムで
素数が導き出せるってのは素敵だな
VB2010なら並列処理も出来るから上手く書けばもっと早くなりそうだ

素数の並び方はプリズムがどうとか言ってたね

素股じゃないのか

ttp://www.nicovideo.jp/watch/sm9039359

これだな

素数だね

vbって2010が出たのかよ・・
おっつかねえよクソッ

なんだっけ、X以下の素数の個数は1/logXに比例するって予想したんだっけガウスは
数学者の自伝読むとわくわくするよな

取りあえず、Excslに計算式入力していって16127までの素数なら出した。
そこで心が折れた。手作業だとこんなもんか…

テキストで貼りつけて拡張子をhtmlで動くようにしてくれよ

このメッセージボックスは画面に収まるのか？

haskellならこれだけで出せる

main = putStrLn $ show $ take 100 $ eratosthenes [2..]

eratosthenes (x:xs) = [x] ++ (eratosthenes $ filter (\y -> mod y x /= 0) xs)
eratosthenes [] = []

そすうってなに？
きすうとちがうの？

全ての奇数は二つの素数の和で表わせて、偶数は3つ以下の素数の和で表わせるんだよね。

じゃんけんぽん
偶数＝グー
奇数＝チョキ
素数＝パー

へー。
よく分からないけど、和は綺麗なかんじだ。
覚えておこう。

落ち着くなあ

６４はすごそう。
おぼえておく。
＞６３，６４ありがとう。

>>63
少なくとも３は２つの素数に分けられないだろ

Dim isPrime As Boolean
dim s As String = ""

For i As integer =2 To 100
isOrime = True
For j As Interger =2 To i - 1
If i Mod j = 0 Then
isOrime = False
exit For
End If
Next
If isPrime Then s &= i & vbCrLf
Next
MsgBox(s)

>>63
１は偶数でも奇数でもない？

素敵やん

誰かVBSでも実行できるように書き直してくれ

>>42
もとかず　って名前いいな(・∀・)

int i, x[25]={2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97};
for(i=0;i<25;i++)
　printf("%d ", x[i]);

VBSにした

For i =2 To 100 
    isPrime = True 
    For j =2 To i - 1 
        If i Mod j = 0 Then 
            isPrime = False 
            exit For 
            End If 
    Next 
    If isPrime Then s &= i & vbCrLf 
Next 
MsgBox s 

素数って規則性発見されてないんだろ？
なのになんでコードで表せるの？

山手線ゲームで「素数」って言った理系がおしぼりぶつけられてた

新規テキスト作って、以下を書き込み保存

s = "素数="

For i = 2 To 100
isPrime = True
For j = 2 To i - 1
If i Mod j = 0 Then
isPrime = False
exit For
End If
Next
If isPrime Then s = s & "," & CSTr(i)
Next
MsgBox(s)

ファイルの拡張子を.vbsに変更
ダブルクリック

>>76
割れる数が見つかるまで延々割り続けてるだけ
規則性と言うか定義に基づいてる

>>76
規則性がある素数以外の数を抜き出すと素数が残るとかけ？

>>77
おしぼりぶつけられるまでがセットのネタなんじゃね？

>>78
97までしかでねえぞ

何かおもしろくなってきた
VB21010ダウンロードしてきて
並列処理させたプログラム書いてみるかな
1万桁ぐらいまでかけるかな？

>>82
元コードが100以下の素数だから。

100のところを自分の環境が許す限り大きくしてみれば？
文字数制限に引っかかるかなｗ

nが素数だったら、n+2が素数であることは結構ある
それと、n^2+1も素数であることが結構ある

どうして0は仲間はずれなの？

ＢＡＳＩＣの練習問題かよ

>>84
だと思って1000000にしたらなかなか出てこなかったわｗ

KEISHU5

ｗ

>>85
結構あるってなんだよナメてんのか

NHKのリーマン予想はおもしろかった
意味不明だけど

>>86
１と自分だけで割り切れる数だっつーのに、
０で割った時点でエラーだろ

エラトステネスのふるいでやると数百万くらいまでなら数秒で計算できるはず

もしスクリプトが動きっぱなしで戻ってこなくなった場合は

Ctrl + Shift + Esc　キーを同時押し
→タスクマネージャーが起動するので
　「プロセス」にタブを切り替え
→「wscript.exe」を見つけて
　右クリック「プロセスの終了」

int i,j,x[101];
for(i=2;i<101;i++)
x[i]=1;
for(i=2;i<101;i++) {
if (x[i] == 1) {
for(j=i+i;j<101;j+=i)
x[j]=0;
}
}
for(i=2;i<101;i++)
if (x[i] == 1)
printf("%d ", i);

mod とか使うと割り算やからくそとろいやないけ

簡単なものに見せたいだけでしょ

マイクロソフトさん
vbsでvbとまったく同じものが動くようにしてくれないかなぁ。

>>14
海外では１も素数に含めるとか聞いたことはあるな

一人ぐらいAppleScriptで書いてくれる奴が出てくるのを期待してる

Excel　vbaでセルに書き出してみた。自分のショボイPCだと時間がかかったな
Sub aaa()
Dim isPrime As Boolean
Dim i As Long
Dim j As Long
Dim l As Long
l = 1

For i = 2 To 1000000
isPrime = True
For j = 2 To i - 1
If i Mod j = 0 Then
isPrime = False
Exit For
End If
Next j
If isPrime Then
Cells(l, 1).Value = i
l = l + 1
End If
Next i
End Sub

ここまでで78498件でした

int i,j,x[101];
for(i=2;i<101;i++)
x[i]=1;
for(i=2;i<101;i++) {
if (x[i] == 1) {
printf("%d ", i);
for(j=i+i;j<101;j+=i)
x[j]=0;
}
}

こうできたな

いったい何をしゃべってんだお前ら

素股

レベル素数ホーリー

呼んだ？

素数って未だになんかの関数かなんかで表せなかったんだ

√i　まで調べればいいんだけどね

√i　＜　i/2　だから　i/2　まででもいい

VS2010インスコ中なんだけど
プログラム作り終わる頃にはこのスレ落ちてそうだな

>>85
n^2+1はほとんど偶数だろふざけんな

>>109
何GB？
VS2005と共存できる？

>>111
MSのサイトで落とせるISOで2.5GBぐらい
今はオンラインインストールしてるから時間かかるかも

今は新規PCにインストールしてるけど
2003と2008が共存できたから共存は出来ると思うよ

>>112
thx
ちょっと興味あるな
文法は殆ど変わってないみたいだが

#define MAX 10000000000000

long long int i,j;
_Bool x[MAX+1];

for(i=2;i<=MAX;i++)
x[i]=1;

for(i=2;i<=MAX;i++) {
if (x[i]) {
printf("%d ", i);
for(j=i+i;j<=MAX;j+=i)
x[j]=0;
}
}

こういう感じが時間最小になるハズ
メモリがあれば

2 3 5 7 9 11 13 までしかわかりません

両親が素数という言葉を発したのを聞いたことは一度もないが、普通に生きてるから、たぶん知らなくても生きていける。

素数を数えて落ち着くスレだと思ったらプログラムを組んでたでござる

>>117
おれもそう思って来たんだけどな・・・・

ただの概念上の記号の羅列だしな

>>14
ググったら素因数分解の一意性の保証のためと書いてあったぞ
本当かどうかは知らん

>>119
そんな羅列が>>54ではなぜか原子核のエネルギー準位と関係してるという

新規テキスト作って、以下を書き込み保存

dim pi

pi=4*(12*ATN(1/49)+32*ATN(1/57)-5*ATN(1/239)+12*ATN(1/110443))

MsgBox("円周率=" & pi)

ファイルの拡張子を.vbsに変更
ダブルクリック

落ち着け

VBScriptに変数宣言はいらないんじゃなかった？

人生って素数だよな



割り切れない

>>124
いらないけどね。

ついダブル指定して　→　エラー　→　あれ？いらないかでAs以降削除　→　この有様

>>125
まあ自分自身で割り切らないとやっていけないな

>>63
ゴールドバッハ予想だっけ
いつ証明されるんでしょうか、これ

>>68>>70
ゴールドバッハの予想
4以上の全ての偶数は、二つの素数の和で表すことができる。
6以上の全ての偶数は、二つの奇素数の和で表すことができる。

わたしに勇気を与えてくれる

UBASIC使えよ。いろいろ捗るぞ

まだVS2010のインストールおわらねぇｗ

それおかしくねえの？

4以上の全ての偶数はふたつの素数の和で表される。
6以上の全ての数は3つの素数の和で表される。

証明できたらほめてやるよ。

ようやく終わった
今ヘルプファイルをインストールしてる
ノートだから遅いのは仕方ないか

>>129
素数って偶数は２しかないんだから
偶数は奇数＋奇数しかないんじゃないの？
上と下の違いがわからないんだけど

そ、素数ちゃうわ！

索敵に見えた。
さようなら

そ、素数ちゃうわ

>>136
2以外のことじゃねーの

メルセンヌ素数を書き出していくと
10分くらいの暇つぶしになるよ。

素数にはロマンを感じる。

http://ja.wikipedia.org/wiki/%E7%B4%A0%E8%82%A1

しかしほんとに不思議だよな
リーマン問題が解かれても俺には意味すらわからないんだろうな

で、このスレで出てきたアルゴリズムでどれが一番早いんだよ

素数に1を含めるかどうかについて気まぐれとかいってる人がいるので。
そう思ってしまうのは

1を除く1と自分自身以外で割りきれない数

という風に教える人がいるからだと思う。
素数の一般的というかわかりやすい定義は、

約数がきっかり2つの数

というもの。8は1,2,4,8だからだめ。
11は1,11だからおｋ。1は1だかだだめぽ。
更にいうと実用上素因数分解(素数の掛け算で表すってことね)をする際に、
素数を1に含めない場合、
24=2*2*2*3
となるが、1を含めちゃうと
24=2*2*2*3=1*2*2*2*3=1*1*2*2*2*3
みたいに色々書けちゃうので面倒。

まぁ1をふくめるなら
素因数分解の際、1*というのを必ずきっかり一個だけ含めるとする
って決めごとつくればたいして困らないが、そしたら今度は素因数分解の定義が不自然になる。

リーマン予想の証明はまだ？

>>78だと全て表示されないので対策版（これも限界あるけど）
新規テキスト作って、以下を書き込み保存

s = "素数="

x = inputbox("計算範囲入力！(10000以上はきつい)",,"100")

For i = 2 To x
isPrime = True
For j = 2 To i - 1
If i Mod j = 0 Then
isPrime = False
exit For
End If
Next
If isPrime Then s = s & "," & CSTr(i)
Next
Wscript.echo(s)

ファイルの拡張子を.vbsに変更
ダブルクリック

キャッシュを超えるからなのか大きな数字で計算するとすっごい時間が掛かるな。
2万で20minぐらいか。

素数やん

やっとVS2010インスコ終わったわ

>>151
どんだけ回線遅いんだよ

VB（猛烈苦笑）

エラトステネスの篩とかイチイチカッコいいんだよ

なんか映画なかった？
立方体の箱みたいなところに突然押し込まれて
脱出方法探しますっていうかんじの

>>155
Cube

>>156
それそれあんがと
また見たくなったなあ

暗号通信用に数百桁の素数を用意してたりするけど
あれってどうやって作ってるんだろ

素数かどうか確認する必要あるだろうし
あの桁数の素数が簡単に作れて確認できたら
暗号自体の強さに影響あるような気がするんだけど

>>155
初代CUBE懐かしいな。あのサヴァンのせいで
数学得意眼鏡っ娘がアイデンティティ失ってカワイソス。

>>158
俺の古い記憶が確かなら >54 の動画内にある。
素数を管理するサーバーは解読の鍵になってるから
厳重な物理的セキュリティで守られてる。

>>158
チミが書いてるとおり
素数判定のアルゴリズムを使うんだろうな

もちろん桁数が多いもんどうしの素数同士の積なら
暗号強度も強くなる

>>27
1で割れる

FFTというソースを見つけて、wavデータから、何かの出力（実数と虚数？）を
得ることができるようになったんだけど、
これから周波数と位相ってどうやって求めるの？
最終的にはマイクから入力した音声に、アイウエオが含まれているかを知りたいんだ。

一千万まで回そうとしたけれど、5354117　,　5354119 の双子素数のあたりでギブアップｗ

>>162
１で割れるのは問題じゃない気がする。
素数同士の積が素数になってしまうのが最大の問題だと思う。
素因数分解しても、最後に１、１、１、１、１、１、１・・・・（>120既出）

>>164
平方根とってループ回数減らせば？
何ていうんだっけ？素数根の一次独立だっけ？

>>160-161
あの番組のサーバーのシーンは見たことあるけど
あれがどっかで作られた素数列を単に管理しているサーバーなのか
素数を自動的にどんどん作って独自のすごいアルゴリズムで
検証しているサーバーなのかよくわからなかった

楕円曲線法で調べると
なにか新しい発見があるかもねｗ

大胆素敵

ばかだな
・配列sosu[]作ってsosu[0]=２を最初にいれてカウンタmax=0をにする
・for idy=1 to 1000000000
・flag=false
・for idx=0 to max
・　if(.・if((idy mod sosu[idx])=0) then
・　　flag=true
・　endif
・next
・if (flag=false) then
・　max=max+1
・　sosu[max]=idy
・endif
・next

説明するつもりがつくっちまった。
配列つかって計算した素数で割り切れなかったら素数を実装してやればいい

>>167
時間かかるのは素因数分解で素数判定自体は大して時間かかんないんじゃないの

javascriptで誰かつくれ

素股に見えてちょっとわくわくしながらスレ開いた

一晩で六千万くらいまで計算して飽きた

索敵に見えた

だれかスクリプトがなにやってるかソースに日本語添えてせつめいしてくれ

>>176
ｘをｘ未満の数で割りまくって割り切れなければ素数
ってとこだろ

>>170
上でも書いたけど、なんで未満の全素数で回そうとすんの？
調べようとする数字の平方根を超えない整数まででいいのに。
101を調べるのに、11以上の素数で割る意味あんの？

偶数を除外しようとしたら
「2」を取得した時点で矛盾が発生して崩壊した

>>177
それをどうやって実装してるかしりたいんだけど

>>2
評価

プッチ神父

>>178 わからん。なんで平方根までなの？

あと実装してやったら動かないでやんの。
デバッグおわって１０００までの答えは出たけど、
さらっと書いて動かなかったんで年くったなぁとおもた。

>>180
ほら日本語で素数判定してるの探してきたやったぞ
http://www26.atwiki.jp/isoroku_be/pages/114.html

2 3 5 7 11 13 17 19 23 29
31 37 41 43 47 53 59 61 67 71
73 79 83 89 97 101 103 107 109 113
127 131 137 139 149 151 157 163 167 173
179 181 191 193 197 199 211 223 227 229
233 239 241 251 257 263 269 271 277 281
283 293 307 311 313 317 331 337 347 349
353 359 367 373 379 383 389 397 401 409
419 421 431 433 439 443 449 457 461 463
467 479 487 491 499 503 509 521 523 541
547 557 563 569 571 577 587 593 599 601
607 613 617 619 631 641 643 647 653 659
661 673 677 683 691 701 709 719 727 733
739 743 751 757 761 769 773 787 797 809
811 821 823 827 829 839 853 857 859 863
877 881 883 887 907 911 919 929 937 941
last ostr= 947 953 967 971 977 983 991 997
max= 167

あってる？

>>183
そら、２つの素数の積に限定してるから
平方根まででいいんじゃない

>>183
Q=P1×P2(P1<=P2)と素因数分解できるものとすると、
P1<=√Q
√Q<=P2が常に成り立ってるはず

P1で割っていればP2もすぐわかるはずだから常に素因数分解の小さい方だけを考えておけばいい
だから√Qより大きい数は無視していい

>>183
ＸＹ＝ＹＸだから

１０万まで計算した。１０万で１万個近くもの素数でてきた
long配列オーバーフローしそうｗｗ

・・・・
99661 99667 99679 99689 99707 99709 99713 99719 99721 99733
99761 99767 99787 99793 99809 99817 99823 99829 99833 99839
99859 99871 99877 99881 99901 99907 99923 99929 99961 99971
last ostr= 99989 99991
max= 9591

誰だよ三年くらい前に2010でVBは死滅するって言ったのは
いまだにC#は単価が高い！とかいうわけのわからん理由でVB書かされる身にもなれよMSは

俺素数

>>183
もし平方根以上の素数で割り切れた場合、
その商は調べようとした整数の平方根以下だからだよ。

例えば、面積が a の長方形がある。
この面積は縦 x と 横 y の積なんだけど、
x と y が整数になるかどうかを調べる為、
x を 2 から1ずつ上げていった場合、
どこで止める？正方形になったら(√a )でしょ？

　　　　　　＿＿＿＿＿_
　　　　 　／　　＼　 　 ／＼
　　　　／ 　し （＞） 　（＜）＼　　落ち着け　落ち着いて　素数を数えろ…
　　　　|　∪　 　 （__人__）　 J |　1、3、5、7、9、11、13
　 　　 ＼　　u　　 ｀⌒´　 　／　| |　　　　　　　　　　|
　　　　ノ　　　　　　　　　　　＼　| |　　　　　　　　　　|

>>192 あ、なんかわかってきたきがする

>>184
これって見つけた素数を暗記しないといけないやつだろ

一番は平方根以下で割りまくるってやつだな
なんでvbsとかわけわからん拡張子なんだよjsでいいだろ

素数なんかもうデフレ状態
完全数でがんばれよ

鉄格子見てたらムラムラしてきた･･･

素数判定プログラム
void sosuu(int n)
{
printf("%dは素数です\n",n);
}

Private Sub コマンド1_Click()
Dim sosu(10000) As Long
Dim max As Integer
Dim flag As Boolean
Dim ostr As String
Dim idy As Long

sosu(0) = 2
max = 0
'For idy = 1 To 1000000000
For idy = 3 To 1000
flag = False
For idx = 0 To max
If (Sqr(idy) < sosu(idx)) Then Exit For
If ((idy Mod sosu(idx)) = 0) Then
flag = True
Exit For
End If
Next
If (flag = False) Then
max = max + 1
sosu(max) = idy
End If
Next

こんでいいのかな？

あと１０兆まで難なく計算できる型と配列がほしいんだがどうすんべー

おっす！オラ素数！

1000万個目の素数を超高速に出力せよ
http://dev.team-lab.com/index.php?itemid=207
＞1位　な、なんとっ　１秒４３　chkfjさんでした
＞彼はまだ学生、、すげー奴が現れたぞ！
＞前回のリベンジおめでと！あ、おめでてー研出身だ。

お前らコイツらに勝てよ

>>199
Sqr(idy) < sosu(idx)

VBSってキャストしなくてもいいんだー。
キャスト不要ならこの判定文、sqr関数の中身みないと何とも言えないけど、
この関数呼び出すより、sosu(idx)を２乗したほうが処理が早くない？

>>203 うお、いわれてみればそうだな

>>202つかＣコンパイラ通して実行形式にしたら早いのは分かりきってるだろそれ

手詰まりだ。うちvb6(笑)だから64bit演算できなさげ　さいならー

>>203いやでもオーバーフロー考えたらキチンとキャストしてSqrとったほうがいいかもね

>>205
別に「おまえらおっそｗｗｗ」みたいなことは言ってないよ

>>208 そうかすまんな。つい意地になってしまうもんなんだこういうのはｗ

#define MAX 10000000000000 /*計算範囲の最大値*/

long long int i,j,k; /*64bitの長い変数を3個取る*/
_Bool x[MAX+1],f; /*2値フラグの配列を最大値+1個, 計算中か否かのフラグを取る*/

/*全配列をまず素数とする*/
for(i=2;i<=MAX;i++)
x[i]=True;
f=True;/*素数計算中*/

for(i=2;i<=MAX;i++) {
if (x[i]) { /* 配列i True ならば素数として出力*/
printf("%d ", i);
if (f) { /*計算中なら下記処理を行う*/
for(j=i+i,k=2;j<=MAX;j+=i,k++)
x[j]=False; /*iの倍数の配列をFalse*/
if (i >= k)/* iのk倍まで計算したkよりiが大きいなら平方根を超えたので計算終了*/
f=False;
}
}
}

>>206
多倍長演算のライブラリ作ってねばらないとｗ

>>211 おれの技量にしろ、vb6(笑)にしろ、どっちも厳しそうなので寝れなくなるわｗ

>>1
sに求めた素数格納するコード書き忘れた

マセマティカを使えば一発だぞ

それ以下の数ぶん回しよりは、計算した素数分で割ってやれば回転回数は激減するよ


代わりに格納する配列分のメモリ＆使用言語のメモリ制限＆メモリ超えた場合のＩＯ処理など
いろいろ別なところで負荷かけることになるだろうけど。

素数専用のクラスを実装すべき

// csc.exe でコンパイル

using System;

namespace prime
{
class MainClass
{
public static void Main(string[] args)
{
int x = 100000;
int i, j; bool isPrime;
for (i = 2; i < x; ++i)
{
isPrime = true;
for (j = 2; j < i - 1; ++j)
{
if ((i % j) == 0)
{
isPrime = false;
break;
}
}
if (isPrime) Console.WriteLine(i.ToString());
}
}
}
}

SIMD化しないと

やっとVB2010でプログラム書いた
>>1を元に桁増やしてマルチコア対応に書き直してみる

リーマン予想を解こうとする数学者は潰れる。

これ豆な。

大学受験のとき素数問題だけはどうしても苦手だった
あればっかりは才能だと思うわ

>>219
できあがったらアップしてくれ。

Sqr(idy)をidxループ前に変数に入れておいた方がmax回分計算しなくてすむからそっちのほうが早いな

ﾏﾀﾞｧ-? （･∀･ ）っ/凵⌒☆ﾁﾝﾁﾝ

すかずってなんだよー

5次元以上の方程式を解く公式を作れない数学者は、全員死んだ方がいいと思うよ。

五右衛門

保守

俺も保守。今Access2000でＤＢ使ってどれほどのもんができるのかやってる

俺、なにげに素数の発生する順番の法則を見つけたんだけど。
これって大発見なのかな？

素数を数えて落ち着いてるやつが思ったより少ないな・・・

>>230 いうてみ

そいつからレスが返ってくることはないと思う

>>232
あ〜！
もっと早く行ってくれればメモ紙シュレッダーにかけなかったのに〜・・・

配列の代わりにDBに素数ぶちこんでく方法つかったら遅かった。
１０万まで：素数約１万個で１０分ぐらいかかったんじゃねぇかな

一定の数までは内部メモリに貯めて
まとめて書き出した方が早いんじゃないの？

現実的には最初の１０００個を配列で保持
それ以上は１０００個単位ぐらいで区切ってファイルで保持
配列の１０００個で素数じゃないと証明できなかったら、ＩＯでそれ以上の素数を引っ張ってきてループ

というのが現実的かなぁ

>>236 AccessはＤＢがファイル単位になってるから全部メモリ上にあるようなもんだし
早いかなとおもってＤＢ試してみたんだ。
まぁなんかやり方おかしかったのかもしれないけど遅かったわ

こういう素数リストって暗号化につかうんだろうけど10兆とか民生用じゃもういらなくね？
MSのCipher.exeで必要十分なのに、意味もなく恐怖を煽って
何の意味もないソフト買わせるような詐欺もどき商法はいい加減にした方がいい

10兆なら多倍長演算不要かぁ
でも数兆個も素数保存するのには、一体何テラバイト必要なんスか・・・

>>158
大きな数をランダムに生成して
それが素数かどうかを確率的な方法（ミラー・ラビン法）でチェックする
一回のチェックで「素数です」という結果が出ても，最悪1/4の確率で間違ってる可能性があるので
チェックを難解か繰り返して全部の結果が「素数です」だった場合に
素数と判定する
ttp://ja.wikipedia.org/wiki/%E3%83%9F%E3%83%A9%E3%83%BC-%E3%83%A9%E3%83%93%E3%83%B3%E7%B4%A0%E6%95%B0%E5%88%A4%E5%AE%9A%E6%B3%95

>>239
10兆＝10の13乗≒40 bit
これで安全性を保つのは無理っす
てか
エラトステネスのふるいで誰もが簡単に求められるような範囲の素数なら
その素数２個を掛け合わせたものは
同じように誰もが簡単に試し割り算で素因数分解できちゃう

>>240計算してみよう
１０万で１万個素数がでてきたから、適当に素数が発生する確率は求める桁数の１/１０程度ということに仮定すると
１０兆内の素数は１兆個でてくるわけで、64bit＝８バイトで表現するとすると
８兆バイト＝８０億ｋバイト＝８００万Ｍバイト＝８０００Ｇバイト＝８Ｔ
実際は１０２４でわっていってるからもう少しすくないとはおもうが

ちなみに１０兆の数字をbit列に換算して１０兆bitを順に1〜１０兆を素数かどうかを判定する真偽値として
データ作成した場合１０兆bit＝１兆２５００億byte＝１２億５０００万ＫＢ＝１２５万ＭＢ＝１２５０ＧＢ＝１.２５ＴＢ

素数かどうかしらべる範囲はその数の平方根まででしらべればいいんだよね？
なら１０兆を素数かどうか調べるには3162277までの素数調べればいいんだ。


で、いくつあんの？

計算する気がしねえｗ

>>244
大体でだけど227674個ぐらいあるのだろう、きっとｗ

64bit=8byte
８byte×２０万＝１６０万byte＝１．６Ｍbyte

あら？メモリ格納できそうじゃね？

素数素敵！

N以下の素数の数π(N)は
N/log_e(N)
で近似できる
N=10兆だとおよそ334072678387 ≒ 3340億個

>>247
充分じゃないのｗ

http://primes.utm.edu/lists/small/millions/
↑ココに素数リストのデータがあるみたいだから

もしよかったら、初期素数リストが必要ならロードして使ってみたら
（ロードを短縮するためにバイナリコンバートしたり工夫してｗ）
毎回、実行するたびに演算して素数リストを作るのも大変だろうしｗ

>>249
おれはリーマン関数を使ったけどｗ
3162277までの素数の個数でいいらしいよ

>>210
#define MAX 10000000000000 /*計算範囲の最大値*/
_Bool x[MAX+1];

こういう大きい配列はコンパイル時にエラーになるのだけど
どうしたらいいの？malloc()で領域とっても実行時にエラーに
なるし。

>>252
mallocの引数intが32bitだからじゃないのｗ

3.141516までなら覚えてる

うおーVBならとりあえずはわかるけどC++は勉強しねえとな・・・

10 cls
20 a=rnd*6+1:circle (rnd*600,rnd*440),rnd*30+10,6,,,,f,a
30 goto 20

素敵に見えた

ア・・アメちゃんくれる？

無理やり配列とループ増やして１００万まで計算してみた７万個素数みつかった
大体３０秒ぐらいかなぁ

999431 999433 999437 999451 999491 999499 999521 999529 999541 999553
999563 999599 999611 999613 999623 999631 999653 999667 999671 999683
999721 999727 999749 999763 999769 999773 999809 999853 999863 999883
last ostr= 999907 999917 999931 999953 999959 999961 999979 999983
max= 78497

エリーマイラブ素数いーと

VBはSetステートメントが分かりにくいから嫌だ

for i=0 to max のmaxが32万以上になったらどうも上手くループがまわらない気がする

Sub Macro1()
Const MAX = 200000000 '計算範囲の最大値
Dim i, j, k As Long '長い変数を3個取る
Dim x(MAX + 1),f As Boolean '2値フラグの配列を最大値+1個計算中か否かのフラグ取る
For i = 2 To MAX x(i) = True Next i'全配列をまず素数とする
f = True '素数計算中
Range("A1").Select
For i = 2 To MAX
If x(i) = True Then '配列i True ならば素数として出力
ActiveCell.FormulaR1C1 = i
If ActiveCell.Row = 65536 Then ActiveCell.Offset(-65535, 1).Activate Else ActiveCell.Offset(1, 0).Activate End If'次のセルへ
If f = True Then '計算中なら下記処理を行う
k = 2
For j = i + i To MAX Step i
x(j) = False 'iの倍数の配列をFalse
k = k + 1
Next j
If i > k Then ' iのk倍まで計算したkよりiが大きいなら平方が範囲を超えたので計算終了
f = False
End If
End If
End If
Next i
End Sub

EXCEL MACRO では2億ぐらいでメモリオーバーだな
ループの中に乗除算関数の重い処理を可能な限り入れない
ベンチマークを競う場合、必須としたい
ただ、プログラムスキルとして評価されるところ日本にあるかわからんけどな

>>262
配列いくらぐらいでとってる？
>>247の20万ぴったりぐらいだと足りないかもしれないね
少し、227674より大きめの数で余裕を持たせないとダメかもｗ

そんじゃ

実行時エラー7
メモリが不足しています

Scala で書いて

Project Eulerスレか

まだやってたのか

まだあったｗ