結果

問題 No.2748 Strange Clock
ユーザー hiikunZhiikunZ
提出日時 2024-04-20 16:59:26
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,541 bytes
コンパイル時間 4,690 ms
コンパイル使用メモリ 244,916 KB
実行使用メモリ 335,816 KB
最終ジャッジ日時 2024-10-12 12:51:02
合計ジャッジ時間 14,337 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 2 ms
5,248 KB
testcase_02 WA -
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 AC 2 ms
5,248 KB
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 AC 14 ms
8,484 KB
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 AC 462 ms
163,324 KB
testcase_28 WA -
testcase_29 AC 1,308 ms
335,668 KB
testcase_30 WA -
testcase_31 AC 126 ms
44,400 KB
testcase_32 AC 123 ms
44,336 KB
testcase_33 AC 397 ms
163,324 KB
testcase_34 WA -
testcase_35 AC 1,174 ms
335,692 KB
testcase_36 AC 1,168 ms
335,816 KB
testcase_37 AC 2 ms
6,816 KB
testcase_38 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <ext/pb_ds/assoc_container.hpp>
using namespace __gnu_pbds;
#pragma GCC optimize("Ofast,unroll-loops")
#include <bits/stdc++.h>
using namespace std;
//#include<atcoder/math>
//using namespace atcoder;
using ll = long long int;
using ull = unsigned long long int;
using ld = long double;
constexpr ll MAX = 2000000000000000000;
constexpr ld PI = 3.14159265358979;
constexpr ll MOD = 0;//2024948111;
ld dotorad(ld K){ return PI * K / 180.0; }
ld radtodo(ld K){ return K * 180.0 / PI; }
mt19937 mt;
void randinit(){ srand((unsigned)time(NULL));mt = mt19937(rand()); }
__int128_t modinv(__int128_t A,__int128_t M){
    //A*r+B*n=gcd(A,B)
    A %= M;
    if(A == 0 || __gcd(A,M) != 1){
        //cout << "Error modinv(" << A << "," << M << ")" << endl;
        return -1;
    }
    __int128_t B = M,U = 1,V = 0;
    while(B){
        __int128_t T = A / B;
        A -= T * B;
        swap(A,B);
        U -= T * V;
        swap(U,V);
    }
    U %= M;
    if(U < 0) U += M;
    return U;
}
int main(){
    ll N,M;
    cin >> N >> M;
    vector<ll> pow3(60,1),pow4(60,1),pow6(60,1),pow2(60,1);
    for(ll i = 1;i < 60;i++){
        pow2[i] = pow2[i - 1] * 2;
        pow3[i] = pow3[i - 1] * 3;
        pow4[i] = pow4[i - 1] * 4;
        pow6[i] = pow6[i - 1] * 6;
    }
    ll l = pow2[N * 2] * pow3[N];
    __int128_t t1 = (modinv(pow3[N],pow2[N * 2]) * pow3[N]) % l;
    __int128_t t2 = (modinv(pow2[N * 2],pow3[N]) * pow2[N * 2]) % l;
    ll ID = 0;
    vector<ll> memo2(pow3[N - 1],0);
    for(ll bb = 0;bb < pow3[N];bb += 3){
        ll c = bb;
        ll p4 = 0,p6 = 0;
        for(ll i = 0;i < N;i++){
            ll x = c % 3;
            c /= 3;
            p4 += x * pow4[i];
            p6 += x * pow6[i];
        }
        ll x = (bb * t2 + p4 * t1) % l;
        ll zurasi = p6 - x;
        zurasi %= pow6[N];
        if(zurasi < 0) zurasi += pow6[N];
        memo2[bb / 3] = (zurasi);
    }
    gp_hash_table<ll,ll> X;

    ll ans = 0;
    for(ll id = 0;id < pow3[N - 1];id++){
        ll bb = 3 * id;
        ll z = memo2[id];
        if(X.find(z) != X.end()){
            ll z1 = X[z];
            ll z2 = bb;
            if(z2 < z1) swap(z1,z2);
            ll d1 = (z2 - (z1 + 3)) % l;
            if(d1 >= M) ans++;
            ll d2 = (z1 - (z2 + 3)) % l;
            if(d2 < 0) d2 += l;
            if(d2 >= M) ans++;

            X.erase(z);
            continue;
        }
        else{
            X[z] = bb;
        }
    }
    for(auto [z,bb] : X){
        if(M <= l + 3) ans++;
    }

    cout << ans << endl;
}
0