#283592 (C++14) No.732 3PrimeCounting

提出ソース

結果

問題	No.732 3PrimeCounting
ユーザー	kei
提出日時	2018-09-07 22:55:44
言語	C++14 (gcc 13.2.0 + boost 1.83.0)
結果	AC
実行時間	2,808 ms / 3,000 ms
コード長	3,957 bytes
コンパイル時間	2,368 ms
コンパイル使用メモリ	168,252 KB
実行使用メモリ	10,932 KB
最終ジャッジ日時	2023-08-19 21:50:13
合計ジャッジ時間	35,229 ms
ジャッジサーバーID （参考情報）	judge12 / judge15

このコードへのチャレンジ（β）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	8 ms 9,524 KB
testcase_01	AC	8 ms 9,272 KB
testcase_02	AC	8 ms 9,516 KB
testcase_03	AC	8 ms 9,408 KB
testcase_04	AC	8 ms 9,276 KB
testcase_05	AC	8 ms 9,316 KB
testcase_06	AC	8 ms 9,268 KB
testcase_07	AC	8 ms 9,404 KB
testcase_08	AC	8 ms 9,276 KB
testcase_09	AC	8 ms 9,520 KB
testcase_10	AC	9 ms 9,404 KB
testcase_11	AC	9 ms 9,276 KB
testcase_12	AC	8 ms 9,524 KB
testcase_13	AC	8 ms 9,280 KB
testcase_14	AC	8 ms 9,480 KB
testcase_15	AC	9 ms 9,588 KB
testcase_16	AC	8 ms 9,272 KB
testcase_17	AC	9 ms 9,404 KB
testcase_18	AC	9 ms 9,396 KB
testcase_19	AC	8 ms 9,404 KB
testcase_20	AC	9 ms 9,548 KB
testcase_21	AC	10 ms 9,464 KB
testcase_22	AC	11 ms 9,276 KB
testcase_23	AC	9 ms 9,524 KB
testcase_24	AC	8 ms 9,280 KB
testcase_25	AC	16 ms 9,392 KB
testcase_26	AC	12 ms 9,596 KB
testcase_27	AC	9 ms 9,376 KB
testcase_28	AC	9 ms 9,316 KB
testcase_29	AC	10 ms 9,480 KB
testcase_30	AC	10 ms 9,596 KB
testcase_31	AC	11 ms 9,472 KB
testcase_32	AC	13 ms 9,496 KB
testcase_33	AC	14 ms 9,468 KB
testcase_34	AC	13 ms 9,316 KB
testcase_35	AC	12 ms 9,396 KB
testcase_36	AC	11 ms 9,480 KB
testcase_37	AC	9 ms 9,312 KB
testcase_38	AC	9 ms 9,452 KB
testcase_39	AC	11 ms 9,416 KB
testcase_40	AC	11 ms 9,524 KB
testcase_41	AC	11 ms 9,456 KB
testcase_42	AC	11 ms 9,424 KB
testcase_43	AC	10 ms 9,472 KB
testcase_44	AC	11 ms 9,320 KB
testcase_45	AC	8 ms 9,428 KB
testcase_46	AC	8 ms 9,376 KB
testcase_47	AC	9 ms 9,268 KB
testcase_48	AC	17 ms 9,320 KB
testcase_49	AC	17 ms 9,268 KB
testcase_50	AC	11 ms 9,268 KB
testcase_51	AC	10 ms 9,404 KB
testcase_52	AC	9 ms 9,404 KB
testcase_53	AC	43 ms 9,308 KB
testcase_54	AC	902 ms 10,024 KB
testcase_55	AC	902 ms 10,028 KB
testcase_56	AC	902 ms 9,900 KB
testcase_57	AC	146 ms 9,764 KB
testcase_58	AC	147 ms 9,716 KB
testcase_59	AC	51 ms 9,540 KB
testcase_60	AC	233 ms 9,716 KB
testcase_61	AC	233 ms 9,716 KB
testcase_62	AC	1,214 ms 10,208 KB
testcase_63	AC	475 ms 9,640 KB
testcase_64	AC	312 ms 9,752 KB
testcase_65	AC	311 ms 9,680 KB
testcase_66	AC	8 ms 9,516 KB
testcase_67	AC	8 ms 9,372 KB
testcase_68	AC	1,040 ms 9,964 KB
testcase_69	AC	1,041 ms 9,880 KB
testcase_70	AC	918 ms 9,908 KB
testcase_71	AC	919 ms 9,876 KB
testcase_72	AC	442 ms 9,612 KB
testcase_73	AC	2,027 ms 10,300 KB
testcase_74	AC	2,029 ms 10,336 KB
testcase_75	AC	15 ms 9,372 KB
testcase_76	AC	1,012 ms 10,104 KB
testcase_77	AC	163 ms 9,620 KB
testcase_78	AC	1,569 ms 10,292 KB
testcase_79	AC	1,057 ms 10,052 KB
testcase_80	AC	1,365 ms 10,164 KB
testcase_81	AC	932 ms 10,120 KB
testcase_82	AC	9 ms 9,412 KB
testcase_83	AC	153 ms 9,756 KB
testcase_84	AC	177 ms 9,808 KB
testcase_85	AC	765 ms 10,008 KB
testcase_86	AC	1,405 ms 10,292 KB
testcase_87	AC	2,808 ms 10,808 KB
testcase_88	AC	2,804 ms 10,932 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include "bits/stdc++.h"
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
const int INF = 1e9;
const ll LINF = 1e18;
template<class S,class T> ostream& operator << (ostream& out,const pair<S,T>& o){ out << "(" << o.first << "," << o.second << ")"; return out; }
template<class T> ostream& operator << (ostream& out,const vector<T> V){ for(int i = 0; i < V.size(); i++){ out << V[i]; if(i!=V.size()-1) out << " ";} return out; }
template<class T> ostream& operator << (ostream& out,const vector<vector<T> > Mat){ for(int i = 0; i < Mat.size(); i++) { if(i != 0) out << endl; out << Mat[i];} return out; }
template<class S,class T> ostream& operator << (ostream& out,const map<S,T> mp){ out << "{ "; for(auto it = mp.begin(); it != mp.end(); it++){ out << it->first << ":" << it->second; if(mp.size()-1 != distance(mp.begin(),it)) out << ", "; } out << " }"; return out; }

