結果

問題 No.1746 Sqrt Integer Segments
ユーザー HIcoderHIcoder
提出日時 2023-07-23 21:31:49
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 709 ms / 2,000 ms
コード長 1,708 bytes
コンパイル時間 1,442 ms
コンパイル使用メモリ 123,528 KB
実行使用メモリ 22,528 KB
最終ジャッジ日時 2024-09-25 00:56:34
合計ジャッジ時間 13,413 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 470 ms
17,280 KB
testcase_03 AC 255 ms
11,264 KB
testcase_04 AC 342 ms
13,696 KB
testcase_05 AC 664 ms
22,272 KB
testcase_06 AC 205 ms
9,856 KB
testcase_07 AC 67 ms
6,944 KB
testcase_08 AC 461 ms
16,896 KB
testcase_09 AC 687 ms
22,272 KB
testcase_10 AC 187 ms
9,216 KB
testcase_11 AC 605 ms
20,608 KB
testcase_12 AC 677 ms
22,528 KB
testcase_13 AC 692 ms
22,528 KB
testcase_14 AC 693 ms
22,400 KB
testcase_15 AC 709 ms
22,400 KB
testcase_16 AC 691 ms
22,400 KB
testcase_17 AC 404 ms
15,232 KB
testcase_18 AC 36 ms
7,936 KB
testcase_19 AC 37 ms
8,064 KB
testcase_20 AC 37 ms
8,064 KB
testcase_21 AC 195 ms
6,944 KB
testcase_22 AC 192 ms
6,944 KB
testcase_23 AC 96 ms
6,940 KB
testcase_24 AC 391 ms
8,064 KB
testcase_25 AC 392 ms
8,064 KB
testcase_26 AC 391 ms
7,936 KB
testcase_27 AC 159 ms
8,064 KB
testcase_28 AC 160 ms
8,064 KB
testcase_29 AC 159 ms
8,064 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

/*
ハッシュ化
*/
#include<iostream>
#include<set>
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<map>
#include<numeric>
#include<queue>
#include<cmath>
#include<random>
//#include<fstream>

using namespace std;
typedef long long ll;
const ll INF=1LL<<60;
typedef pair<ll,ll> P;
typedef pair<int,P> PP;
const ll MOD=998244353;

vector<pair<ll,ll>> prime_factor(ll n){
    vector<pair<ll,ll>> res;

    for(ll i=2;i*i<=n;i++){
        if(n%i==0){
            pair<ll,ll> tmp;
            tmp.first=i;
            tmp.second=0;

            while(n%i==0){
                tmp.second++;
                n/=i;
            }

            res.push_back(tmp);
        }
    }

    if(n>1){
        res.emplace_back(n,1);
        n/=n;
    }

    return res;
}

int main(){
    int N;
    cin>>N;
    vector<ll> a(N);
    for(int i=0;i<N;i++){
        cin>>a[i];
    }

    std::random_device seed_gen;
    mt19937 engine(seed_gen());

    uniform_int_distribution dist(1LL,1LL<<60);//[1,1LL<<60]までの乱数

    vector<ll> b(N),c(N+1);
    map<ll,ll> mp;//mp[v]=vに対しての乱数

    for(int i=0;i<N;i++){
        auto c=prime_factor(a[i]);

        for(auto [v,num]:c){
            if(num%2==1){
                if(!mp.count(v)){
                    mp[v]=dist(engine);//ハッシュ
                }

                b[i]^=mp[v];
            }
        }
    }

    //累積和のxorバージョン

    for(int i=0;i<N;i++){
        c[i+1]=c[i]^b[i];
    }
    map<ll,ll> mp2;
    for(int i=0;i<N+1;i++){
        mp2[c[i]]++;
    }

    ll ans=0;
    for(auto [v,num]:mp2){
        //vが同じになる区間
        ans+=num*(num-1)/2;
    }

    cout<<ans<<endl;

}
0