結果

問題 No.144 エラトステネスのざる
ユーザー isofuraBonjourisofuraBonjour
提出日時 2015-02-06 10:28:03
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 2,240 bytes
コンパイル時間 1,267 ms
コンパイル使用メモリ 74,408 KB
実行使用メモリ 32,936 KB
最終ジャッジ日時 2023-09-05 14:22:30
合計ジャッジ時間 2,767 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 WA -
testcase_09 AC 2 ms
4,376 KB
testcase_10 WA -
testcase_11 AC 1 ms
4,380 KB
testcase_12 AC 2 ms
4,376 KB
testcase_13 AC 87 ms
31,716 KB
testcase_14 AC 101 ms
31,560 KB
testcase_15 WA -
testcase_16 AC 99 ms
31,260 KB
testcase_17 AC 94 ms
32,132 KB
testcase_18 AC 94 ms
32,272 KB
testcase_19 AC 103 ms
32,132 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <vector>
#include <string>
#include <map>
#include <set>
#include <cmath>

#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define RFOR(i,a,b) for (int i=(b)-1;i>=(a);i--)
#define REP(i,n) for (int i=0;i<(n);i++)
#define RREP(i,n) for (int i=(n)-1;i>=0;i--)
#define PB push_back
#define ALL(a) (a).begin(),(a).end()
#define PII pair<int,int>
#define VS vector<string>
#define VD vector<double>
#define VL vector<long>
#define V(type) vector<type>
#define DEBUG(x) cout<<#x<<": "<<x<<endl

using namespace std;


VD Eratosthenes(int n, double p, double& result_p){
    VD koho, dst;
    
    FOR(i, 2, n+1){
        koho.push_back(i);
    }
    VD backup_koho = koho;
    //    if ( p != 1.0 ){
    REP(i, koho.size()){
        dst.push_back(koho[i]);
        
        int j = 0;
        //        if ( koho[i] > 1 ){
        double piv = backup_koho[i];
        while ( i + j * piv <= koho.size() - 1 ) {
            if ( koho[i + j * piv] >= 1 ) {
                koho[i + j * piv] = 1 - p;
                ++j;
                continue;
            }
            koho[i + j * piv] *= ( 1 - p );
            ++j;
        }
    }
    //      }
    //        }
    //    } else {
    //        REP(i, koho.size()){
    //            dst.push_back(koho[i]);
    //            int j = 0;
    //            if ( koho[i] > 1 ){
    //                double piv = backup_koho[i];
    //                while ( i + j * piv <= koho.size() - 1 ) {
    //                    koho[i + j * piv] = 0;
    //                    ++j;
    //                }
    //                //      }
    //            }
    //        }
    //
    //    }
    return dst;
}

int main(int argc, const char * argv[]) {
    
    int N;
    double P;
    VD result;
    double result_p = 0;
    cin >> N >> P;
    result = Eratosthenes(N, P, result_p);
    
    //    cout << result_p << endl;
    
    double sum = 0;
    FOR(i, 0, result.size()){
        if ( result[i] >= 1 ){
            sum++;
        } else {
            sum += result[i];
        }
        
    }
    
    //    if ( P == 0 ){
    //        cout << result.size() << endl;
    //    } else {
    cout << sum << endl;
    //    }
    
    return 0;
}
0