結果

問題 No.2748 Strange Clock
コンテスト
ユーザー ponjuice
提出日時 2024-04-19 13:43:43
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 2,531 bytes
コンパイル時間 2,457 ms
コンパイル使用メモリ 208,232 KB
最終ジャッジ日時 2025-02-21 03:17:54
ジャッジサーバーID
（参考情報） 		judge4 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 2
other AC * 28 TLE * 1 -- * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using ll = long long;

#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,avx2,tune=native")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")

ll ll_pow(ll a, ll n) {
    ll res = 1;
    while (n > 0) {
        if (n & 1) res *= a;
        a = a * a;
        n >>= 1;
    }
    return res;
}

ll n, m;

array<int, 15> itoar(ll i) {
    array<int,15>res;
    for(int j = 0; j < 15; j++) {
        res[14 - j] = i % 3;
        i /= 3;
    }
    return res;
}

ll artoi(array<int, 15> a, ll e) {
    ll res = 0;
    for(int j = 0; j < 15; j++) {
        res *= e;
        res += a[j];
    }
    return res;
}

pair<ll, ll> extgcd(ll a, ll b) {
    ll s = a, xs = 1, ys = 0, t = b, xt = 0, yt = 1;
    while (true) {
        ll temp = s / t;
        if(temp * t == s)break;
        ll u = s - t * temp;
        ll xu = xs - xt * temp;
        ll yu = ys - yt * temp;

        s = t;
        xs = xt;
        ys = yt;
        t = u;
        xt = xu;
        yt = yu;
    }

    return {xt, yt};
}

ll extgcd(ll a, ll b, ll &p, ll &q) {  
    if (b == 0) { 
        p = 1;
        q = 0;
        return a;
    }  
    ll d = extgcd(b, a % b, q, p);  
    q -= a / b * p;  
    return d;  
}

ll crt(ll a, ll mod1, ll b, ll mod2) {
    auto [p, q] = extgcd(mod1, mod2);
    ll mod = (mod1 * mod2);
    return ((a + mod1 * ((b - a) * p % mod2)) % mod + mod) % mod;
}

int main(){
    cin >> n >> m;

    const ll pow3 = ll_pow(3, n);
    const ll pow4 = ll_pow(4, n);
    const ll pow6 = ll_pow(6, n);
    const ll pow12 = ll_pow(12, n);

    vector<pair<ll,ll>> res;

    for(int i = 0; i < pow3; i += 3) {
        // i
        array<int,15> ret = itoar(i);
        ll i3 = i;
        ll i4 = artoi(ret, 4);
        ll i6 = artoi(ret, 6);

        ll t = crt(i3, pow3, i4, pow4);

        i6 = ((i6 - t) % pow6 + pow6) % pow6;

        res.emplace_back(i6, t);

        // i + 1
        res.emplace_back(i6, t + 1);

        // i + 2
        res.emplace_back(i6, t + 2);
    }

    sort(res.begin(), res.end());
    ll ans = 0;

    ll nw = -1, lt = -1;
    for(int i = 0; i < (int)res.size(); i++) {
        if(nw != res[i].first){
            nw = res[i].first;
            lt = res[i].second + pow12;
        }

        if(i + 1 == res.size() || res[i].first != res[i + 1].first) {
            ans += (lt - res[i].second > m) ? 1 : 0;
        }else{
            ans += (res[i+1].second - res[i].second > m) ? 1 : 0;
        }
    }

    cout << ans << endl;
}
0