結果
| 問題 | No.719 Coprime |
| コンテスト | |
| ユーザー |
vjudge1
|
| 提出日時 | 2025-12-26 23:35:53 |
| 言語 | C++17 (gcc 15.2.0 + boost 1.90.0) |
| 結果 |
AC
|
| 実行時間 | 3 ms / 3,000 ms |
| + 286µs | |
| コード長 | 894 bytes |
| 記録 | |
| コンパイル時間 | 181 ms |
| コンパイル使用メモリ | 54,236 KB |
| 実行使用メモリ | 5,888 KB |
| 最終ジャッジ日時 | 2026-07-19 20:43:31 |
| 合計ジャッジ時間 | 2,732 ms |
|
ジャッジサーバーID (参考情報) |
judge3_0 / judge2_0 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 61 |
ソースコード
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=2055;
int dp[N],g[N],pr[300],n,cnt,isin[N];
bool vis[N];
int main() {
scanf("%d",&n);
for(int i=2;i<=n;i++) {
if(!vis[i]) pr[++cnt]=i;
for(int j=1;pr[j]*i<=n;j++) {
vis[pr[j]*i]=1;
if(i%pr[j]==0) break;
}
}
for(int i=2;i<=n;i++) isin[i]=-1;
for(int i=1;i<=min(cnt,11);i++) {
for(int j=1;j*pr[i]<=n;j++) {
if(isin[j]!=-1) isin[j*pr[i]]=(isin[j]|(1<<i-1));
}
}
for(int i=2;i<=n;i++) {
if(isin[i]!=-1) {
for(int j=0;j<2048;j++) {
if((j&isin[i])==0) dp[j|isin[i]]=max(dp[j|isin[i]],dp[j]+i);
}
}
}
for(int i=12;i<=cnt;i++) {
for(int j=1;pr[i]*j<=n;j++) {
if(isin[j]!=-1) {
for(int k=0;k<2048;k++) {
if((k&isin[j])==0) g[k|isin[j]]=max(g[k|isin[j]],dp[k]+pr[i]*j);
}
}
}
for(int j=0;j<2048;j++) dp[j]=g[j];
}
printf("%d\n",dp[(1<<min(cnt,11))-1]);
return 0;
}
vjudge1