結果

問題 No.187 中華風 (Hard)
ユーザー suisanka12suisanka12
提出日時 2023-03-23 22:36:38
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 207 ms / 3,000 ms
コード長 3,057 bytes
コンパイル時間 1,732 ms
コンパイル使用メモリ 176,812 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-18 19:42:26
合計ジャッジ時間 5,352 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 163 ms
4,348 KB
testcase_03 AC 160 ms
4,348 KB
testcase_04 AC 206 ms
4,348 KB
testcase_05 AC 206 ms
4,348 KB
testcase_06 AC 206 ms
4,348 KB
testcase_07 AC 206 ms
4,348 KB
testcase_08 AC 142 ms
4,348 KB
testcase_09 AC 142 ms
4,348 KB
testcase_10 AC 142 ms
4,348 KB
testcase_11 AC 207 ms
4,348 KB
testcase_12 AC 206 ms
4,348 KB
testcase_13 AC 65 ms
4,348 KB
testcase_14 AC 65 ms
4,348 KB
testcase_15 AC 160 ms
4,348 KB
testcase_16 AC 163 ms
4,348 KB
testcase_17 AC 2 ms
4,348 KB
testcase_18 AC 2 ms
4,348 KB
testcase_19 AC 1 ms
4,348 KB
testcase_20 AC 157 ms
4,348 KB
testcase_21 AC 2 ms
4,348 KB
testcase_22 AC 206 ms
4,348 KB
testcase_23 AC 2 ms
4,348 KB
testcase_24 AC 1 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

namespace crt{
  long long gcd(long long p, long long q){
    if(p < q)swap(p, q);
    if(q == 0)return p;
    return gcd(q, p%q);
  }
  long long extgcd(long long p,long long q,long long &s,long long &t){
    if(q==0){s=1;t=0;return p;}
    long long par;
    par = extgcd(q,p%q,t,s);
    t-=p/q*s;
    return par;
  }
  pair<long long,long long> solve2CRT(pair<long long,long long>s,pair<long long,long long>t){long long p,q,r,x,ans;r=extgcd(s.second,t.second,p,q);if((s.first-t.first)%r!=0)return make_pair(-1,-1);x=s.second/r*t.second;ans=p*(s.first-t.first)/r;return make_pair(s.second*ans+s.first,x);}
  pair<long long,long long> solveCRT(vector<pair<long long,long long>>x){auto p=make_pair(1LL,1LL);for(auto it:x){p=solve2CRT(p,it);if(p.first==-1)return p;}return p;}
  //mod mにおけるaの逆元を返す。　aとmが互いに素でないときの動作はこわい。
  long long inv(long long a,long long m){long long p,q;extgcd(a,m,p,q);return (p%m+m)%m;}
  // (b, n)でx≡b(mod n)を表し、そのvectorを受け取ってその連立を満たす最小の組(s mod mod, t mod mod)を返す。
  //要件　このアルゴリズムを使うときはすべての法が互いに素でなくてはならない
  pair<long long,long long> garner(vector<pair<long long, long long>> lis, long long mod){
    int n = lis.size();
    lis.emplace_back(0, mod);
    vector<long long> ta(n+1, 0), mooda(n+1, 1);
    for(int i = 0; i < n; i++){
      long long m = lis[i].second;
      long long a = ((lis[i].first-ta[i])*inv(mooda[i], m)%m+m)%m;
      for(int j = i+1; j < n+1; j++){
        ta[j] += a*mooda[j];
        ta[j] %= lis[j].second;
        mooda[j] *= m;
        mooda[j] %= lis[j].second;
      }
    }
    return make_pair(ta[n], mooda[n]);
  }
  //互いに素でない場合でも使えてほしいgarner
  pair<long long, long long> super_garner(vector<pair<long long, long long>> &v, long long mod){
    int n = v.size();
    for(int i = 0; i < n; i++){
      for(int j = 0; j < i; j++){
        long long g = gcd(v[i].second, v[j].second);
        if((v[i].first-v[j].first)%g != 0) return make_pair(-1, -1);
        v[i].second /= g;
        v[j].second /= g;
        long long gi = gcd(v[i].second, g);
        long long gj = g/gi;
        do{
          g = gcd(gi, gj);
          gi *= g;
          gj /= g;
        }while(g != 1);
        v[i].second *= gi;
        v[j].second *= gj;
        v[i].first %= v[i].second;
        v[j].first %= v[j].second;
      }
    }
    long long lcm = 1;
    for(int i = 0; i < n; i++){
      lcm *= v[i].second;
      lcm %= mod;
    }

    return garner(v, mod);
  }

}
int main(){
  int n;
  cin >> n;
  vector<pair<long long, long long>> v;
  bool fl = true;
  for(int i = 0; i < n; i++){
    long long x, y;
    cin >> x >> y;
    v.emplace_back(x, y);
    if(x) fl = false;
  }
  auto p = crt::super_garner(v, 1000000007);
  if(p.first == -1) cout << -1 << endl;
  else if(!fl) cout << p.first << endl;
  else cout << p.second << endl;
}
0