結果

問題 No.187 中華風 (Hard)
ユーザー n_vipn_vip
提出日時 2015-05-24 19:02:52
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 1,949 ms / 3,000 ms
コード長 4,121 bytes
コンパイル時間 940 ms
コンパイル使用メモリ 111,152 KB
実行使用メモリ 8,236 KB
最終ジャッジ日時 2023-09-20 11:24:43
合計ジャッジ時間 20,978 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 17 ms
8,068 KB
testcase_01 AC 17 ms
8,068 KB
testcase_02 AC 982 ms
8,084 KB
testcase_03 AC 985 ms
7,952 KB
testcase_04 AC 999 ms
8,156 KB
testcase_05 AC 1,004 ms
7,980 KB
testcase_06 AC 1,003 ms
8,000 KB
testcase_07 AC 998 ms
8,156 KB
testcase_08 AC 976 ms
8,116 KB
testcase_09 AC 974 ms
7,920 KB
testcase_10 AC 975 ms
8,024 KB
testcase_11 AC 1,000 ms
7,976 KB
testcase_12 AC 999 ms
7,968 KB
testcase_13 AC 948 ms
8,152 KB
testcase_14 AC 948 ms
8,212 KB
testcase_15 AC 1,947 ms
8,008 KB
testcase_16 AC 1,949 ms
8,144 KB
testcase_17 AC 16 ms
8,236 KB
testcase_18 AC 17 ms
8,088 KB
testcase_19 AC 16 ms
7,944 KB
testcase_20 AC 872 ms
8,180 KB
testcase_21 AC 17 ms
8,136 KB
testcase_22 AC 1,003 ms
8,000 KB
testcase_23 AC 21 ms
7,892 KB
testcase_24 AC 17 ms
7,940 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 long long 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 ll MOD=1e9+7;

ll modpow(ll r,ll n,ll m=MOD){
  ll re=1,d=r%m;
  for(;n;n/=2){
    if(n&1)(re*=d)%=m;
    (d*=d)%=m;
  }
  return re;
}

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;
}

//m[i]=p[i]^t : 逆元求める用 modpow書き換える
ll garner(vector<ll> a,vector<ll> m,vector<ll> p,int f,ll mod=MOD){
  int n=a.size();
  m.pb(mod); p.pb(mod);
  vector<ll> d(n+1,1),sum(n+1);
  rep(i,n){
    //    sum[i]+x*d[i]==a[i] mod m[i];
    ll x=(a[i]-sum[i])%m[i]*modpow(d[i],m[i]/p[i]*(p[i]-1)-1,m[i])%m[i];
    //cout<<a[i]<<"-"<<sum[i]<<","<<d[i]<<":";//d[i]おかしい
    if(x<0) x+=m[i];
    //cout<<x<<",";cout<<"["<<d[i]<<"*"<<modpow(d[i],m[i]/p[i]*(p[i]-1)-1,m[i])<<"="<<d[i]*modpow(d[i],m[i]/p[i]*(p[i]-1)-1,m[i])%m[i]<<endl;
    if(f && x==0 && i==0)x=m[i];
    for(int j=i+1;j<=n;++j){
      (sum[j]+=x%m[j]*d[j]%m[j])%=m[j];
      (d[j]*=m[i])%=m[j];
    }
    //cout<<sum<<d;
    //cout<<sum.back()<<endl;
    //rep(j,n)cout<<sum.back()%m[j]<<",";cout<<endl;
  }
  return (sum.back()%MOD+MOD)%MOD;
}

//中国剰余定理
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,int f=0){
  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,pp;
  for(auto temp:mp){
    num=temp.Y.X; p=temp.Y.Y;
    aa.pb((num%p+p)%p);
    mm.pb(p);
    pp.pb(temp.X);
  }
  //cout<<aa<<mm;
  return garner(aa,mm,pp,f);
}

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;
  //cout<<re<<endl;
  rep(i,n)
    lcm*=m[i]/__gcd(lcm,m[i]);
  //  (re+=lcm)%=lcm;
  if(!re)re=china(a,m,1);
  ll d=1e18;
  long long res=re/d,ult=re%d;
  if(res)
    cout<<res<<ult<<endl;
  else
    cout<<ult<<endl;
  return 0;
}
0