結果

問題 No.551 夏休みの思い出（２）
ユーザー miyo2580miyo2580
提出日時 2024-02-21 00:54:34
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 3,002 bytes
コンパイル時間 2,444 ms
コンパイル使用メモリ 219,276 KB
実行使用メモリ 88,888 KB
最終ジャッジ日時 2024-02-21 00:55:47
合計ジャッジ時間 65,785 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge16
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 48 ms
6,676 KB
testcase_01 AC 52 ms
6,676 KB
testcase_02 AC 80 ms
6,676 KB
testcase_03 AC 97 ms
6,676 KB
testcase_04 AC 205 ms
6,676 KB
testcase_05 AC 433 ms
6,676 KB
testcase_06 AC 1,032 ms
6,676 KB
testcase_07 AC 55 ms
6,676 KB
testcase_08 AC 59 ms
6,676 KB
testcase_09 AC 68 ms
6,676 KB
testcase_10 AC 57 ms
6,676 KB
testcase_11 AC 67 ms
6,676 KB
testcase_12 AC 66 ms
6,676 KB
testcase_13 AC 63 ms
6,676 KB
testcase_14 AC 70 ms
6,676 KB
testcase_15 AC 59 ms
6,676 KB
testcase_16 AC 67 ms
6,676 KB
testcase_17 AC 137 ms
8,592 KB
testcase_18 AC 125 ms
6,676 KB
testcase_19 AC 199 ms
16,436 KB
testcase_20 AC 123 ms
6,676 KB
testcase_21 AC 200 ms
13,620 KB
testcase_22 AC 122 ms
6,676 KB
testcase_23 AC 120 ms
6,676 KB
testcase_24 AC 140 ms
8,464 KB
testcase_25 AC 133 ms
8,464 KB
testcase_26 AC 165 ms
13,236 KB
testcase_27 AC 1,753 ms
42,416 KB
testcase_28 AC 2,785 ms
88,888 KB
testcase_29 TLE -
testcase_30 AC 2,843 ms
88,888 KB
testcase_31 TLE -
testcase_32 AC 2,845 ms
88,888 KB
testcase_33 AC 2,763 ms
88,888 KB
testcase_34 AC 2,921 ms
88,888 KB
testcase_35 AC 3,961 ms
88,888 KB
testcase_36 AC 3,106 ms
88,888 KB
testcase_37 AC 2,643 ms
88,888 KB
testcase_38 AC 2,675 ms
88,888 KB
testcase_39 AC 2,736 ms
88,888 KB
testcase_40 AC 2,571 ms
88,888 KB
testcase_41 AC 2,674 ms
88,888 KB
testcase_42 AC 2,549 ms
88,888 KB
testcase_43 AC 2,585 ms
88,888 KB
testcase_44 AC 2,573 ms
88,888 KB
testcase_45 AC 2,627 ms
88,888 KB
testcase_46 AC 2,554 ms
88,888 KB
testcase_47 AC 47 ms
6,676 KB
testcase_48 AC 47 ms
6,676 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