結果

問題 No.186 中華風 (Easy)
コンテスト
ユーザー legosuke
提出日時 2019-06-05 11:09:37
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 2,220 bytes
コンパイル時間 1,740 ms
コンパイル使用メモリ 170,720 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 18:34:55
合計ジャッジ時間 2,131 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define int long long
#define pii pair<int,int>
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(c) (c).begin(),(c).end()
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
#define MINF(a) memset(a,0x3f,sizeof(a))
#define POW(n) (1LL<<(n))
#define IN(i,a,b) (a <= i && i <= b)
using namespace std;
template <typename T> inline bool CHMIN(T& a,T b) { if(a>b) { a=b; return 1; } return 0; }
template <typename T> inline bool CHMAX(T& a,T b) { if(a<b) { a=b; return 1; } return 0; }
template <typename T> inline void SORT(T& a) { sort(ALL(a)); }
template <typename T> inline void REV(T& a) { reverse(ALL(a)); }
template <typename T> inline void UNI(T& a) { sort(ALL(a)); a.erase(unique(ALL(a)),a.end()); }
const int MOD = 1000000007;
const int INF = 0x3f3f3f3f3f3f3f3f;
const double EPS = 1e-10;
/* ---------------------------------------------------------------------------------------------------- */

// 拡張gcd
int ext_gcd(int a, int b, int &x, int &y) {
  if (b == 0) {
    x = 1; y = 0; return a;
  }
  int q = a/b;
  int g = ext_gcd(b,a-q*b,x,y);
  int z = x-q*y;
  x = y; y = z;
  return g;
}

// 逆元 | a と m は互いに素
int inv_mod(int a, int m) {
  int x,y;
  ext_gcd(a,m,x,y);
  x %= m;
  if (x < 0) x += m;
  return x;
}

// 中国剰余定理
// x ≡ r (mod m) の (r,m) を返す．m = lcm(m_1,m_2,...)
// 解なしの場合 (0,-1) を返す．
pair<int,int> crt(const vector<int> &b, const vector<int> &m) {
  int r = 0, M = 1;
  for (int i = 0; i < (int)b.size(); i++) {
    int p,q;
    int d = ext_gcd(M,m[i],p,q);
    if ((b[i]-r)%d != 0) return make_pair(0,-1);
    int tmp = (b[i]-r)/d*p%(m[i]/d);
    r += M*tmp;
    M *= m[i]/d;
  }
  r %= M;
  if (r < 0) r += M;
  return make_pair(r,M);
}

signed main() {
  cin.tie(0);
  ios_base::sync_with_stdio(false);
  cout << fixed << setprecision(10);

  vector<int> b(3),m(3);
  bool all_zero = true;
  for (int i = 0; i < 3; i++) {
    cin >> b[i] >> m[i];
    if (b[i]) all_zero = false;
  }
  pii ans = crt(b,m);
  if (ans.second < 0) cout << -1 << endl;
  else if (all_zero) cout << ans.second << endl;
  else cout << ans.first << endl;

  return 0;
}
0