結果

問題 No.186 中華風 (Easy)
コンテスト
ユーザー hokuto
提出日時 2021-08-30 16:13:01
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,171 bytes
コンパイル時間 832 ms
コンパイル使用メモリ 96,496 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-23 22:33:34
合計ジャッジ時間 1,588 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 17 WA * 6
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <utility>
#include <vector>
#include <string>
#include <cmath>
#include <algorithm>
#include <functional>
#include <iomanip>
#include <map>
#include <set>
#include <queue>
using namespace std;
using ll = long long int;
#define chmax(x,y)  x = (x > (y)) ? x : (x = (y));
#define chmin(x,y)  x = (x < (y)) ? x : (x = (y));

ll mod(ll a, ll m) {
    return (a % m + m) % m;
}

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

pair<ll,ll> crt(vector<ll> &b, vector<ll> &m)
{
  ll r = 0, M = 1;
  for(int i=0;i<b.size();++i)
  {
    ll p,q;
    ll d = extgcd(M,m[i],p,q);
    if( (b[i]-r)%d != 0 )return {0, -1};
    ll tmp = (b[i]-r)/d * p  %(m[i]/d);
    r += M*tmp;
    M *= m[i]/d;
  }
  return {mod(r,M),M};
}

int main(int argc, char const *argv[])
{
  vector<ll> B,M;
  int zero = 0;
  for(int i=0;i<3;++i)
  {
    ll a,b; cin >> a >> b;
    B.emplace_back(a);
    M.emplace_back(b);
    if(a == 0)zero+=1;
  }
  auto answer = crt(B,M);
  if(zero)std::cout << answer.second << std::endl;
  else std::cout << answer.first << std::endl;
}
0