結果
| 問題 |
No.1255 ハイレーツ・オブ・ボリビアン
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2020-10-10 01:47:51 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,171 bytes |
| コンパイル時間 | 1,805 ms |
| コンパイル使用メモリ | 199,460 KB |
| 最終ジャッジ日時 | 2025-01-15 06:08:27 |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 14 WA * 1 |
コンパイルメッセージ
main.cpp: In function ‘void solve()’:
main.cpp:37:21: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
37 | int n; scanf("%d",&n);
| ~~~~~^~~~~~~~~
main.cpp: In function ‘int main()’:
main.cpp:71:21: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
71 | int q; scanf("%d",&q); rep(_,q) solve();
| ~~~~~^~~~~~~~~
ソースコード
// OEIS A002326
#include <bits/stdc++.h>
#define rep(i,n) for(int i=0;i<(n);i++)
using namespace std;
using lint=long long;
lint modpow(lint a,lint k,int m){
lint r=1%m;
for(lint x=a%m;k>0;k>>=1,x=x*x%m) if(k&1) r=r*x%m;
return r;
}
vector<long long> divisors(long long a){
vector<long long> res;
for(long long i=1;i*i<=a;i++) if(a%i==0) {
res.emplace_back(i);
if(i*i<a) res.emplace_back(a/i);
}
sort(res.begin(),res.end());
return res;
}
lint phi(lint a){
lint res=a;
for(lint p=2;p*p<=a;p++) if(a%p==0) {
res=res/p*(p-1);
do{ a/=p; }while(a%p==0);
}
if(a>1) res=res/a*(a-1);
return res;
}
void solve(){
int n; scanf("%d",&n);
for(lint d:divisors(phi(2*n-1))){
if(modpow(2,d,2*n-1)==1){
printf("%lld\n",d);
break;
}
}
}
int main(){
// experiment
/*
for(int n=1;n<30;n++){
auto f=[&](vector<int> p){
vector<int> q(2*n);
rep(i,n){
q[2*i]=p[i];
q[2*i+1]=p[n+i];
}
return q;
};
vector<int> p(2*n);
iota(p.begin(),p.end(),0);
auto p0=p;
int cnt=0;
while(1){
p=f(p);
cnt++;
if(p==p0) break;
}
printf("n = %2d: %d\n",n,cnt);
}
*/
int q; scanf("%d",&q); rep(_,q) solve();
return 0;
}