結果

問題 No.186 中華風 (Easy)
ユーザー n_vipn_vip
提出日時 2015-05-14 20:10:04
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 3,084 bytes
コンパイル時間 1,205 ms
コンパイル使用メモリ 112,156 KB
実行使用メモリ 39,368 KB
最終ジャッジ日時 2024-07-06 03:43:47
合計ジャッジ時間 3,000 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 48 ms
29,200 KB
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 44 ms
26,416 KB
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 AC 44 ms
25,924 KB
testcase_20 AC 41 ms
22,536 KB
testcase_21 AC 65 ms
39,020 KB
testcase_22 AC 61 ms
38,604 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <string>
#include <vector>
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<stack>
#include<queue>
#include<cmath>
#include<algorithm>
#include<functional>
#include<list>
#include<deque>
#include<bitset>
#include<set>
#include<map>
#include<unordered_map>
#include<cstring>
#include<sstream>
#include<complex>
#include<iomanip>
#include<numeric>
#define X first
#define Y second
#define pb push_back
#define rep(X,Y) for (int (X) = 0;(X) < (Y);++(X))
#define rrep(X,Y) for (int (X) = (Y-1);(X) >=0;--(X))
#define repe(X,Y) for ((X) = 0;(X) < (Y);++(X))
#define peat(X,Y) for (;(X) < (Y);++(X))
#define all(X) (X).begin(),(X).end()
#define rall(X) (X).rbegin(),(X).rend()

using namespace std;
typedef __int128 ll;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
template<class T> using vv=vector<vector<T>>;
template<class T> ostream& operator<<(ostream &os, const vector<T> &t) {
os<<"{"; rep(i,t.size()) {os<<t[i]<<",";} os<<"}"<<endl; return os;}
template<class S, class T> ostream& operator<<(ostream &os, const pair<S,T> &t) { return os<<"("<<t.first<<","<<t.second<<")";}


const int PRIME_MAX=1123456;
vector<int> primes;
int prime[PRIME_MAX];
void findPrime(){
  for(int i=2;i<PRIME_MAX;i++)if(!prime[i]){
      primes.pb(i);
      for(int j=i*2;j<PRIME_MAX;j+=i)
	prime[j]=1;
  }
}
//extGCD
ll extgcd(ll a,ll b,ll &x,ll &y){
  ll g=a; x=1,y=0;
  if(b){
    g=extgcd(b,a%b,y,x);
    y-=(a/b)*x;
  }
  return g;
}

//中国剰余定理
ll china2(const vector<ll> &a,const vector<ll> &m){
  ll p=1,re=0,pp,x,y;
  rep(i,m.size())
    p*=m[i];
  //cout<<p<<"!!!!!!!"<<endl;
  rep(i,m.size()){
    pp=p/m[i];
    extgcd(m[i],pp,x,y);
    (re+=pp*y%p*a[i]%p)%=p;
  }
  return (re%p+p)%p;
}
int updChina(map<ll,pll> &mp,ll a,ll m,ll p){
  //cout<<a<<","<<m<<","<<p<<endl;
  if(mp.count(p)){
    pll tmp=mp[p];
    if((tmp.X-a)%min(tmp.Y,m))
      return 1;
    if(tmp.Y<m)
      mp[p]=pll(a,m);
  }else{
    mp[p]=pll(a,m);
  }
  return 0;
}
//互いに素じゃなくてもよい。解なしは-1
ll china(const vector<ll> &a,const vector<ll> &m){
  map<ll,pll> mp;  //PRIME -> REM * MOD
  pll tmp;
  ll num,p=1;
  rep(i,a.size()){
    num=m[i];
    rep(j,primes.size())if(m[i]%primes[j]==0){
      p=1;
      while(num%primes[j]==0){
	p*=primes[j];
	num/=primes[j];
      }
      if(updChina(mp,a[i],p,primes[j]))
	return -1;
    }
    if(num>1)
      if(updChina(mp,a[i],num,num))
	return -1;
  }
  //for(auto t:mp)cout<<t<<",";cout<<endl;
  vector<ll> aa,mm;
  for(auto temp:mp){
    num=temp.Y.X; p=temp.Y.Y;
    aa.pb((num%p+p)%p);
    mm.pb(p);
  }
  //cout<<aa<<mm;
  return china2(aa,mm);
}

int main(){
  ios_base::sync_with_stdio(false);
  cout<<fixed<<setprecision(0);
  int i,j,k;
  findPrime();
  int n;
  cin>>n;
  vector<ll> a(n),m(n);
  long long te,mp;
  rep(i,n){
    cin>>te>>mp;
    a[i]=te; m[i]=mp;
  }
  ll re=china(a,m),lcm=1;
  rep(i,n)
    lcm*=m[i]/__gcd(lcm,m[i]);
  //  (re+=lcm)%=lcm;
  if(!re)re+=lcm;
  ll d=1e18;
  long long res=re/d,ult=re%d;
  if(res)
    cout<<res<<ult<<endl;
  else
    cout<<ult<<endl;
  return 0;
}
0