結果

問題 No.3093 Please GCD
ユーザー silv723silv723
提出日時 2022-04-01 22:54:04
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 2,000 bytes
コンパイル時間 2,256 ms
コンパイル使用メモリ 206,452 KB
実行使用メモリ 24,964 KB
平均クエリ数 6.40
最終ジャッジ日時 2024-04-30 15:20:31
合計ジャッジ時間 6,034 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 26 ms
24,836 KB
testcase_01 AC 25 ms
24,964 KB
testcase_02 AC 31 ms
24,964 KB
testcase_03 TLE -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using ll=long long;
ll mpow(ll n,ll k,ll m){
    ll ret=1;
    n%=m;
    while(k>0){
        if(k&1)ret*=n;
        n*=n;
        k>>=1;
        ret%=m;
        n%=m;
    }
    return ret;
}
bool miller_rabin(ll n){
    if(n==2) return true;
    if(n==1) return false;
    if(n%2==0) return false;
    ll m=n-1;
    ll s=0;
    ll d=0;
    while(m%2==0){
        s++;
        m/=2;
        d=m;
    }
    vector<ll> primecheck={2,325,9375,28178,450775,9780504,1795265022};
    for(int i=0;i<7;i++){
        if(n<=primecheck[i]) break;
        if(mpow(primecheck[i],d,n)!=1){
            ll cnt=d;
            ll y=mpow(primecheck[i],cnt,n);
            while(cnt!=n-1&&y!=1&&y!=n-1){
                y*=y;
                y%=n;
                cnt*=2;
            }
            if(y!=n-1&&cnt%2==0) return false;
        }
    }
    return true;
}
int main(){
	ll n;
	cin>>n;
	vector<ll> v;
	for(ll i=2;i<=n;i++){
		if(miller_rabin(i)){
			v.push_back(i);
		}
	}
  ll now=0;
  ll cnt=1;
  vector<ll> ans;
	while(now<v.size()){
      if(cnt*v[now]<=n){
        cnt*=v[now];
        now++;
      }
      else{
        cout<<"? "<<cnt<<endl;
        ll k;
        cin>>k;
        if(miller_rabin(k)) ans.push_back(k);
        for(ll i=2;i*i<=k;i++){
          if(k%i==0){
            ans.push_back(i);
            ans.push_back(k/i);
          }
        }
        cnt=1;
      }
    }
  if(cnt!=1){
    cout<<"? "<<cnt<<endl;
        ll k;
        cin>>k;
        if(miller_rabin(k)) ans.push_back(k);
        for(ll i=2;i*i<=k;i++){
          if(k%i==0){
            ans.push_back(i);
            ans.push_back(k/i);
          }
        }
  }
  ll fans=1;
  for(int i=0;i<ans.size();i++){
    fans*=ans[i];
  }
  for(int i=0;i<ans.size();i++){
    if(ans[i]*fans>n) continue;
    ll x=ans[i];
    while(x<=n){
      x*=ans[i];
    }
    x/=ans[i];
    cout<<"? "<<x<<endl;
    ll z;
    cin>>z;
    fans*=z;
    fans/=ans[i];
  }
  cout<<"! "<<fans<<endl;
}
0