結果

問題 No.186 中華風 (Easy)
ユーザー tubo28tubo28
提出日時 2015-10-23 19:19:05
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,599 bytes
コンパイル時間 553 ms
コンパイル使用メモリ 66,260 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 18:27:22
合計ジャッジ時間 1,349 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 1 ms
5,376 KB
testcase_02 AC 1 ms
5,376 KB
testcase_03 AC 1 ms
5,376 KB
testcase_04 AC 1 ms
5,376 KB
testcase_05 AC 1 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 1 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 1 ms
5,376 KB
testcase_13 AC 2 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 1 ms
5,376 KB
testcase_19 AC 1 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 1 ms
5,376 KB
testcase_22 AC 1 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

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;
    // k *= (a2-a1)/g;
    // lll x = (a1 + mul(k,m1,l)) % l;
    // return make_pair(x,l);
    lll xx = a1 + mul(mul((a2 - a1) / g, k, l) , m1, l);
    return make_pair(xx,l);
  }
}

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

pair<lll,lll> solve(){
  lll mod = 1, ans = 0;
  for(int i = 0; i < 3; 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(1){
    for(int i = 0; i < 3; 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