結果
問題 | No.187 中華風 (Hard) |
ユーザー | ふっぴー |
提出日時 | 2018-06-09 16:12:47 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 3,027 bytes |
コンパイル時間 | 1,647 ms |
コンパイル使用メモリ | 170,604 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-30 12:36:40 |
合計ジャッジ時間 | 2,472 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | AC | 3 ms
5,376 KB |
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 | WA | - |
testcase_11 | WA | - |
testcase_12 | WA | - |
testcase_13 | WA | - |
testcase_14 | WA | - |
testcase_15 | WA | - |
testcase_16 | WA | - |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 2 ms
5,376 KB |
testcase_20 | WA | - |
testcase_21 | AC | 2 ms
5,376 KB |
testcase_22 | WA | - |
testcase_23 | AC | 2 ms
5,376 KB |
testcase_24 | WA | - |
ソースコード
#include "bits/stdc++.h" using namespace std; #define DEBUG(x) cerr<<#x<<": "<<x<<endl; #define DEBUG_VEC(v) cerr<<#v<<":";for(int i=0;i<v.size();i++) cerr<<" "<<v[i]; cerr<<endl typedef long long ll; #define vi vector<int> #define vl vector<ll> #define vii vector< vector<int> > #define vll vector< vector<ll> > #define vs vector<string> #define pii pair<int,int> #define pis pair<int,string> #define psi pair<string,int> #define pll pair<ll,ll> template<class S, class T> pair<S, T> operator+(const pair<S, T> &s, const pair<S, T> &t) { return pair<S, T>(s.first + t.first, s.second + t.second); } template<class S, class T> pair<S, T> operator-(const pair<S, T> &s, const pair<S, T> &t) { return pair<S, T>(s.first - t.first, s.second - t.second); } template<class S, class T> ostream& operator<<(ostream& os, pair<S, T> p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define fi first #define se second #define rep(i,n) for(ll i=0;i<(ll)(n);i++) #define rep1(i,n) for(ll i=1;i<=(ll)(n);i++) #define rrep(i,n) for(ll i=(ll)(n)-1;i>=0;i--) #define rrep1(i,n) for(ll i=(ll)(n);i>0;i--) #define REP(i,a,b) for(ll i=(ll)a;i<(ll)b;i++) #define in(x, a, b) (a <= x && x < b) #define all(v) v.begin(),v.end() #define UNIQUE(v) v.erase(unique(all(v)), v.end()) template<class T> bool chmax(T &a, const T &b) { if (a<b) { a = b; return 1; } return 0; } template<class T> bool chmin(T &a, const T &b) { if (a>b) { a = b; return 1; } return 0; } const ll inf = 1000000001; const ll INF = 2e18; const ll MOD = 1000000007; const double pi = 3.14159265358979323846; #define Sp(p) cout<<setprecision(15)<< fixed<<p<<endl; int dx[4] = { 1,0, -1,0 }, dy[4] = { 0,1,0,-1 }; int dx2[8] = { 1,1,0,-1,-1,-1,0,1 }, dy2[8] = { 0,1,1,1,0,-1,-1,-1 }; #define fio() cin.tie(0); ios::sync_with_stdio(false); //ax + by = gcd(a, b) となるx, yをもとめる ll extgcd(ll a, ll b, ll& x, ll& y) { ll d = a; if (b != 0) { d = extgcd(b, a%b, y, x); y -= (a / b)*x; } else { x = 1; y = 0; } return d; } //中国剰余定理 //am[i].first:a_i, am[i].second:m_i //return (x, y) y = lcm(m1, m2, ..., m_n), x%m_i = a_i typedef __int128 lll; pll crt(vector<pll> am) { rep(i, am.size() - 1) { ll g, x, y, z; am[i].first = (am[i].first % am[i].second + am[i].second) % am[i].second; am[i + 1].first = (am[i + 1].first % am[i + 1].second + am[i + 1].second) % am[i + 1].second; g = extgcd(am[i].second, am[i + 1].second, x, y); if (am[i].first%g != am[i].second%g) return pll(-1, 0); //解なし lll lcm = (lll)am[i].second*(am[i + 1].second / g); if (lcm < am[i].second) return pll(-2, 0); //オーバーフロー ll mo = am[i + 1].second = am[i].second * am[i + 1].second; lll v = am[i].first + ((lll)(am[i].second / g) * x % mo * (am[i + 1].first - am[i].first + mo) % mo); am[i + 1].first = (v%mo + mo) % mo; } return am.back(); } int main() { int n; cin >> n; vector<pll> am(n); rep(i, n) { cin >> am[i].first >> am[i].second; } auto ans = crt(am); cout << ans.first % MOD << endl; }