結果

問題 No.186 中華風 (Easy)
コンテスト
ユーザー tubo28
提出日時 2015-10-23 19:10:16
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 2,120 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 551 ms
コンパイル使用メモリ 65,496 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 18:27:16
合計ジャッジ時間 1,345 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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 k1,k2;
  lll g = extgcd(m1,m2,k1,k2);
  printf(" g = %lld, m1 = %lld, m2 = %lld\n", (ll)g, (ll)m1, (ll)m2);
  printf(" a1 = %lld, a2 = %lld, g = %lld, abs(a1-a2) = %lld, abs(a1-a2) mod g = %lld\n",
	 (ll)a1, (ll)a2, (ll)(g), (ll)abs(a1-a2), (ll)(abs(a1-a2)%g));
  if(abs(a1-a2)%g != 0){
    return make_pair(-1,-1);
  } else {
    lll l = m1 / g * m2;
    k1 *= (a2-a1)/g;
    k2 *= (a2-a1)/(-g);
    assert((a2-a1)%g == 0);
    assert((a2-a1)%(-g) == 0);
    assert(m1*k1 - m2*k2 == a2-a1);
    printf(" %lld %lld %lld %lld\n", (ll)m1,(ll)k1,(ll)m2,(ll)k2);
    lll x1 = (a1 + mul(k1,m1,l)) % l;
    lll x2 = (a2 + mul(k2,m2,l)) % l;
    printf(" x1=%lld x2=%lld\n", (ll)x1,(ll)x2);
    assert(x1 == x2);
    return make_pair(x1,l);
  }
}

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

pair<lll,lll> solve(){
  printf("---------------\n");
  lll mod = 1, ans = 0;
  for(int i = 0; i < 3; i++){
    tie(ans,mod) = crt(ans, mod, x[i], y[i]);
    printf(" ans = %lld, mod = %lld\n", (ll)ans, (ll)mod);
    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