結果

問題 No.186 中華風 (Easy)
ユーザー snteasntea
提出日時 2017-03-09 01:37:03
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 5,819 bytes
コンパイル時間 2,161 ms
コンパイル使用メモリ 196,116 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 18:30:43
合計ジャッジ時間 2,859 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 1 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 1 ms
5,376 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 1 ms
5,376 KB
testcase_13 AC 1 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 1 ms
5,376 KB
testcase_19 AC 2 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 2 ms
5,376 KB
testcase_22 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:11:21: warning: overflow in conversion from 'double' to 'int' changes value from '1.0e+18' to '2147483647' [-Woverflow]
   11 | const int mod = 1e18+7;
      |                 ~~~~^~

ソースコード

diff # 

#ifdef LOCAL111
	#define _GLIBCXX_DEBUG
#else
	#define NDEBUG
#endif
#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
const int INF = 1e9;
using namespace std;
template<typename T, typename U> ostream& operator<< (ostream& os, const pair<T,U>& p) { cout << '(' << p.first << ' ' << p.second << ')'; return os; }
const int mod = 1e18+7;

#define endl '\n'
#define ALL(a)  (a).begin(),(a).end()
#define SZ(a) int((a).size())
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define RFOR(i,a,b) for (int i=(b)-1;i>=(a);i--)
#define REP(i,n)  FOR(i,0,n)
#define RREP(i,n) for (int i=(n)-1;i>=0;i--)
#ifdef LOCAL111
	#define DEBUG(x) cout<<#x<<": "<<(x)<<endl
	template<typename T> void dpite(T a, T b){ for(T ite = a; ite != b; ite++) cout << (ite == a ? "" : " ") << *ite; cout << endl;}
#else
	#define DEBUG(x) true
	template<typename T> void dpite(T a, T b){ return; }
#endif
#define F first
#define S second
#define SNP string::npos
#define WRC(hoge) cout << "Case #" << (hoge)+1 << ": "
template<typename T> void pite(T a, T b){ for(T ite = a; ite != b; ite++) cout << (ite == a ? "" : " ") << *ite; cout << endl;}
template<typename T> bool chmax(T& a, T b){if(a < b){a = b; return true;} return false;}
template<typename T> bool chmin(T& a, T b){if(a > b){a = b; return true;} return false;}

template <typename T>
long long gcd(T x, T y){
	return y==0 ? x : gcd(y, x%y);
}

template <typename T>
long long lcm(T x, T y){
	return x/gcd(x,y)*y;
}

using Integer = long long;

//res.first*a+res.second*b == 1 となるresを返す (a,bは互いに素)
pair<Integer, Integer> extgcd(Integer a,Integer b)
{
	if(b==1){
		return pair<Integer, Integer>(0,1);
	}
	pair<Integer, Integer> t=extgcd(b,a%b);
	return pair<Integer, Integer>(t.second,t.first-a/b*t.second);
}

//modの逆元を返す
Integer inverse(Integer a,Integer modl)
{
	return (extgcd(modl,a).second+modl)%modl;
}


// x = second (mod first) 0以上の最小解
Integer garner_gen(const vector<pair<Integer, Integer>>& ex) {
	int n = ex.size();
	Integer res = 0;
	Integer k = 1;
	for(int i = 0; i < n; ++i) {
		Integer x, m;
		tie(m,x) = ex[i];
		x = (x%m+m)%m;
		Integer g = gcd(k,m);
		if((x-res)%g != 0) return -1;
		Integer inv = inverse(k/g,m/g);
		Integer v = ((x-res)/g%m+m)%m*inv%m;
		Integer nk = k/g*m;
		res = (res+v*k%nk)%nk;
		k = k/g*m;
	}
	return res;
}

Integer garner_gen(const vector<Integer>& x, const vector<Integer>& mod){
	int n = x.size();
	vector<pair<Integer, Integer>> v(n);
	for(int i = 0; i < n; ++i) {
		v[i] = {mod[i],x[i]};
	}
	return garner_gen(v);
}


//あんま検証してないよ
Integer garner(const vector<pair<Integer, Integer>>& ex) {
	int n = ex.size();
	Integer res = 0;
	Integer k = 1;
	for(int i = 0; i < n; ++i) {
		Integer x, m;
		tie(m,x) = ex[i];
		// x = (x%m+m)%m;
		Integer inv = inverse(k,m);
		Integer v = ((x-res)%m+m)%m*inv%m;
		res = (res+v*k);
		k = k*m;
	}
	return res;
}

