結果

問題 No.1528 Not 1
コンテスト
ユーザー funera1
提出日時 2021-06-05 03:46:12
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,091 bytes
コンパイル時間 1,787 ms
コンパイル使用メモリ 178,460 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-20 19:26:24
合計ジャッジ時間 3,607 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 18 WA * 1
権限があれば一括ダウンロードができます

ソースコード

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