結果

問題 No.2324 Two Countries within UEC
ユーザー KoseT
提出日時 2023-05-28 14:44:14
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 98 ms / 2,000 ms
コード長 1,226 bytes
コンパイル時間 2,022 ms
コンパイル使用メモリ 195,076 KB
最終ジャッジ日時 2025-02-13 11:58:17
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
 
using namespace std;

int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  long long n, m, p;
  int q;
  cin >> n >> m >> p >> q;
  vector<pair<long long, long long>> query(q);
  for (int i=0;i<q;++i) cin >> query[i].first >> query[i].second;
  /**
   * query[i].first * y (mod p) == query[i].second
   * 1 <= y <= m
  */
  /**
   * x * y = cp + f
   * x * y'= c'p+ f
   * x * (y-y') = (c-c') * p
   * 
   * if (x == d * p) => d * (y-y') = (c-c')
   * else            => x * (p *(a-a')) = (c-c') * p
  */
  for (int i=0;i<q;++i) {
    long long x = query[i].first;
    long long f = query[i].second;
    if (x % p == 0) {
      if (f == 0) cout << m << '\n';
      else cout << 0 << '\n';
    } else {
      // 1つめの解
      // y = x^(p-2) (mod p)
      // x * y = 1 (mod p)
      // x * (f*y) = f (mod p)
      long long d {x};
      long long y {1};
      for (long long ct = p-2;ct != 0;ct >>= 1) {
        if (ct & 1) y *= d;
        y %= p;
        d *= d;
        d %= p;
      }
      y *= f;
      y %= p;
      //cout << y << ':';
      long long ans {};
      if (y <= m) {
        ans = (m - y + p)/p;
        if (f == 0) --ans;
      }
      cout << ans << '\n';
    }
  }
}
0