結果

問題 No.719 Coprime
コンテスト
ユーザー vjudge1
提出日時 2025-12-26 23:35:53
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
結果
AC  
実行時間 3 ms / 3,000 ms
コード長 894 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 332 ms
コンパイル使用メモリ 43,716 KB
実行使用メモリ 7,852 KB
最終ジャッジ日時 2025-12-26 23:35:56
合計ジャッジ時間 2,503 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 61
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:8:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
    8 |         scanf("%d",&n);
      |         ~~~~~^~~~~~~~~

ソースコード

diff #
raw source code 

#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;
}
0