#62224 (C++11（廃止可能性あり）) No.308 素数は通れません

• 問題: No.308 素数は通れません
• ユーザー: btk
• 提出日時: 2015-12-03 10:46:15
• 言語: C++11（廃止可能性あり） (gcc 13.3.0)
• 結果: WA (最新)
  AC
  (最初)
• 実行時間: -
• コード長: 2,282 bytes
• コンパイル時間: 1,745 ms
• コンパイル使用メモリ: 168,536 KB
• 実行使用メモリ: 6,944 KB
• 最終ジャッジ日時: 2024-09-14 08:25:53
• 合計ジャッジ時間: 4,847 ms
• ジャッジサーバーID （参考情報）: judge1 / judge3

ファイルパターン	結果
other	AC * 106 WA * 1

ソースコード

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
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
#include <bits/stdc++.h>
using namespace std;
vector<bool> prime(10000,true);
typedef __int128 Int;
bool bfs(int w,int n){
    vector<bool> used(30000,false);
    used[1]=true;
    used[0]=true;
    queue<int> que;
    que.push(1);
    while(que.empty()==false){
        int v=que.front();que.pop();
        if(v%w!=1&&used[v-1]==false&&prime[v-1]==false){
            que.push(v-1);
            used[v-1]=true;
        }
        if(v%w!=0&&used[v+1]==false&&v+1<=n&&prime[v+1]==false){
            que.push(v+1);
            used[v+1]=true;
        }
        if(v+w<=n&&used[v+w]==false&&prime[v+w]==false){
            que.push(v+w);
            used[v+w]=true;
        }
        if(v-w>=1&&used[v-w]==false&&prime[v-w]==false){
            que.push(v-w);
            used[v-w]=true;
        }
    }
    return used[n];
}
void eratosthenes(){
    for(int i=2;i<10000;i++){
        if(prime[i]){
            for(int j=i*2;j<10000;j+=i){
                    prime[j]=false;
            }
        }
    }
}
Int powmod(Int a,Int n,Int Mod){
    Int res=1;
    while(n>0){
        if((n&1)>0){
            res*=a;
            res%=Mod;
        }
        a=a*a;
        a%=Mod;
        n>>=1;
    }
    return res;
}
const Int mask=(1<<30);
bool miller_rabin(Int n){
    Int d=n-1,s=0;
    while(d%2==0){
        d/=2;
        s+=1;
    }
    for(int k=0;k<10000;k++){
        Int a=0;
        for(int i=0;i<2;i++){
            a<<=30;
            Int rnd=rand();
            a+=rnd%mask;
        }
        a=a%(n-1)+1;
        //cout<<(int)a<<endl;
        if(powmod(a,d,n)!=1){
            bool b=true;
            for(Int r=0,_r=1;r<s;r++,_r<<=1){
                Int _d=d*_r;
                b&=(powmod(a,_d,n)!=n-1);
            }
            if(b)return false;
        }
    }
    return true;
}
int main(){
    string S;
    cin>>S;
    if(S.size()<4){
        eratosthenes();
        int a=stoi(S);
        for(int w = 2; w < a; w++){
            if(bfs(w,a)){
                cout<<w<<endl;
                return 0;
            }
        }
    }
    else {
        Int a=0;
        for(auto& it : S){
            a*=10;
            a+=(it-'0');
        }
        a-=8;
        if((a%8==1)&&miller_rabin(a))cout<<14<<endl;
        else cout<<8<<endl;
    }
    return 0;
}
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX