結果

問題 No.551 夏休みの思い出（２）
ユーザー miyo2580miyo2580
提出日時 2024-02-21 00:54:34
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 3,002 bytes
コンパイル時間 2,630 ms
コンパイル使用メモリ 218,472 KB
実行使用メモリ 88,892 KB
最終ジャッジ日時 2024-09-29 04:05:24
合計ジャッジ時間 72,380 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 48 ms
6,816 KB
testcase_01 AC 51 ms
6,820 KB
testcase_02 AC 78 ms
6,816 KB
testcase_03 AC 96 ms
6,820 KB
testcase_04 AC 200 ms
6,820 KB
testcase_05 AC 428 ms
6,816 KB
testcase_06 AC 991 ms
6,820 KB
testcase_07 AC 55 ms
6,816 KB
testcase_08 AC 60 ms
6,820 KB
testcase_09 AC 66 ms
6,820 KB
testcase_10 AC 58 ms
6,816 KB
testcase_11 AC 66 ms
6,816 KB
testcase_12 AC 66 ms
6,820 KB
testcase_13 AC 61 ms
6,816 KB
testcase_14 AC 69 ms
6,816 KB
testcase_15 AC 57 ms
6,816 KB
testcase_16 AC 67 ms
6,816 KB
testcase_17 AC 148 ms
8,580 KB
testcase_18 AC 127 ms
6,816 KB
testcase_19 AC 282 ms
16,312 KB
testcase_20 AC 126 ms
6,816 KB
testcase_21 AC 240 ms
13,496 KB
testcase_22 AC 122 ms
6,816 KB
testcase_23 AC 118 ms
6,820 KB
testcase_24 AC 151 ms
8,340 KB
testcase_25 AC 145 ms
8,340 KB
testcase_26 AC 206 ms
13,112 KB
testcase_27 AC 2,189 ms
42,304 KB
testcase_28 AC 3,092 ms
88,636 KB
testcase_29 TLE -
testcase_30 AC 3,131 ms
88,768 KB
testcase_31 TLE -
testcase_32 AC 3,074 ms
88,768 KB
testcase_33 AC 2,997 ms
88,764 KB
testcase_34 AC 3,146 ms
88,764 KB
testcase_35 TLE -
testcase_36 AC 3,314 ms
88,768 KB
testcase_37 AC 2,894 ms
88,764 KB
testcase_38 AC 2,956 ms
88,636 KB
testcase_39 AC 2,991 ms
88,764 KB
testcase_40 AC 2,897 ms
88,760 KB
testcase_41 AC 3,086 ms
88,764 KB
testcase_42 AC 2,977 ms
88,764 KB
testcase_43 AC 3,110 ms
88,764 KB
testcase_44 AC 2,986 ms
88,632 KB
testcase_45 AC 2,926 ms
88,764 KB
testcase_46 AC 2,859 ms
88,764 KB
testcase_47 AC 47 ms
6,820 KB
testcase_48 AC 47 ms
6,816 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define repd(i,a,b) for (ll i=(a);i<(b);i++)
#define rep(i,n) repd(i,0,n)
#define all(x) (x).begin(),(x).end()
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; }
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; }
typedef long long ll;
typedef pair<ll,ll> P;
typedef vector<ll> vec;
using Graph = vector<vector<ll>>;
const long long INF = 1LL<<60;
const long long MOD = 1000000007;

// auto mod int
// https://youtu.be/L8grWxBlIZ4?t=9858
// https://youtu.be/ERZuLAxZffQ?t=4807 : optimize
// https://youtu.be/8uowVvQ_-Mo?t=1329 : division
int mod = 1000000007;
struct mint {
  ll x; // typedef long long ll;
  mint(ll x=0):x((x%mod+mod)%mod){}
  mint operator-() const { return mint(-x);}
  mint& operator+=(const mint a) {
    if ((x += a.x) >= mod) x -= mod;
    return *this;
  }
  mint& operator-=(const mint a) {
    if ((x += mod-a.x) >= mod) x -= mod;
    return *this;
  }
  mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this;}
  mint operator+(const mint a) const { return mint(*this) += a;}
  mint operator-(const mint a) const { return mint(*this) -= a;}
  mint operator*(const mint a) const { return mint(*this) *= a;}
  mint pow(ll t) const {
    if (!t) return 1;
    mint a = pow(t>>1);
    a *= a;
    if (t&1) a *= *this;
    return a;
  }

  // for prime mod
  mint inv() const { return pow(mod-2);}
  mint& operator/=(const mint a) { return *this *= a.inv();}
  mint operator/(const mint a) const { return mint(*this) /= a;}
};
istream& operator>>(istream& is, mint& a) { return is >> a.x;}
ostream& operator<<(ostream& os, const mint& a) { return os << a.x;}

int main()
{  
    ios::sync_with_stdio(false);
    cin.tie(0);
    ll p,r;
    cin>>mod>>r;
    ll q;cin>>q;
    unordered_map<ll,ll> mp;
    mint memo=1;
    rep(i,500)memo*=r;
    mint now=1;
    rep(i,2000000){
        now*=memo;
        if(mp.count(now.x))continue;
        mp[now.x]=i+1;
    }
    rep(qi,q){
        ll A,B,C;
        mint a,b,c;
        cin>>A>>B>>C;
        a=A;
        b=B;
        c=C;
        b/=a;
        c/=a;
        mint x=b;
        x*=b;
        x/=4;
        x-=c;
        if(x.x==0){
            mint ans=0;
            ans-=b/2;
            cout<<ans<<'\n';
            continue;
        }
        vec ans;
        repd(i,1,501){
            x*=r;
            if(mp.count(x.x)){
                ll X=mp[x.x]*500;
                X-=i;
                if(X%2!=0)continue;
                ans.push_back(X/2);
                X+=mod-1;
                ans.push_back(X/2);
            }
        }
        map<ll,ll> xs;
        rep(i,ans.size()){
            mint g=mint(r).pow(ans[i]);
            g-=b/2;
            xs[g.x]=1;
        }
        if(xs.size()==0){
            cout<<-1<<'\n';
            continue;
        }
        for(auto[X,_]:xs){
            cout<<X<<' ';
        }
        cout<<'\n';
    }
    return 0;
}
0