結果

問題 No.551 夏休みの思い出（２）
コンテスト
ユーザー char134217728
提出日時 2017-08-15 09:35:11
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 1,628 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,469 ms
コンパイル使用メモリ 171,672 KB
実行使用メモリ 71,424 KB
最終ジャッジ日時 2024-10-13 10:55:45
合計ジャッジ時間 26,588 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 2
other AC * 27 WA * 10 TLE * 1 -- * 9
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:28:1: warning: ISO C++ forbids declaration of ‘main’ with no type [-Wreturn-type]
   28 | main(){
      | ^~~~

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define FORR(i,a,b) for (int i=(a);i>=(b);i--)
#define pb push_back

using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef vector<int> vi;
typedef set<int> si;
const int inf = 1e9;
const int mod = 1e9+7;

ll p, q, r;
map<ll, ll> m;
ll qu[3][10000];
vector<pair<ll, int> > li;
ll powMod(ll a, ll n, ll m){
  ll p = 1;
  while(n > 0){
    if(n & 1) p = p * a % m;
    a = a * a % m;
    n >>= 1;
  }
  return p;
}

main(){
  cin.tie(0);
  ios::sync_with_stdio(false);

  cin >> p >> r >> q;
  ll inv2 = powMod(2, p - 2, p), invr = powMod(r, p - 2, p);
  FOR(i, 0, q){
    ll a, b, c, d;
    cin >> a >> b >> c;
    d = powMod(a, p - 2, p);
    b = b * d % p * inv2 % p;
    c = c * d % p;
    d = b * b % p;
    c = (d - c + p) % p;
    qu[0][i] = b;
    qu[1][i] = c;
    qu[2][i] = -1;
    if(c != 0){
      li.pb(pair<ll, int>(c, i));
    }
  }
  ll a = 1, e = min(p-1, 1000000ll);
  FOR(i, 0, e){
    m[a] = ll(i);
    a = a * r % p;
  }
  FOR(i, 0, 1+p/e){
    FOR(j, 0, li.size()){
      if(m.find(li[j].first) != m.end()){
        qu[2][li[j].second] = i*e + m[li[j].first];
      }
      li[j].first = li[j].first * a % p;
    }
  }
  FOR(i, 0, q){
    if(qu[1][i] == 0){
      cout << (p - qu[0][i])%p << endl;
    }else if(qu[2][i] % 2 != 0){
      cout << -1 << endl;
    }else{
      ll ans1, ans2;
      ans1 = (powMod(r, qu[2][i] / 2, p) - qu[0][i] + p)%p;
      ans2 = (powMod(r, qu[2][i] / 2 + (p-1) / 2, p) - qu[0][i] + p)%p;
      if(ans1 > ans2)swap(ans1, ans2);
      cout << ans1 << " " << ans2 << endl;
    }
  }
}
0