結果

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

ソースコード

diff # 

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

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

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

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, ll p, ll q) {
    ll mod = (mod1 * mod2);
    ll res = (a + mod1 * ((b - a) * p % mod2)) % mod;
    if(res < 0) res += mod;
    return res;
}

int main(){
    ll n, m;
    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);
    ll p, q;
    extgcd(pow3, pow4, p, q);

    vector<pair<ll,ll>> res;
    res.reserve(pow3 / 3);

    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, p, q);

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

        res.emplace_back(i6, t);
    }

    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 + 2) > m) ? 1 : 0;
        }else{
            ans += (res[i+1].second - (res[i].second + 2) > m) ? 1 : 0;
        }
    }

    cout << ans << endl;
}
0