結果

問題 No.719 Coprime
ユーザー mugen_1337mugen_1337
提出日時 2020-10-15 01:02:09
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 10 ms / 3,000 ms
コード長 2,170 bytes
コンパイル時間 1,920 ms
コンパイル使用メモリ 204,312 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-28 01:01:26
合計ジャッジ時間 4,339 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 2 ms
4,376 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 3 ms
4,376 KB
testcase_11 AC 2 ms
4,376 KB
testcase_12 AC 2 ms
4,376 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 2 ms
4,376 KB
testcase_15 AC 2 ms
4,376 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 2 ms
4,380 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 2 ms
4,376 KB
testcase_20 AC 3 ms
4,376 KB
testcase_21 AC 2 ms
4,376 KB
testcase_22 AC 2 ms
4,376 KB
testcase_23 AC 3 ms
4,376 KB
testcase_24 AC 2 ms
4,380 KB
testcase_25 AC 2 ms
4,376 KB
testcase_26 AC 2 ms
4,376 KB
testcase_27 AC 2 ms
4,380 KB
testcase_28 AC 2 ms
4,376 KB
testcase_29 AC 2 ms
4,376 KB
testcase_30 AC 2 ms
4,376 KB
testcase_31 AC 2 ms
4,380 KB
testcase_32 AC 2 ms
4,376 KB
testcase_33 AC 3 ms
4,380 KB
testcase_34 AC 2 ms
4,384 KB
testcase_35 AC 2 ms
4,380 KB
testcase_36 AC 2 ms
4,376 KB
testcase_37 AC 2 ms
4,380 KB
testcase_38 AC 2 ms
4,380 KB
testcase_39 AC 2 ms
4,380 KB
testcase_40 AC 2 ms
4,376 KB
testcase_41 AC 2 ms
4,376 KB
testcase_42 AC 3 ms
4,376 KB
testcase_43 AC 3 ms
4,380 KB
testcase_44 AC 2 ms
4,380 KB
testcase_45 AC 2 ms
4,380 KB
testcase_46 AC 3 ms
4,376 KB
testcase_47 AC 3 ms
4,376 KB
testcase_48 AC 3 ms
4,376 KB
testcase_49 AC 3 ms
4,380 KB
testcase_50 AC 3 ms
4,380 KB
testcase_51 AC 3 ms
4,380 KB
testcase_52 AC 5 ms
4,376 KB
testcase_53 AC 6 ms
4,376 KB
testcase_54 AC 8 ms
4,376 KB
testcase_55 AC 7 ms
4,376 KB
testcase_56 AC 8 ms
4,380 KB
testcase_57 AC 9 ms
4,376 KB
testcase_58 AC 10 ms
4,380 KB
testcase_59 AC 10 ms
4,380 KB
testcase_60 AC 10 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define ALL(x) begin(x),end(x)
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define mod 1000000007
using ll=long long;
const int INF=1000000000;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}

struct IOSetup{
    IOSetup(){
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout<<fixed<<setprecision(12);
    }
} iosetup;
 
template<typename T>
ostream &operator<<(ostream &os,const vector<T>&v){
    for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");
    return os;
}
template<typename T>
istream &operator>>(istream &is,vector<T>&v){
    for(T &x:v)is>>x;
    return is;
}
bool isPrime(ll num){
    if(num<2)           return false;
    else if(num==2)     return true;
    else if(num%2==0)   return false;
    for(ll i=3;i*i<=num;i+=2) if(num%i==0) return false;
    return true;
}



signed main(){
    // 2*3*(クソデカPrime)で遷移した後，5*(クソデカPrime)で遷移して，
    // (クソデカPrime)がダブる可能性があってかなりアなので
    // 最大素因数で集合ごと分けてシバくのがお気持ち


    ll n;cin>>n;
    vector<ll> p;
    rep(i,1300)if(isPrime(i)) p.push_back(i);
    
    vector<vector<ll>> g(p.size());
    for(ll i=2;i<=n;i++){
        int j;
        for(j=(int)p.size()-1;j>=0;j--){
            if(i%p[j]==0) break;
        }
        g[j].push_back(i);
    }

    const int m=12;
    ll dp[1<<m];
    rep(i,1<<m)dp[i]=0;

    for(auto &w:g){
        ll ndp[1<<m];
        rep(i,1<<m) ndp[i]=dp[i];
        for(auto &x:w){
            int mask=0;
            rep(i,m)if(x%p[i]==0) mask|=(1<<i);
            rep(bit,1<<m){
                if((bit&mask)==0) chmax(ndp[bit|mask],dp[bit]+x);
            }
        }
        rep(i,1<<m) dp[i]=ndp[i];
    }

    cout<<*max_element(dp,dp+(1<<m))<<endl;
    return 0;
}
0