結果

問題 No.1544 [Cherry 2nd Tune C] Synchroscope
ユーザー abc3
提出日時 2021-06-11 22:26:53
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 355 ms / 2,000 ms
コード長 2,335 bytes
コンパイル時間 2,246 ms
コンパイル使用メモリ 197,292 KB
最終ジャッジ日時 2025-01-22 06:13:03
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 48
権限があれば一括ダウンロードができます

ソースコード

diff # 

    #include <bits/stdc++.h>
    using namespace std;
    #include <math.h>
    #include <iomanip>
    #include <cstdint>
    #include <string>
    #include <sstream>
    template<class T> inline bool chmax(T& a, T b) { if (a <= b) { a = b; return 1; } return 0; }
    template<class T> inline bool chmin(T& a, T b) { if (a >= b) { a = b; return 1; } return 0; }
    #define rep(i,n) for (int i = 0; i < (n); ++i)
    typedef long long ll;
    using P=pair<ll,ll>;
    const int INF=1001001001;
    //const int mod=998244353;
    
    // 負の数にも対応した mod
    // 例えば -17 を 5 で割った余りは本当は 3 (-17 ≡ 3 (mod. 5))
    // しかし単に -17 % 5 では -2 になってしまう
    inline long long mod(long long a, long long m) {
        return (a % m + m) % m;
    }
    
    // 拡張 Euclid の互除法
    // ap + bq = gcd(a, b) となる (p, q) を求め、d = gcd(a, b) をリターンします
    long long extGcd(long long a, long long b, long long &p, long long &q) {  
        if (b == 0) { p = 1; q = 0; return a; }  
        long long d = extGcd(b, a%b, q, p);  
        q -= a/b * p;  
        return d;  
    }
    
    // 中国剰余定理
    // リターン値を (r, m) とすると解は x ≡ r (mod. m)
    // 解なしの場合は (0, -1) をリターン
    pair<long long, long long> ChineseRem(long long b1, long long m1, long long b2, long long m2) {
      long long p, q;
      long long d = extGcd(m1, m2, p, q); // p is inv of m1/d (mod. m2/d)
      if ((b2 - b1) % d != 0) return make_pair(0, -1);
      long long m = m1 * (m2/d); // lcm of (m1, m2)
      long long tmp = (b2 - b1) / d * p % (m2/d);
      long long r = mod(b1 + m1 * tmp, m);
      return make_pair(r, m);
    }

    void solve(){
        int n,m;
        cin>>n>>m;
        vector<int>a(n),b(m);
        rep(i,n){cin>>a[i];}
        rep(i,m){cin>>b[i];}
        ll mn=1e18;
        rep(i,n){
            rep(j,m){
                if(a[i]!=b[j]){continue;}
                auto t=ChineseRem(i,n,j,m);
                if(t.first==0&&t.second==-1){continue;}
                chmin(mn,t.first+1);
            }
        }
        if(mn==1e18){mn=-1;}
        cout<<mn<<endl;
    }
    int main(){
        ios_base::sync_with_stdio(false);
        cin.tie(NULL);
        solve();
        return 0;
    }
0