結果

問題 No.2798 Multiple Chain
ユーザー GOTKAKOGOTKAKO
提出日時 2024-06-29 00:03:32
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 2,425 bytes
コンパイル時間 2,449 ms
コンパイル使用メモリ 213,656 KB
実行使用メモリ 82,764 KB
最終ジャッジ日時 2024-06-29 00:03:42
合計ジャッジ時間 9,127 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 6 ms
6,944 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 26 ms
6,940 KB
testcase_06 AC 38 ms
8,080 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 2 ms
6,944 KB
testcase_09 AC 2 ms
6,944 KB
testcase_10 AC 2 ms
6,944 KB
testcase_11 AC 2 ms
6,940 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 2 ms
6,940 KB
testcase_14 AC 4 ms
6,944 KB
testcase_15 AC 5 ms
6,944 KB
testcase_16 AC 6 ms
6,944 KB
testcase_17 AC 7 ms
6,944 KB
testcase_18 AC 6 ms
6,944 KB
testcase_19 AC 7 ms
6,944 KB
testcase_20 AC 5 ms
6,940 KB
testcase_21 AC 4 ms
6,944 KB
testcase_22 AC 5 ms
6,944 KB
testcase_23 AC 5 ms
6,944 KB
testcase_24 AC 5 ms
6,944 KB
testcase_25 AC 5 ms
6,944 KB
testcase_26 AC 5 ms
6,940 KB
testcase_27 AC 3 ms
6,944 KB
testcase_28 AC 4 ms
6,940 KB
testcase_29 AC 4 ms
6,944 KB
testcase_30 AC 5 ms
6,940 KB
testcase_31 AC 4 ms
6,944 KB
testcase_32 AC 11 ms
6,940 KB
testcase_33 AC 4 ms
6,944 KB
testcase_34 AC 9 ms
6,944 KB
testcase_35 AC 4 ms
6,944 KB
testcase_36 AC 5 ms
6,940 KB
testcase_37 AC 6 ms
6,940 KB
testcase_38 AC 4 ms
6,944 KB
testcase_39 AC 3 ms
6,940 KB
testcase_40 AC 249 ms
22,220 KB
testcase_41 AC 1,095 ms
45,540 KB
testcase_42 TLE -
testcase_43 AC 1,074 ms
47,116 KB
testcase_44 AC 401 ms
24,572 KB
testcase_45 AC 1 ms
6,944 KB
testcase_46 AC 2 ms
6,940 KB
testcase_47 AC 2 ms
6,948 KB
testcase_48 AC 2 ms
6,944 KB
testcase_49 AC 2 ms
6,944 KB
testcase_50 AC 2 ms
6,944 KB
testcase_51 AC 5 ms
6,940 KB
testcase_52 AC 5 ms
6,944 KB
testcase_53 AC 3 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    long long N; cin >> N;
    long long memoN = N;
    vector<int> pf;
    {
        int e2 = 0,e3 = 0;
        while(N%2 == 0) e2++,N /= 2;
        while(N%3 == 0) e3++,N /= 3;
        if(e2 >= 2) pf.push_back(e2);
        if(e3 >= 2) pf.push_back(e3);
    }
    {
        long long i=5;
        while(i*i*i <= memoN){
            int e = 0;
            while(N%i == 0) e++,N /= i;
            if(e >= 2) pf.push_back(e);
            e = 0; i += 2;
            while(N%i == 0) e++,N /= i;
            if(e >= 2) pf.push_back(e);
            i += 4;
        }
    }
    long long low = 0,high = 2e9;
    while(high-low > 1){
        long long mid = (high+low)/2;
        if(mid*mid >= N) high = mid;
        else low = mid;
    }
    if(high*high == N){
        bool prime = true;
        for(int i=2; i*i<=high; i++){
            if(high%i == 0) prime = false;
            int e = 0;
            while(high%i == 0) e++,high /= i;
            pf.push_back(e*2);
        }
        if(high != 1) pf.push_back(2);
    }

    int answer = 0,n = pf.size();
    vector<int> multi(n+1,1);
    for(int i=1; i<=n; i++){
        int now = pf.at(i-1)+1;
        multi.at(i) = multi.at(i-1)*now;
    }
    swap(pf,multi); //まちがえた.

    unordered_map<long long,long long> M;
    auto dfs = [&](auto dfs,int left,int need) -> long long {
        long long &ret = M[1001001001LL*left+need];
        if(ret){answer += ret; return ret;}

        for(int i=0; i<n; i++){
            int le = left%pf.at(i+1)/pf.at(i),ne = need%pf.at(i+1)/pf.at(i);
            if(ne > le) return 0;
            left -= ne*pf.at(i);
        }
        answer++; ret++;

        int next = left;
        if(need) ret += dfs(dfs,left,need);

        while(true){
            int pos = 0;
            while(pos != n){
                int now = next%pf.at(pos+1)/pf.at(pos);
                if(now == 0){
                    next += left%pf.at(pos+1)/pf.at(pos)*pf.at(pos);
                    need -= left%pf.at(pos+1)/pf.at(pos)*pf.at(pos);
                    pos++;
                }
                else{next -= pf.at(pos),need += pf.at(pos); break;}
            }
            if(pos == n) break;
            ret += dfs(dfs,next,need);
        }
        return ret;
    };
    dfs(dfs,pf.at(n)-1,0);
    cout << answer << endl;
}
0