結果

問題 No.458 異なる素数の和
ユーザー Yuu EguciYuu Eguci
提出日時 2020-12-05 14:07:16
言語 PyPy3
(7.3.15)
結果
MLE  
実行時間 -
コード長 4,778 bytes
コンパイル時間 376 ms
コンパイル使用メモリ 87,296 KB
実行使用メモリ 526,516 KB
最終ジャッジ日時 2023-10-13 21:33:45
合計ジャッジ時間 11,946 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge13
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 84 ms
75,804 KB
testcase_01 AC 490 ms
191,436 KB
testcase_02 AC 593 ms
219,052 KB
testcase_03 AC 170 ms
96,588 KB
testcase_04 AC 201 ms
108,712 KB
testcase_05 AC 1,200 ms
465,908 KB
testcase_06 AC 575 ms
218,256 KB
testcase_07 AC 89 ms
75,844 KB
testcase_08 AC 1,208 ms
468,528 KB
testcase_09 AC 118 ms
84,296 KB
testcase_10 AC 78 ms
71,316 KB
testcase_11 MLE -
testcase_12 AC 234 ms
71,552 KB
testcase_13 AC 74 ms
71,660 KB
testcase_14 AC 76 ms
71,688 KB
testcase_15 AC 76 ms
71,552 KB
testcase_16 AC 138 ms
84,236 KB
testcase_17 AC 80 ms
75,788 KB
testcase_18 AC 81 ms
75,808 KB
testcase_19 AC 74 ms
71,308 KB
testcase_20 AC 81 ms
75,752 KB
testcase_21 AC 73 ms
71,164 KB
testcase_22 AC 75 ms
71,220 KB
testcase_23 AC 80 ms
75,860 KB
testcase_24 AC 81 ms
75,816 KB
testcase_25 AC 75 ms
71,468 KB
testcase_26 AC 78 ms
75,564 KB
testcase_27 AC 547 ms
206,916 KB
testcase_28 MLE -
testcase_29 AC 99 ms
76,484 KB
testcase_30 AC 371 ms
157,772 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import sqrt
# from pprint import pprint
# import cProfile
# import pstats


def prime(N: int):

    # 素数一覧。
    primes = get_primes_until_n(N)
    # (value, weight) 化します。
    # それぞれの数字の価値は全部 1 です。なるべく多く入ってほしい = 全価値同じ、ってことだから。
    # HACK: 速度のため、わかりやすく (value, weight) にしていたのをやめます。
    # primes = [(1, _) for _ in primes_]
    # print(primes)

    # 表のいちばん左の列。行の数は、冒頭の例でいえば 18+1。重量 0 のときの行があるからひとつ多い。
    # None は、「その重さは作れない」ことを意味します。だって 1 とか素数足しても作れないでしょ?
    # dp = [[None] * (N + 1)]
    # だけど重さ 0 が、素数 0 個選択で作れることは分かってるので 0 に更新しときます。
    # dp[0][0] = 0
    # HACK: 速度のために、一番最初に全部作ることにしました。
    # HACK: 速度のために None -> -1 にしました。
    dp = [[0] + [-1] * N for i in range(len(primes) + 1)]
    # print(sys.getsizeof(dp))

    # 素数の個数分計算を行います。
    for i in range(len(primes)):

        # i の周回では primes[i] まで選べるときの計算をします。
        # 表でいえば i+1 列目ね。
        # value = primes[i][0]
        # weight = primes[i][1]
        # HACK: 速度のため、わかりやすく (value, weight) にしていたのをやめます。

        # i+1 列目を作っておきます。
        # dp.append([None] * (N + 1))
        # dp[i + 1][0] = 0
        # HACK: 速度のために、一番最初に全部作ることにしました。

        for w in range(N + 1):

            # HACK: 速度のために、下の if をひとつにします。
            if dp[i][w] != -1 and w + primes[i] <= N:
                dp[i + 1][w + primes[i]] = dp[i][w] + 1

            # dp[i][w] のマスに、 primes[i] を足したらどうなる? って検証をするんだけど、
            # そもそもこのマスが作成不可マス……どう素数を足しても作れない重さ……だったら計算する意味がないです。
            # if dp[i][w] == -1:
                # continue

            # このマス右隣(dp[i + 1][w])と比較します。
            # 現在のマス dp[i][w] には、 primes[i] まで選べるときの最大値が入っていて、
            # 右隣には、 primes[i+1] を足した場合の値が入ってます。
            # それが不適切(今のままのほうが大きい)場合があるから、更新します。
            # HACK: 速度のため max -> if にします。
            dp[i + 1][w] = max(dp[i][w], dp[i + 1][w])

            # じゃあこのマスに primes[i] を足してみます。
            # ひとつ足すわけだから、格納する列は i+1 になります。
            # 格納する行は、 w + primes[i] の重さ になります。
            # ただ N を超えるようならどうだっていいのでスキップです。
            # if w + primes[i] > N:
            #     continue
            # 格納する値は、いまのところからひとつ足すので +1 です。今回では value が全部 1。
            # dp[i + 1][w + primes[i]] = dp[i][w] + 1

    # dp テーブル作り終わりました。
    # この表には、素数が i 個選べるとき w をぴったり作れる素数の個数が入っています。
    # pprint(dp, width=200)
    # 今回知りたいのは、 N をぴったり作れるときの最大数です。
    # N をぴったり作れる個数の一覧がこれ。
    possible_values = [i[N] for i in dp]
    # その中で最大のものを返します。あっ、 None が混ざっているので None は -1 に変換します。
    # possible_values = [_ if _ else -1 for _ in possible_values]
    return max(possible_values)


def get_primes_until_n(n: int):
    S = []
    isprime = [True] * (n + 1)
    isprime[0] = False
    isprime[1] = False
    for i in range(2, int(sqrt(n)) + 1):
        if isprime[i]:
            j = i + i
            while j <= n:
                isprime[j] = False
                j = j + i

    for i in range(n + 1):
        if isprime[i]:
            S.append(i)
    return S


# print(prime(18) == 3)
# print(prime(4) == -1)
# print(prime(3) == 1)
# print(prime(1) == -1)
# print(prime(3344) == 41)
print(prime(int(input())))


# def test():
#     return prime(int(input()))


# profiler = cProfile.Profile()
# profiler.runcall(test)
# stats = pstats.Stats(profiler)
# stats.strip_dirs()
# stats.sort_stats('cumulative')
# stats.print_stats()
0