結果

問題 No.719 Coprime
ユーザー startcppstartcpp
提出日時 2016-07-02 22:27:46
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 128 ms / 3,000 ms
コード長 2,051 bytes
コンパイル時間 872 ms
コンパイル使用メモリ 65,200 KB
実行使用メモリ 6,068 KB
最終ジャッジ日時 2023-08-02 06:25:35
合計ジャッジ時間 4,318 ms
ジャッジサーバーID
(参考情報)
judge11 / judge12
このコードへのチャレンジ(β)

テストケース

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

ソースコード

diff #

#include <iostream>
#include <vector>
using namespace std;

typedef vector<int> Union;

bool isPrime[2526];	//値→素数か?
int toId[2526];		//素数→素数Id
int toNum[2526];	//素数Id→素数
vector<Union> uni;	//uni[i] = toNum[i]を最大の素因数とする値の集合
int dp[500][1<<15];
int n;

//usedの下位iビット目が1なら、toNum[i]を素因数とする値を使用済み (ただし、toNum[i] >= √Nとなるiについては未定義)
//primeId…最大の素因数がtoNum[primeId]の値から探索を開始する (primeId昇順に探索する…こうすると、常に素因数toNum[primeId]が未使用な状態からスタートできる)
//返り値…ここから選べる値の和の最大値
int dfs(int primeId, int used) {
	if (primeId == uni.size()) return 0;
	if (dp[primeId][used] != 0) return dp[primeId][used];
	
	int ret = dfs(primeId + 1, used);
	for (int i = 0; i < uni[primeId].size(); i++) {
		int value = uni[primeId][i];
		bool isOk = true;
		for (int j = 0; j < primeId && toNum[j] * toNum[j] < n; j++) {
			if ((used >> j) % 2 == 1 && value % toNum[j] == 0) {
				isOk = false;
				break;
			}
		}
		if (!isOk) continue;
		
		//値valueを使用する
		int newUsed = used;
		for (int j = 0; j <= primeId && toNum[j] * toNum[j] < n; j++) {
			if (value % toNum[j] == 0) {
				newUsed |= (1 << j);
			}
		}
		ret = max(ret, value + dfs(primeId + 1, newUsed));
	}
	return (dp[primeId][used] = ret);
}

void prepare() {
	int i, j;
	
	for (i = 0; i <= n; i++) isPrime[i] = true; isPrime[0] = false; isPrime[1] = false;
	for (i = 2; i <= n; i++) {
		if (isPrime[i]) {
			for (j = i * 2; j <= n; j += i) {
				isPrime[j] = false;
			}
		}
	}
	int id = 0;
	for (i = 2; i <= n; i++) {
		if (isPrime[i]) {
			toId[i] = id;
			toNum[id] = i;
			id++;
		}
	}
	
	uni.resize(id);
	for (i = 2; i <= n; i++) {
		for (j = i; j >= 2; j--) {
			if (i % j == 0 && isPrime[j]) {
				uni[toId[j]].push_back(i);
				break;
			}
		}
	}
}

int main() {
	cin >> n;
	prepare();
	cout << dfs(0, 0) << endl;
	return 0;
}
0