結果

問題 No.186 中華風 (Easy)
コンテスト
ユーザー koyumeishi
提出日時 2015-04-20 00:36:10
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 20 ms / 2,000 ms
コード長 1,949 bytes
コンパイル時間 811 ms
コンパイル使用メモリ 78,792 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 18:23:31
合計ジャッジ時間 1,651 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

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


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

	long long ans = -1;
	long long ub = lcm(y[0], y[1]);
	for(long long val = x[0]; val<ub; val += y[0]){
		if(val % y[1] == x[1]){
			long long ub_ = lcm(ub, y[2]);
			for(long long val_ = val; val_<ub_; val_ += ub){
				if(val_ % y[2] == x[2]){
					if(ans < 0) ans = val_;
					else ans = min(ans, val_);
				}
			}
		}
	}

	if(ans == 0) ans = lcm(y[0], lcm(y[1],y[2]));
	cout << ans << endl;

	return 0;
}
0