/*
 <url:>
 問題文============================================================
 =================================================================
 解説=============================================================
 ================================================================
 */

const ll MAX_PRIME = 500050;
vector<int> primes;
vector<int> is_prime(MAX_PRIME + 1,true);
void init_primes(){
    is_prime[0] = is_prime[1] = false;
    for(int i = 2; i <= MAX_PRIME;i++){
        if(is_prime[i]){
            primes.push_back(i);
            for(int j = i*2; j <= MAX_PRIME; j+=i) is_prime[j] = false;
        }
    }
}

int dp[4][MAX_PRIME];
int dp1[4][MAX_PRIME];
int dp2[4][MAX_PRIME];

ll solve(){
    init_primes();
    ll res = 0;
    ll N; cin >> N;
    dp[0][0] = 1;
    
    for(int i = 0; i < primes.size();i++){
        if(N < primes[i]) continue;
        
        int j = 2;
        for(int k = 2*primes[i]; k >= 0;k--){
            if(dp[j][k] == 0) continue;
            dp[j+1][k+primes[i]] += dp[j][k];
        }
        j = 1;
        for(int k = primes[i]; k >= 0; k--){
            if(dp[j][k] == 0) continue;
            dp[j+1][k+primes[i]] += dp[j][k];
        }
        j = 0;
        dp[j+1][primes[i]] += dp[j][0];
//        for(int j = 2; j >= 0;j--){
//            for(int k = 3*primes[i]; k >= 0;k--){
//                if(dp[j][k] == 0) continue;
//                dp[j+1][k+primes[i]] += dp[j][k];
//            }
//        }
    }
    for(int i = 0; i <= 3*N;i++){
        if(is_prime[i]){
            res += dp[3][i];
        }
    }

    
    //    ll Lim = upper_bound(primes.begin(),primes.end(),N) - primes.begin();
    //    dp1[0][0] = 1;
    //    for(ll i = 0; i < Lim/2;i++){
    //        //if(N < primes[i]) continue;
    //        for(int j = 2; j >= 0;j--){
    //            for(int k = 3*primes[i]; k >= 0;k--){
    //                if(dp1[j][k] == 0) continue;
    //                dp1[j+1][k+primes[i]] += dp1[j][k];
    //            }
    //        }
    //    }
    //    cout << "ok" << endl;
    //
    //    dp2[0][0] = 1;
    //    for(ll i = Lim/2; i < Lim;i++){
    //        for(int j = 2; j >= 0;j--){
    //            for(int k = 3*primes[i]; k >= 0;k--){
    //                if(dp2[j][k] == 0) continue;
    //                dp2[j+1][k+primes[i]] += dp2[j][k];
    //            }
    //        }
    //    }
    //    cout << "ok2" << endl;
    //    for(ll k1 = 0; k1 <= 3*primes[Lim/2];k1++){
    //        for(ll k2 = 3*primes[Lim/2]; k2 <= 3*N;k2++){
    //            if(is_prime[k1+k2]){
    //                res += dp1[1][k1]*dp2[2][k2];
    //                res += dp1[2][k1]*dp2[1][k2];
    //            }
    //        }
    //    }
    //
    //    for(int i = 0; i <= 3*N;i++){
    //        if(is_prime[i]){
    //            res += dp1[3][i];
    //            res += dp2[3][i];
    //        }
    //    }
    //
    return res;
}
int main(void) {
    cin.tie(0); ios_base::sync_with_stdio(false);
    cout << solve() << endl;
    return 0;
}

yukicoder

結果

テストケース

ソースコード