結果

問題 No.3056 量子コンピュータで素因数分解 Easy
ユーザー koprickykopricky
提出日時 2020-02-06 04:55:01
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 157 ms / 2,000 ms
コード長 2,692 bytes
コンパイル時間 11,688 ms
コンパイル使用メモリ 403,568 KB
実行使用メモリ 24,360 KB
平均クエリ数 3.04
最終ジャッジ日時 2023-08-30 06:54:02
合計ジャッジ時間 12,572 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 46 ms
23,928 KB
testcase_01 AC 51 ms
23,244 KB
testcase_02 AC 92 ms
23,928 KB
testcase_03 AC 92 ms
23,304 KB
testcase_04 AC 45 ms
23,424 KB
testcase_05 AC 46 ms
23,496 KB
testcase_06 AC 48 ms
23,496 KB
testcase_07 AC 44 ms
24,060 KB
testcase_08 AC 48 ms
23,256 KB
testcase_09 AC 48 ms
23,244 KB
testcase_10 AC 50 ms
23,712 KB
testcase_11 AC 66 ms
23,484 KB
testcase_12 AC 88 ms
23,388 KB
testcase_13 AC 49 ms
24,360 KB
testcase_14 AC 157 ms
23,712 KB
testcase_15 AC 157 ms
23,304 KB
testcase_16 AC 58 ms
23,844 KB
testcase_17 AC 60 ms
23,268 KB
testcase_18 AC 62 ms
24,264 KB
testcase_19 AC 154 ms
24,252 KB
testcase_20 AC 84 ms
23,628 KB
testcase_21 AC 61 ms
24,228 KB
testcase_22 AC 70 ms
23,304 KB
testcase_23 AC 72 ms
23,640 KB
testcase_24 AC 66 ms
23,292 KB
testcase_25 AC 41 ms
23,424 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <boost/multiprecision/cpp_int.hpp>
#define ll long long
#define INF 1000000005
#define MOD 1000000007
#define EPS 1e-10
#define rep(i,n) for(int i=0;i<(int)(n);++i)
#define rrep(i,n) for(int i=(int)(n)-1;i>=0;--i)
#define srep(i,s,t) for(int i=(int)(s);i<(int)(t);++i)
#define each(a,b) for(auto& (a): (b))
#define all(v) (v).begin(),(v).end()
#define len(v) (int)(v).size()
#define zip(v) sort(all(v)),v.erase(unique(all(v)),v.end())
#define cmx(x,y) x=max(x,y)
#define cmn(x,y) x=min(x,y)
#define fi first
#define se second
#define pb push_back
#define show(x) cout<<#x<<" = "<<(x)<<endl
#define sar(a,n) {cout<<#a<<":";rep(pachico,n)cout<<" "<<a[pachico];cout<<endl;}

using namespace std;

template<typename S,typename T>auto&operator<<(ostream&o,pair<S,T>p){return o<<"{"<<p.fi<<","<<p.se<<"}";}
template<typename T>auto&operator<<(ostream&o,set<T>s){for(auto&e:s)o<<e<<" ";return o;}
template<typename S,typename T,typename U>
auto&operator<<(ostream&o,priority_queue<S,T,U>q){while(!q.empty())o<<q.top()<<" ",q.pop();return o;}
template<typename K,typename T>auto&operator<<(ostream&o,map<K,T>&m){for(auto&e:m)o<<e<<" ";return o;}
template<typename T>auto&operator<<(ostream&o,vector<T>v){for(auto&e:v)o<<e<<" ";return o;}
void ashow(){cout<<endl;}template<typename T,typename...A>void ashow(T t,A...a){cout<<t<<" ";ashow(a...);}
template<typename S,typename T,typename U>
struct TRI{S fi;T se;U th;TRI(){}TRI(S f,T s,U t):fi(f),se(s),th(t){}
bool operator<(const TRI&_)const{return(fi==_.fi)?((se==_.se)?(th<_.th):(se<_.se)):(fi<_.fi);}};
template<typename S,typename T,typename U>
auto&operator<<(ostream&o,TRI<S,T,U>&t){return o<<"{"<<t.fi<<","<<t.se<<","<<t.th<<"}";}

typedef pair<int, int> P;
typedef pair<ll, ll> pll;
typedef TRI<int, int, int> tri;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef vector<vi> vvi;
typedef vector<vl> vvl;
typedef vector<P> vp;
typedef vector<double> vd;
typedef vector<string> vs;

const int MAX_N = 100005;

using MPI = boost::multiprecision::cpp_int;

MPI mod_pow(MPI a, MPI b, const MPI& mod){
	assert(b >= 0);
	MPI res(1);
	while(b){
		if(b % 2){
			res = res * a % mod;
		}
		a = a * a % mod;
		b /= 2;
	}
	return res;
}

MPI gcd(MPI a, MPI b){
    MPI tmp;
    while(b) tmp = a, a = b, b = tmp % b;
    return a;
}

int main()
{
	MPI n;
	cin >> n;
	MPI a = 1, res, hoge, p, q;
	while(true){
		++a;
		if(n % a == 0){
			cout << "! " << a << " " << n / a << endl;
			break;
		}
		cout << "? " << a << endl;
		cin >> res;
		if(res % 2) continue;
		hoge = mod_pow(a, res / 2, n) + 1;
		if(hoge == n) continue;
		p = gcd(n, hoge);
		cout << "! " << p << " " << n / p << endl;
		break;
	}
    return 0;
}
0