結果
問題 | No.551 夏休みの思い出(2) |
ユーザー | btk |
提出日時 | 2017-04-05 21:18:30 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2,091 ms / 4,000 ms |
コード長 | 2,457 bytes |
コンパイル時間 | 2,991 ms |
コンパイル使用メモリ | 173,412 KB |
実行使用メモリ | 48,572 KB |
最終ジャッジ日時 | 2023-09-22 20:25:40 |
合計ジャッジ時間 | 30,173 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,380 KB |
testcase_01 | AC | 2 ms
4,380 KB |
testcase_02 | AC | 2 ms
4,380 KB |
testcase_03 | AC | 2 ms
4,380 KB |
testcase_04 | AC | 3 ms
4,380 KB |
testcase_05 | AC | 5 ms
4,384 KB |
testcase_06 | AC | 5 ms
4,376 KB |
testcase_07 | AC | 2 ms
4,380 KB |
testcase_08 | AC | 1 ms
4,376 KB |
testcase_09 | AC | 1 ms
4,380 KB |
testcase_10 | AC | 1 ms
4,380 KB |
testcase_11 | AC | 2 ms
4,384 KB |
testcase_12 | AC | 1 ms
4,380 KB |
testcase_13 | AC | 2 ms
4,384 KB |
testcase_14 | AC | 2 ms
4,380 KB |
testcase_15 | AC | 2 ms
4,380 KB |
testcase_16 | AC | 2 ms
4,380 KB |
testcase_17 | AC | 38 ms
8,640 KB |
testcase_18 | AC | 186 ms
18,908 KB |
testcase_19 | AC | 174 ms
18,132 KB |
testcase_20 | AC | 184 ms
18,756 KB |
testcase_21 | AC | 3 ms
4,376 KB |
testcase_22 | AC | 10 ms
5,480 KB |
testcase_23 | AC | 4 ms
4,380 KB |
testcase_24 | AC | 17 ms
6,096 KB |
testcase_25 | AC | 32 ms
8,384 KB |
testcase_26 | AC | 80 ms
12,692 KB |
testcase_27 | AC | 979 ms
48,364 KB |
testcase_28 | AC | 961 ms
48,288 KB |
testcase_29 | AC | 800 ms
48,356 KB |
testcase_30 | AC | 948 ms
48,192 KB |
testcase_31 | AC | 845 ms
48,072 KB |
testcase_32 | AC | 863 ms
48,440 KB |
testcase_33 | AC | 972 ms
48,188 KB |
testcase_34 | AC | 1,052 ms
48,340 KB |
testcase_35 | AC | 891 ms
48,344 KB |
testcase_36 | AC | 862 ms
48,248 KB |
testcase_37 | AC | 2,091 ms
48,508 KB |
testcase_38 | AC | 1,962 ms
48,356 KB |
testcase_39 | AC | 1,067 ms
48,296 KB |
testcase_40 | AC | 1,396 ms
48,212 KB |
testcase_41 | AC | 982 ms
48,364 KB |
testcase_42 | AC | 1,784 ms
48,392 KB |
testcase_43 | AC | 1,817 ms
47,936 KB |
testcase_44 | AC | 1,673 ms
48,432 KB |
testcase_45 | AC | 1,211 ms
48,472 KB |
testcase_46 | AC | 1,952 ms
48,572 KB |
testcase_47 | AC | 1 ms
4,376 KB |
testcase_48 | AC | 2 ms
4,380 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; typedef long long LL; #define CIN_ONLY if(1) struct cww{cww(){ CIN_ONLY{ ios::sync_with_stdio(false);cin.tie(0); } }}star; #define fin "\n" #define FOR(i,bg,ed) for(int i=(bg);i<(ed);i++) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() #define fi first #define se second #define pb push_back #define DEBUG if(0) #define REC(ret, ...) std::function<ret (__VA_ARGS__)> template <typename T>inline bool chmin(T &l,T r) {bool a=l>r;if(a)l=r;return a;} template <typename T>inline bool chmax(T &l,T r) {bool a=l<r;if(a)l=r;return a;} template <typename T> istream& operator>>(istream &is,vector<T> &v){ for(auto &it:v)is>>it; return is; } const int H=1101010; const int h=1024; unordered_map<LL,LL> giant[1024]; void giant_step(LL g,int p){ LL gg=1; REP(i,h)gg=gg*g%p; LL ggg=1; REP(i,H){ LL key=ggg; LL val=h*i%(p-1); LL id=key&1023; if(giant[id].count(key)==1)return; giant[id][key]=val; ggg=ggg*gg%p; } } LL baby_step(LL key,LL r,int p){ REP(i,h){ LL id=key&1023; if(giant[id].count(key)){ LL ah=giant[id][key]; return (ah-i+p-1)%(p-1); } key=key*r%p; } return -1; } // a x + b y = gcd(a, b) // O(log (a+b) ) LL extgcd(LL a, LL b, LL &x, LL &y) { LL g = a; x = 1; y = 0; if (b != 0) g = extgcd(b, a % b, y, x), y -= (a / b) * x; return g; } // mを法とするaの逆元 // O(log a) LL invMod(LL a,int MOD) { LL x, y; if (extgcd(a, MOD, x, y) == 1)return (x + MOD) % MOD; else return 0; // unsolvable } LL pow_mod(LL a,LL n,int MOD){ LL res=1; LL b=1; a%=MOD; while(n>=b){ if(n&b) res=(res*a)%MOD; a=(a*a)%MOD; b<<=1; } return res; } LL getX(LL g,int p,LL XX){ if(XX==0)return 0; LL x=baby_step(XX,g,p); if(x%2==0)return pow_mod(g,x/2,p); else return -1; } int main(){ int P,Q,R; cin>>P>>R>>Q; giant_step(R,P); REP(i,Q){ LL a,b,c; cin>>a>>b>>c; LL sq=getX(R,P,(b*b%P+P-4*a*c%P)%P); if(sq==-1)cout<<-1<<fin; else { LL iv=invMod(2*a%P,P); LL x=((-b-sq+P+P)%P)*iv%P; LL y=((-b+sq+P+P)%P)*iv%P; if(x==y)cout<<x<<fin; if(x>y)cout<<y<<" "<<x<<fin; if(x<y)cout<<x<<" "<<y<<fin; } } return 0; }