結果

問題 No.1528 Not 1
ユーザー funera1funera1
提出日時 2021-06-05 03:46:12
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,091 bytes
コンパイル時間 1,909 ms
コンパイル使用メモリ 177,676 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-04-30 21:25:11
合計ジャッジ時間 3,646 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 13 ms
6,940 KB
testcase_02 AC 6 ms
6,940 KB
testcase_03 AC 7 ms
6,940 KB
testcase_04 AC 13 ms
6,940 KB
testcase_05 AC 4 ms
6,940 KB
testcase_06 AC 9 ms
6,944 KB
testcase_07 AC 14 ms
6,944 KB
testcase_08 AC 8 ms
6,944 KB
testcase_09 AC 5 ms
6,940 KB
testcase_10 AC 8 ms
6,944 KB
testcase_11 WA -
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 2 ms
6,944 KB
testcase_14 AC 2 ms
6,944 KB
testcase_15 AC 2 ms
6,944 KB
testcase_16 AC 2 ms
6,940 KB
testcase_17 AC 2 ms
6,940 KB
testcase_18 AC 15 ms
6,940 KB
testcase_19 AC 15 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
//#include <atcoder/all>
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#define rep(i,n)for(int i=0;(i)<(int)(n);i++)
#define REP(i,a,b)for(int i=(int)(a);(i)<=(int)(b);i++)
#define ALL(a)  (a).begin(),(a).end()
#define pb push_back
#define fi first
#define se second
#define sz(x) ((int)x.size())
 
using namespace std;
//using namespace atcoder;
using ld = long double;
using ll = long long;
using P = pair<ll, ll>;
 
template<typename T> bool chmin(T& a, const T& b) { if(a >= b){ a = b; return 1;} return 0; }
template<typename T> bool chmax(T& a, const T& b) { if(a <= b){ a = b; return 1;} return 0; }

const ll MOD = 1e9 + 7;
const ll INF = LLONG_MAX;

//エラトステネスの篩
vector<bool> IsPrime;

void sieve(size_t max){
    if(max+1 > IsPrime.size()){     // resizeで要素数が減らないように
        IsPrime.resize(max+1,true); // IsPrimeに必要な要素数を確保
    } 
    IsPrime[0] = false; // 0は素数ではない
    IsPrime[1] = false; // 1は素数ではない

    for(size_t i=2; i*i<=max; ++i) // 0からsqrt(max)まで調べる
        if(IsPrime[i]) // iが素数ならば
            for(size_t j=2; i*j<=max; ++j) // (max以下の)iの倍数は
                IsPrime[i*j] = false;      // 素数ではない
}

int main(){
    int n;
    cin >> n;
    vector<int> a;
    sieve(n + 1); 
    vector<int> p, used(n + 1, 0);
    int lim = (n + 1) / 2;
    rep(i, n + 1)if(IsPrime[i])p.pb(i);

    rep(i, sz(p)){
        if(sz(a) == lim)break;
        for(int j = p[i]; j <= n; j += p[i]){
            if(i + 1 < sz(p) && p[i] * p[i + 1] == j)continue;
            if(used[j])continue;
            a.pb(j);
            used[j] = 1;
            if(sz(a) == lim)break;
        }
        if(sz(a) == lim)break;

        if(i + 1 < sz(p) && (ll)p[i] * (ll)p[i + 1] <= n){
            a.pb(p[i] * p[i + 1]);
            used[p[i] * p[i + 1]] = 1;
        }
        else break;
    }
    if(sz(a) < lim)cout << -1 << endl;
    else for(auto ai : a)cout << ai << " " ;cout << endl;
}
0