結果

問題 No.187 中華風 (Hard)
コンテスト
ユーザー tubo28
提出日時 2015-10-23 19:24:12
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,512 bytes
コンパイル時間 555 ms
コンパイル使用メモリ 65,148 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-22 11:44:09
合計ジャッジ時間 1,423 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 5 WA * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cassert>
#include <tuple>
#include <cstdio>
#include <algorithm>
using namespace std;

typedef long long ll;
typedef __int128_t lll;

#define abs(a) (a > 0 ? (a) : -(a))
#define printf(...)

lll extgcd(lll a, lll b, lll &x, lll &y) {
  for (lll u = y = 1, v = x = 0; a; ) {
    lll q = b / a;
    swap(x -= q * u, u);
    swap(y -= q * v, v);
    swap(b -= q * a, a);
  }
  return b;
}

lll gcd(lll a, lll b){
  return b ? gcd(b,a%b) : a;
}

lll lcm(lll a, lll b){
  return a/gcd(a,b)*a;
}

lll modinv(lll x, lll m) {
  lll s, t;
  extgcd(x, m, s, t);
  return (m+s) % m;
}

lll mul(lll a, lll b, lll m){
  if(a < 0) a = m - abs(a)%m;
  else a %= m;
  if(b < 0) b = m - abs(b)%m;
  else b %= m;
  return a*b%m;
}

// x, mod
pair<lll,lll> crt(lll a1, lll m1, lll a2, lll m2){
  lll k,k2;
  lll g = extgcd(m1,m2,k,k2);
  if(abs(a1-a2)%g != 0){
    return make_pair(-1,-1);
  } else {
    lll l = m1 / g * m2;
    lll x = a1 + mul(mul((a2 - a1) / g, k, l) , m1, l);
    return make_pair(x,l);
  }
}

int n;
lll x[1010], y[1010];

pair<lll,lll> solve(){
  lll mod = 1, ans = 0;
  for(int i = 0; i < n; i++){
    tie(ans,mod) = crt(ans, mod, x[i], y[i]);
    if(ans == -1) return make_pair(-1,-1);
  }
  return make_pair(ans,mod);
}

int main(){
  while(cin >> n){
    for(int i = 0; i < n; i++){
      ll xx,yy;
      if(!(cin >> xx >> yy)) exit(0);
      x[i] = xx, y[i] = yy;
    }
    lll ans,m;
    tie(ans,m) = solve();
    if(ans == 0) ans = m;
    cout << (ll)ans << endl;
  }
}
0