結果

問題 No.187 中華風 (Hard)
コンテスト
ユーザー imgry22
提出日時 2015-04-26 20:50:10
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,073 bytes
コンパイル時間 1,271 ms
コンパイル使用メモリ 164,620 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-07-05 02:51:16
合計ジャッジ時間 2,156 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 5 WA * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

typedef long long int ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<pair<int, int> > vii;
#define rrep(i, m, n) for(int (i)=(m); (i)<(n);  (i)++)
#define erep(i, m, n) for(int (i)=(m); (i)<=(n); (i)++)
#define  rep(i, n)    for(int (i)=0; (i)<(n);  (i)++)
#define rrev(i, m, n) for(int (i)=(n)-1; (i)>=(m); (i)--)
#define erev(i, m, n) for(int (i)=(n); (i)>=(m); (i)--)
#define  rev(i, n)    for(int (i)=(n)-1; (i)>=0; (i)--)
#define vrep(i, c)    for(__typeof((c).begin())i=(c).begin(); i!=(c).end(); i++)
#define  ALL(v)       (v).begin(), (v).end()
#define mp            make_pair
#define pb            push_back
template<class T, class S> inline bool minup(T& m, S x){ return m>(T)x ? (m=(T)x, true) : false; }
template<class T, class S> inline bool maxup(T& m, S x){ return m<(T)x ? (m=(T)x, true) : false; }

const int    INF = 1000000000;
const ll     MOD = 1000000007LL;
const double EPS = 1E-12;

template<class T> T gcd(T a, T b)
{
  T tmp;

  while(b){
    tmp = a;
    a = b;
    b = tmp%b;
  }

  return a;
}

template<class T> inline T lcm(T a, T b){ return a * (b / gcd(a, b)); }

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

template<typename T> T modInverse(T a, T m)
{
  T x, y;
  extgcd(a, m, x, y);
  return (m + x % m) % m;
}

template<typename T> pair<T, T> linearCongruence(const vector<T>& A, const vector<T>& B, const vector<T>& M)
{
  T x = (T)0;
  T m = (T)1;

  rep(i, A.size()){
    T a = A[i] * m;
    T b = B[i] - A[i] * x;
    T d = gcd(M[i], a);
    if(b % d) return mp((T)-1, (T)-1);
    T t = b / d * modInverse(a / d, M[i] / d) % (M[i] / d);
    x = x + m * t;
    m *= M[i] / d;
  }
  return mp(x % m, m);
}

int N;
ll x, y;
vector<ll> A, B, M;

int main()
{
  cin >> N;

  rep(i, N){
    cin >> x >> y;
    A.pb(1LL);
    B.pb(x);
    M.pb(y);
  }

  cout << linearCongruence(A, B, M).first << endl;

  return 0;
}
0