結果
問題 | No.719 Coprime |
ユーザー |
![]() |
提出日時 | 2019-01-22 11:44:36 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 107 ms / 3,000 ms |
コード長 | 1,585 bytes |
コンパイル時間 | 839 ms |
コンパイル使用メモリ | 100,152 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-15 13:31:31 |
合計ジャッジ時間 | 3,262 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 61 |
ソースコード
#include <cstdio> #include <cstring> #include <iostream> #include <string> #include <cmath> #include <bitset> #include <vector> #include <map> #include <set> #include <queue> #include <deque> #include <algorithm> #include <complex> #include <unordered_map> #include <unordered_set> #include <random> #include <cassert> using namespace std; typedef long long int ll; typedef pair<ll, int> P; const int MAX=1262; bool isprime[MAX]; vector<int> prime; void sieve(){ for(int i=3; i<MAX; i+=2){ isprime[i]=1; } isprime[2]=1; prime.push_back(2); for(int i=3; i<MAX; i++){ if(isprime[i]){ prime.push_back(i); for(int j=2*i; j<MAX; j+=i){ isprime[j]=0; } } } return; } int main() { sieve(); int n; cin>>n; int b=(int)sqrt((double)n)+1; int t=lower_bound(prime.begin(), prime.end(), b)-prime.begin(); int c=0; bool ok[1263]={}; int dp[2][1<<11]={}; for(int i=t; i<prime.size(); i++){ if(prime[i]>n) break; for(int l=0; l<(1<<t); l++) dp[c^1][l]=dp[c][l]; for(int l=0; l<(1<<t); l++){ for(int j=1; j*prime[i]<=n; j++){ ok[j*prime[i]]=1; int p=0; for(int k=0; k<t; k++){ if(j%prime[k]==0) p^=(1<<k); } if(l&p) continue; dp[c^1][l^p]=max(dp[c][l]+j*prime[i], dp[c^1][l^p]); } } c^=1; } for(int i=2; i<=n; i++){ if(ok[i]) continue; for(int l=0; l<(1<<t); l++) dp[c^1][l]=dp[c][l]; for(int l=0; l<(1<<t); l++){ int p=0; for(int k=0; k<t; k++){ if(i%prime[k]==0) p^=(1<<k); } if(l&p) continue; dp[c^1][l^p]=max(dp[c][l]+i, dp[c^1][l^p]); } c^=1; } cout<<dp[c][(1<<t)-1]<<endl; return 0; }