結果

問題 No.1255 ハイレーツ・オブ・ボリビアン
ユーザー fura
提出日時 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();
      |                ~~~~~^~~~~~~~~

ソースコード

diff # 

// 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;
}
0