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

提出ソース

結果

問題	No.732 3PrimeCounting
ユーザー	kei
提出日時	2018-09-07 22:44:59
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	TLE
実行時間	-
コード長	3,564 bytes
コンパイル時間	2,327 ms
コンパイル使用メモリ	166,792 KB
実行使用メモリ	20,160 KB
最終ジャッジ日時	2023-08-19 21:12:47
合計ジャッジ時間	45,502 ms
ジャッジサーバーID （参考情報）	judge15 / judge13

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	8 ms 20,160 KB
testcase_01	AC	7 ms 9,480 KB
testcase_02	AC	7 ms 9,480 KB
testcase_03	AC	7 ms 9,624 KB
testcase_04	AC	7 ms 9,276 KB
testcase_05	AC	8 ms 9,396 KB
testcase_06	AC	7 ms 9,416 KB
testcase_07	AC	7 ms 9,536 KB
testcase_08	AC	7 ms 9,280 KB
testcase_09	AC	8 ms 9,300 KB
testcase_10	AC	7 ms 9,540 KB
testcase_11	AC	7 ms 9,448 KB
testcase_12	AC	7 ms 9,316 KB
testcase_13	AC	7 ms 9,300 KB
testcase_14	AC	7 ms 9,476 KB
testcase_15	AC	7 ms 9,296 KB
testcase_16	AC	7 ms 9,400 KB
testcase_17	AC	7 ms 9,276 KB
testcase_18	AC	7 ms 9,540 KB
testcase_19	AC	7 ms 9,400 KB
testcase_20	AC	8 ms 9,440 KB
testcase_21	AC	11 ms 9,304 KB
testcase_22	AC	10 ms 9,504 KB
testcase_23	AC	7 ms 9,404 KB
testcase_24	AC	7 ms 9,540 KB
testcase_25	AC	19 ms 9,500 KB
testcase_26	AC	13 ms 9,408 KB
testcase_27	AC	8 ms 9,520 KB
testcase_28	AC	7 ms 9,284 KB
testcase_29	AC	9 ms 9,544 KB
testcase_30	AC	10 ms 9,372 KB
testcase_31	AC	11 ms 9,592 KB
testcase_32	AC	16 ms 9,544 KB
testcase_33	AC	16 ms 9,276 KB
testcase_34	AC	15 ms 9,556 KB
testcase_35	AC	12 ms 9,536 KB
testcase_36	AC	12 ms 9,372 KB
testcase_37	AC	8 ms 9,536 KB
testcase_38	AC	7 ms 9,276 KB
testcase_39	AC	12 ms 9,316 KB
testcase_40	AC	11 ms 9,300 KB
testcase_41	AC	11 ms 9,604 KB
testcase_42	AC	12 ms 9,560 KB
testcase_43	AC	10 ms 9,400 KB
testcase_44	AC	10 ms 9,316 KB
testcase_45	AC	8 ms 9,328 KB
testcase_46	AC	8 ms 9,368 KB
testcase_47	AC	9 ms 9,496 KB
testcase_48	AC	22 ms 9,524 KB
testcase_49	AC	22 ms 9,484 KB
testcase_50	AC	11 ms 9,560 KB
testcase_51	AC	11 ms 9,316 KB
testcase_52	AC	8 ms 9,300 KB
testcase_53	AC	62 ms 9,504 KB
testcase_54	AC	1,396 ms 10,060 KB
testcase_55	AC	1,394 ms 9,884 KB
testcase_56	AC	1,514 ms 10,164 KB
testcase_57	AC	222 ms 9,704 KB
testcase_58	AC	223 ms 9,792 KB
testcase_59	AC	71 ms 9,416 KB
testcase_60	AC	354 ms 9,644 KB
testcase_61	AC	354 ms 9,680 KB
testcase_62	AC	1,885 ms 10,328 KB
testcase_63	AC	731 ms 9,640 KB
testcase_64	AC	477 ms 9,708 KB
testcase_65	AC	478 ms 9,640 KB
testcase_66	AC	9 ms 9,608 KB
testcase_67	AC	8 ms 9,276 KB
testcase_68	AC	1,613 ms 10,048 KB
testcase_69	AC	1,613 ms 10,060 KB
testcase_70	AC	1,426 ms 10,180 KB
testcase_71	AC	1,427 ms 9,888 KB
testcase_72	AC	679 ms 9,828 KB
testcase_73	TLE	-
testcase_74	TLE	-
testcase_75	AC	20 ms 9,540 KB
testcase_76	AC	1,571 ms 10,048 KB
testcase_77	AC	245 ms 9,704 KB
testcase_78	AC	2,409 ms 10,152 KB
testcase_79	AC	1,641 ms 9,892 KB
testcase_80	AC	2,120 ms 10,204 KB
testcase_81	AC	1,442 ms 9,972 KB
testcase_82	AC	9 ms 9,372 KB
testcase_83	AC	234 ms 9,928 KB
testcase_84	AC	268 ms 9,784 KB
testcase_85	AC	1,187 ms 9,904 KB
testcase_86	AC	2,179 ms 10,392 KB
testcase_87	TLE	-
testcase_88	--	-

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

ソースコード

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

結果

テストケース

ソースコード