結果
問題 | No.719 Coprime |
ユーザー |
![]() |
提出日時 | 2020-10-15 01:02:09 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 11 ms / 3,000 ms |
コード長 | 2,170 bytes |
コンパイル時間 | 1,782 ms |
コンパイル使用メモリ | 201,204 KB |
最終ジャッジ日時 | 2025-01-15 07:21:33 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 61 |
ソースコード
#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; }