#62971 (Python3) No.308 素数は通れません

提出ソース

結果

問題	No.308 素数は通れません
ユーザー	rpy3cpp
提出日時	2015-12-05 10:48:58
言語	Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果	AC
実行時間	44 ms / 1,000 ms
コード長	998 bytes
コンパイル時間	476 ms
コンパイル使用メモリ	10,844 KB
実行使用メモリ	9,004 KB
最終ジャッジ日時	2023-08-18 13:22:05
合計ジャッジ時間	6,164 ms
ジャッジサーバーID （参考情報）	judge11 / judge13

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	19 ms 8,828 KB
testcase_01	AC	19 ms 8,876 KB
testcase_02	AC	19 ms 8,864 KB
testcase_03	AC	19 ms 8,804 KB
testcase_04	AC	20 ms 8,876 KB
testcase_05	AC	19 ms 8,948 KB
testcase_06	AC	19 ms 8,804 KB
testcase_07	AC	19 ms 8,828 KB
testcase_08	AC	19 ms 8,996 KB
testcase_09	AC	19 ms 8,848 KB
testcase_10	AC	19 ms 8,804 KB
testcase_11	AC	19 ms 8,864 KB
testcase_12	AC	19 ms 8,952 KB
testcase_13	AC	19 ms 8,772 KB
testcase_14	AC	19 ms 8,932 KB
testcase_15	AC	19 ms 8,904 KB
testcase_16	AC	19 ms 8,948 KB
testcase_17	AC	19 ms 8,932 KB
testcase_18	AC	19 ms 8,904 KB
testcase_19	AC	19 ms 8,808 KB
testcase_20	AC	19 ms 8,804 KB
testcase_21	AC	19 ms 8,776 KB
testcase_22	AC	19 ms 8,876 KB
testcase_23	AC	18 ms 8,848 KB
testcase_24	AC	19 ms 9,004 KB
testcase_25	AC	19 ms 8,928 KB
testcase_26	AC	19 ms 8,840 KB
testcase_27	AC	19 ms 8,812 KB
testcase_28	AC	20 ms 8,876 KB
testcase_29	AC	19 ms 8,864 KB
testcase_30	AC	19 ms 8,832 KB
testcase_31	AC	19 ms 8,828 KB
testcase_32	AC	19 ms 8,864 KB
testcase_33	AC	19 ms 8,812 KB
testcase_34	AC	19 ms 8,784 KB
testcase_35	AC	20 ms 8,952 KB
testcase_36	AC	19 ms 8,832 KB
testcase_37	AC	19 ms 8,932 KB
testcase_38	AC	19 ms 8,984 KB
testcase_39	AC	19 ms 8,932 KB
testcase_40	AC	18 ms 8,812 KB
testcase_41	AC	19 ms 8,804 KB
testcase_42	AC	19 ms 8,908 KB
testcase_43	AC	19 ms 8,828 KB
testcase_44	AC	19 ms 8,840 KB
testcase_45	AC	19 ms 8,776 KB
testcase_46	AC	19 ms 8,784 KB
testcase_47	AC	20 ms 8,832 KB
testcase_48	AC	19 ms 8,888 KB
testcase_49	AC	19 ms 8,828 KB
testcase_50	AC	19 ms 8,844 KB
testcase_51	AC	19 ms 8,732 KB
testcase_52	AC	19 ms 8,772 KB
testcase_53	AC	19 ms 8,932 KB
testcase_54	AC	19 ms 8,828 KB
testcase_55	AC	19 ms 8,812 KB
testcase_56	AC	19 ms 8,816 KB
testcase_57	AC	19 ms 8,932 KB
testcase_58	AC	19 ms 8,828 KB
testcase_59	AC	19 ms 8,836 KB
testcase_60	AC	19 ms 8,828 KB
testcase_61	AC	19 ms 8,760 KB
testcase_62	AC	19 ms 8,932 KB
testcase_63	AC	19 ms 8,936 KB
testcase_64	AC	19 ms 8,864 KB
testcase_65	AC	19 ms 8,996 KB
testcase_66	AC	20 ms 8,908 KB
testcase_67	AC	19 ms 8,804 KB
testcase_68	AC	19 ms 8,732 KB
testcase_69	AC	19 ms 8,928 KB
testcase_70	AC	19 ms 8,908 KB
testcase_71	AC	19 ms 8,780 KB
testcase_72	AC	19 ms 8,736 KB
testcase_73	AC	19 ms 8,780 KB
testcase_74	AC	19 ms 8,908 KB
testcase_75	AC	19 ms 8,744 KB
testcase_76	AC	22 ms 8,788 KB
testcase_77	AC	18 ms 8,836 KB
testcase_78	AC	19 ms 8,856 KB
testcase_79	AC	19 ms 8,812 KB
testcase_80	AC	23 ms 8,828 KB
testcase_81	AC	28 ms 8,944 KB
testcase_82	AC	19 ms 8,844 KB
testcase_83	AC	19 ms 8,832 KB
testcase_84	AC	31 ms 8,952 KB
testcase_85	AC	19 ms 8,808 KB
testcase_86	AC	33 ms 8,932 KB
testcase_87	AC	36 ms 8,860 KB
testcase_88	AC	19 ms 8,944 KB
testcase_89	AC	19 ms 8,936 KB
testcase_90	AC	42 ms 8,984 KB
testcase_91	AC	19 ms 8,784 KB
testcase_92	AC	19 ms 8,888 KB
testcase_93	AC	19 ms 8,948 KB
testcase_94	AC	20 ms 8,936 KB
testcase_95	AC	20 ms 8,872 KB
testcase_96	AC	19 ms 8,888 KB
testcase_97	AC	19 ms 8,736 KB
testcase_98	AC	19 ms 8,988 KB
testcase_99	AC	19 ms 8,836 KB
testcase_100	AC	20 ms 8,896 KB
testcase_101	AC	44 ms 8,828 KB
testcase_102	AC	19 ms 8,936 KB
testcase_103	AC	20 ms 8,984 KB
testcase_104	AC	19 ms 8,804 KB
testcase_105	AC	19 ms 8,912 KB
testcase_106	AC	19 ms 8,860 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

import random

def MillerRabin(q, k = 1000):
    '''
    素数判定の ミラー・ラビン テスト 以下のコードをもとに作成。
    http://d.hatena.ne.jp/pashango_p/20090704/1246692091
    '''
    if q == 2:
        return True
    if q == 1 or not (q & 1):
        return False
    d = (q - 1) >> 1
    while not (d & 1):
        d >>= 1
    for i in range(k):
        a = random.randint(1, q - 1)
        t = d
        y = pow(a, t, q)
        while t != q - 1 and y != 1 and y != q - 1: 
            y = pow(y, 2, q)
            t <<= 1
        if y != q - 1 and not (t & 1):
            return False
    return True

def solve(N):
    dic = {4:3, 6:5, 8:7, 9:7, 10:7, 12:11, 14:13, 15:7, 16:7, 18:8, 20:19, 21:19, 22:7, 24:23, 25:23, 26:8}
    if N in dic:
        return dic[N]
    if N == 49:
        return 14
    if N < 57:
        return 8
    if N % 8 == 1:
        M = N - 8
        if MillerRabin(M):
            return 14
    return 8

N = int(input())
print(solve(N))

yukicoder

結果

テストケース

ソースコード