結果
問題 | No.551 夏休みの思い出(2) |
ユーザー | beet |
提出日時 | 2018-11-23 15:52:08 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,205 bytes |
コンパイル時間 | 2,790 ms |
コンパイル使用メモリ | 207,436 KB |
実行使用メモリ | 12,344 KB |
最終ジャッジ日時 | 2023-08-26 04:15:52 |
合計ジャッジ時間 | 17,960 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | WA | - |
testcase_02 | WA | - |
testcase_03 | WA | - |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | AC | 3 ms
4,380 KB |
testcase_08 | AC | 4 ms
4,380 KB |
testcase_09 | WA | - |
testcase_10 | AC | 4 ms
4,384 KB |
testcase_11 | WA | - |
testcase_12 | AC | 3 ms
4,384 KB |
testcase_13 | WA | - |
testcase_14 | WA | - |
testcase_15 | WA | - |
testcase_16 | AC | 3 ms
4,380 KB |
testcase_17 | AC | 50 ms
8,164 KB |
testcase_18 | AC | 51 ms
8,008 KB |
testcase_19 | AC | 52 ms
8,084 KB |
testcase_20 | AC | 51 ms
8,008 KB |
testcase_21 | AC | 52 ms
8,072 KB |
testcase_22 | AC | 45 ms
7,288 KB |
testcase_23 | AC | 51 ms
8,016 KB |
testcase_24 | AC | 51 ms
8,084 KB |
testcase_25 | AC | 50 ms
8,008 KB |
testcase_26 | AC | 54 ms
8,020 KB |
testcase_27 | AC | 1,312 ms
8,024 KB |
testcase_28 | AC | 1,303 ms
8,084 KB |
testcase_29 | AC | 228 ms
8,016 KB |
testcase_30 | AC | 1,007 ms
8,164 KB |
testcase_31 | AC | 238 ms
8,020 KB |
testcase_32 | AC | 992 ms
8,040 KB |
testcase_33 | AC | 1,300 ms
8,008 KB |
testcase_34 | AC | 786 ms
8,228 KB |
testcase_35 | AC | 269 ms
8,152 KB |
testcase_36 | AC | 504 ms
8,020 KB |
testcase_37 | TLE | - |
testcase_38 | -- | - |
testcase_39 | -- | - |
testcase_40 | -- | - |
testcase_41 | -- | - |
testcase_42 | -- | - |
testcase_43 | -- | - |
testcase_44 | -- | - |
testcase_45 | -- | - |
testcase_46 | -- | - |
testcase_47 | -- | - |
testcase_48 | -- | - |
ソースコード
#include<bits/stdc++.h> using namespace std; using Int = long long; template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;} template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;} Int MOD = 1000000007; template<typename T> struct Mint{ T v; Mint():v(0){} Mint(signed v):v(v){} Mint(long long t){v=t%MOD;if(v<0) v+=MOD;} Mint pow(long long k){ Mint res(1),tmp(v); while(k){ if(k&1) res*=tmp; tmp*=tmp; k>>=1; } return res; } Mint inv(){return pow(MOD-2);} Mint& operator+=(Mint a){v+=a.v;if(v>=MOD)v-=MOD;return *this;} Mint& operator-=(Mint a){v+=MOD-a.v;if(v>=MOD)v-=MOD;return *this;} Mint& operator*=(Mint a){v=1LL*v*a.v%MOD;return *this;} Mint& operator/=(Mint a){return (*this)*=a.inv();} Mint operator+(Mint a) const{return Mint(v)+=a;}; Mint operator-(Mint a) const{return Mint(v)-=a;}; Mint operator*(Mint a) const{return Mint(v)*=a;}; Mint operator/(Mint a) const{return Mint(v)/=a;}; Mint operator-(){return v?MOD-v:v;} bool operator==(const Mint a)const{return v==a.v;} bool operator!=(const Mint a)const{return v!=a.v;} bool operator <(const Mint a)const{return v <a.v;} // find x s.t. a^x = b static T log(Mint a,Mint b){ static T sq = 100000; static map<Mint, T> dp; if(dp.empty()){ Mint res(1); for(int r=0;r<sq;r++){ if(!dp.count(res)) dp[res]=r; res*=a; } } Mint p=a.inv().pow(sq),res=b; for(int q=0;q<=MOD/sq;q++){ if(dp.count(res)){ T idx=q*sq+dp[res]; if(idx>0) return idx; } res*=p; } return T(-1); } }; //INSERT ABOVE HERE signed main(){ using M = Mint<int>; int r; cin>>MOD>>r; M R(r); int q; cin>>q; for(int i=0;i<q;i++){ int a,b,c; cin>>a>>b>>c; M A(a),B(b),C(c); M D=(A.inv()*B*M(2).inv()).pow(2)-A.inv()*C; M E=-A.inv()*B*M(2).inv(); if(D.v==0){ cout<<E.v<<endl; continue; } int k=M::log(R,D); if(k&1){ cout<<-1<<endl; continue; } M F=R.pow(k/2); int x=(E+F).v,y=(E-F).v; if(x>y) swap(x,y); cout<<x<<" "<<y<<endl; } return 0; }