Integer garner(const vector<Integer>& x, const vector<Integer>& mod){
	int n = x.size();
	vector<pair<Integer, Integer>> v(n);
	for(int i = 0; i < n; ++i) {
		v[i] = {mod[i],x[i]};
	}
	return garner(v);
}

Integer garner_mod(vector<pair<Integer, Integer>> ex, Integer mod) {
	int n = ex.size();
	Integer res = 0;
	// Integer k = 1;
	vector<pair<Integer,Integer>> memo(n);
	for(int i = 0; i < n; ++i) {
		Integer x, m;
		ex[i].S %= ex[i].F;
		tie(m,x) = ex[i];
		DEBUG(m); DEBUG(x);
		// x = (x%m+m)%m;
		Integer k = 1;
		Integer kk = 1;
		Integer ress = 0;
		for(int j = 0; j < i; ++j) {
			k = k*(ex[j].F%m)%m;
			ress = (ress+memo[j].F*kk%m)%m;
			kk = kk*(ex[j].F%mod)%mod;
		}
		Integer inv = inverse(k,m);
		Integer v = (x-ress+m)*inv%m;
		DEBUG(v); DEBUG(k);
		res = (res+v*kk%mod)%mod;
		memo[i] = {v,kk};
		DEBUG(res);
	}
	return res;
}

Integer garner_mod(const vector<Integer>& x, const vector<Integer>& mod, Integer modi){
	int n = x.size();
	vector<pair<Integer, Integer>> v(n);
	for(int i = 0; i < n; ++i) {
		v[i] = {mod[i],x[i]};
	}
	return garner_mod(v,modi);
}

template<typename T>
unordered_map<T,int> primeFactorize(T x){
	unordered_map<T,int> res;
	for(T i = 2; i*i <= x; i++){
		while(x%i == 0){
			res[i]++;
			x /= i;
		}
	}
	if(x != 1) res[x]++;
	return res;
}

typedef long long int LL;
typedef unsigned long long ULL;
typedef pair<int,int> P;

void ios_init(){
	//cout.setf(ios::fixed);
	//cout.precision(12);
#ifdef LOCAL111
	return;
#endif
	ios::sync_with_stdio(false); cin.tie(0);	
}

int main()
{
	ios_init();
	int n = 3;
	// while(cin >> n){
	{
		using P = pair<Integer,Integer>;
		vector<Integer> x(n),y(n);
		REP(i,n) cin >> x[i] >> y[i];
		
		vector<unordered_map<LL,int>> facs(n);
		unordered_map<LL,vector<P>> mp;
		REP(i,n){
			facs[i] = primeFactorize(y[i]);
			for(auto&& e : facs[i]) {
				mp[e.F].emplace_back(e.S,i);
			}
		}
		bool f = true;
		for(auto&& e : mp) {
			DEBUG(e.F);
			sort(ALL(e.S));
			int mod = pow(e.F,e.S.back().F);
			int val = x[e.S.back().S];
			REP(i,SZ(e.S)-1){
				DEBUG(y[e.S[i].S]);
				int modt = pow(e.F,e.S[i].F);
				int valt = x[e.S[i].S]%modt;
				if(val%modt != valt){
					DEBUG(i);
					f = false;
					// break;
				}
				y[e.S[i].S] /= modt;
				x[e.S[i].S] %= y[e.S[i].S];
			}
		}
		dpite(ALL(x)); dpite(ALL(y));
		
		if(!f){
			cout << -1 << endl;
		}else{
			bool zf = true;
			for(int i = 0; i < n; ++i) {
				if(x[i] != 0) zf = false;
			}
			LL ans = 1;
			if(zf){
				for(int i = 0; i < n; ++i) {
					ans *= y[i];
					// ans %= mod;
				}
			}else{
				vector<P> ex(n);
				REP(i,n) ex[i] = {y[i],x[i]};
				sort(ALL(ex));
				ex.erase(unique(ALL(ex)),ex.end());
				if(binary_search(ALL(ex),P({1,0}))){
					ex.erase(lower_bound(ALL(ex),P(0,1)));
				}
				if(binary_search(ALL(ex),P({1,1}))){
					ex.erase(lower_bound(ALL(ex),P(1,1)));
				}
				dpite(ALL(ex));
				ans = garner(ex);
			}
			cout << ans << endl;
		}
	}
	return 0;
}
0