結果

問題 No.187 中華風 (Hard)
コンテスト
ユーザー koyumeishi
提出日時 2015-04-20 03:16:11
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 2,509 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 677 ms
コンパイル使用メモリ 80,208 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-07-04 16:46:04
合計ジャッジ時間 5,900 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 22 WA * 3
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <vector>
#include <cstdio>
#include <sstream>
#include <map>
#include <string>
#include <algorithm>
#include <queue>
#include <cmath>
#include <set>
using namespace std;

long long gcd(long long a, long long b){
	if(b==0) return a;
	return gcd(b, a%b);
}

long long lcm(long long a, long long b){
	if(a<b) swap(a,b);
	if(b==1) return a;
	return a * (b/gcd(a,b));
}


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

long long mod_inverse(long long a, long long m){
	long long x,y;
	extgcd(a,m,x,y);
	return (m+x%m)%m;
}


// return x where x mod p_i = value_i
// v[i] = pair{ value_i , p_i }
long long chinese_remainder_theorem(const vector<pair<long long,long long>>& v){
	long long M = 1;
	for(int i=0; i<v.size(); i++){
		M = lcm(M,v[i].second);
	}

	long long ret = 0;
	for(int i=0; i<v.size(); i++){
		long long M_i = (M/v[i].second);
		ret += (v[i].first * M_i * mod_inverse( M_i%v[i].second, v[i].second ))%M;
		ret = ret%M;
	}
	return ret%M;
}

// return x where x mod p_i = value_i
// v[i] = pair{ value_i , p_i }
long long chinese_remainder_theorem(const vector<long long>& value, const vector<long long>& mod){
	vector<pair<long long,long long>> v(value.size());
	for(int i=0; i<v.size(); i++){
		v[i] = {value[i], mod[i]};
	}
	return chinese_remainder_theorem(v);
}

#define MOD 1000000007

int main(){
	int N = 3;
	cin >> N;
	vector<long long> x(N), y(N);
	for(int i=0; i<N; i++){
		cin >> x[i] >> y[i];
	}

	bool valid = true;

	for(int i=0; i<N; i++){
		for(int j=i+1; j<N; j++){
			if(i == j) continue;
			long long g = gcd(y[i], y[j]);

			if( x[i]%g != x[j]%g ) valid = false;

			if(g != 1){
				y[i] /= g; y[j] /= g;
				long long g_ = gcd(y[i], g);
				while(g_ != 1){
					y[i] *= g_;
					g /= g_;
					g_ = gcd(y[i], g);
				}
				y[j] *= g;

				x[i] %= y[i];
				x[j] %= y[j];
			}
		}
	}

	if(!valid){
		cout << -1 << endl;
		return 0;
	}

	vector<long long> z(N);
	for(int i=0; i<N; i++){
		z[i] = x[i];
		for(int j=0; j<i; j++){
			z[i] = mod_inverse(y[j], y[i])%y[i] * (z[i] - z[j])%y[i];

			z[i] = (z[i]+y[i])%y[i];
		}
	}

	long long ans = 0;
	long long tmp = 1;
	for(int i=0; i<N; i++){
		ans = (ans + z[i] * tmp)%MOD;
		tmp = (tmp * y[i])%MOD;
	}

	if(ans == 0){
		long long lcm_ = 1;
		for(int i=0; i<N; i++){
			lcm_ = lcm(lcm_, y[i])%MOD;
		}
		ans = lcm_/*%MOD*/;
	}

	cout << ans << endl;
	return 0;
}
0