結果

問題 No.186 中華風 (Easy)
ユーザー moti
提出日時 2015-04-20 00:34:00
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,179 bytes
コンパイル時間 578 ms
コンパイル使用メモリ 65,488 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-04 16:05:16
合計ジャッジ時間 1,330 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15 WA * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

#define REP(i,a,b) for(int i=a;i<(int)b;i++)
#define rep(i,n) REP(i,0,n)

typedef long long ll;

ll gcd(ll a, ll b) {
  if(b == 0) { return a; }
  return gcd(b, a%b);
}

ll extgcd(ll a, ll b, ll& x, ll& y)
{
  ll d = a;
  if(b != 0) {
    d = extgcd(b, a % b, y, x);
    y -= (a / b) * x;
  }
  else {
    x = 1; y = 0;
  }
  return d;
}

ll mod_inverse(ll a, ll m)
{
  ll x, y;
  extgcd(a, m, x, y);
  return (m + x % m) % m;
}

pair<ll, ll> linear_congruence(const vector<ll>& B, const vector<ll>& M)
{
  ll x = 0, m = 1;
  rep(i, B.size()) {
    ll a = m, b = B[i] - x, d = gcd(M[i],a);
    if(b % d != 0) { return std::make_pair(0,-1); }
    ll t = b / d * mod_inverse(a / d, M[i] / d) % (M[i] / d);
    x = x + m * t;
    m *= M[i] / d;
  }
  return std::make_pair(x % m, m);
}

int main() {

    vector<ll> B(3),M(3);
    rep(i, 3) {
      cin >> B[i] >> M[i];
    }
    
    auto res = linear_congruence(B, M);
    if(res.second == -1) { cout << -1 << endl; }
    else {
      if(res.first == 0) { cout << res.second << endl; }
      else { cout << res.first << endl; }
    }

  return 0;
}
0