結果
問題 | No.385 カップ麺生活 |
ユーザー |
|
提出日時 | 2017-07-25 12:37:25 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,143 bytes |
コンパイル時間 | 877 ms |
コンパイル使用メモリ | 115,796 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-12 21:09:36 |
合計ジャッジ時間 | 1,943 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 32 |
ソースコード
import std.algorithm, std.conv, std.range, std.stdio, std.string;void main(){auto m = readln.chomp.to!int;auto n = readln.chomp.to!size_t;auto c = readln.split.to!(int[]);auto dp = new int[](m+1);dp[1..$] = int.min;foreach (i; 0..n)foreach (j; c[i]..m+1)dp[j] = max(dp[j], dp[j-c[i]] + 1);auto p = new bool[](m+1);foreach (pi; primes(m)) p[pi] = true;auto r = 0L, ma = 0L;foreach (i; 0..m+1)if (p[i] && dp[m-i] > 0)r += dp[m-i];writeln(r + dp.maxElement);}pure T[] primes(T)(T n){import std.algorithm, std.bitmanip, std.conv, std.range;auto sieve = BitArray();sieve.length((n + 1) / 2);sieve = ~sieve;foreach (p; 1..((nsqrt(n) - 1) / 2 + 1))if (sieve[p])for (auto q = p * 3 + 1; q < (n + 1) / 2; q += p * 2 + 1)sieve[q] = false;auto r = sieve.bitsSet.map!(to!T).map!("a * 2 + 1").array;r[0] = 2;return r;}pure T nsqrt(T)(T n){import std.algorithm, std.conv, std.range, core.bitop;if (n <= 1) return n;T m = 1 << (n.bsr / 2 + 1);return iota(1, m).map!"a * a".assumeSorted!"a <= b".lowerBound(n).length.to!T;}