#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
実行時間	64 ms / 1,000 ms
コード長	998 bytes
コンパイル時間	222 ms
コンパイル使用メモリ	12,800 KB
実行使用メモリ	11,264 KB
最終ジャッジ日時	2024-11-27 18:40:18
合計ジャッジ時間	6,185 ms
ジャッジサーバーID （参考情報）	judge1 / judge3

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	33 ms 11,136 KB
testcase_01	AC	33 ms 11,136 KB
testcase_02	AC	33 ms 11,264 KB
testcase_03	AC	33 ms 11,008 KB
testcase_04	AC	33 ms 11,264 KB
testcase_05	AC	33 ms 11,008 KB
testcase_06	AC	32 ms 11,008 KB
testcase_07	AC	32 ms 11,008 KB
testcase_08	AC	33 ms 11,136 KB
testcase_09	AC	36 ms 11,264 KB
testcase_10	AC	32 ms 11,136 KB
testcase_11	AC	33 ms 11,008 KB
testcase_12	AC	32 ms 11,136 KB
testcase_13	AC	32 ms 11,136 KB
testcase_14	AC	32 ms 11,264 KB
testcase_15	AC	32 ms 11,136 KB
testcase_16	AC	32 ms 11,136 KB
testcase_17	AC	33 ms 11,136 KB
testcase_18	AC	33 ms 11,008 KB
testcase_19	AC	32 ms 11,008 KB
testcase_20	AC	32 ms 11,008 KB
testcase_21	AC	33 ms 11,008 KB
testcase_22	AC	32 ms 11,264 KB
testcase_23	AC	32 ms 11,136 KB
testcase_24	AC	32 ms 11,008 KB
testcase_25	AC	32 ms 11,136 KB
testcase_26	AC	33 ms 11,136 KB
testcase_27	AC	33 ms 11,008 KB
testcase_28	AC	32 ms 11,136 KB
testcase_29	AC	33 ms 11,008 KB
testcase_30	AC	32 ms 11,008 KB
testcase_31	AC	33 ms 11,136 KB
testcase_32	AC	33 ms 11,136 KB
testcase_33	AC	32 ms 11,136 KB
testcase_34	AC	32 ms 11,008 KB
testcase_35	AC	32 ms 11,136 KB
testcase_36	AC	32 ms 11,136 KB
testcase_37	AC	33 ms 11,136 KB
testcase_38	AC	33 ms 11,008 KB
testcase_39	AC	32 ms 11,136 KB
testcase_40	AC	32 ms 11,008 KB
testcase_41	AC	33 ms 11,136 KB
testcase_42	AC	33 ms 11,136 KB
testcase_43	AC	33 ms 11,008 KB
testcase_44	AC	32 ms 11,136 KB
testcase_45	AC	33 ms 11,136 KB
testcase_46	AC	32 ms 11,136 KB
testcase_47	AC	32 ms 11,008 KB
testcase_48	AC	32 ms 11,136 KB
testcase_49	AC	32 ms 11,264 KB
testcase_50	AC	32 ms 11,008 KB
testcase_51	AC	32 ms 11,008 KB
testcase_52	AC	33 ms 11,136 KB
testcase_53	AC	32 ms 11,008 KB
testcase_54	AC	33 ms 11,008 KB
testcase_55	AC	32 ms 11,008 KB
testcase_56	AC	32 ms 11,136 KB
testcase_57	AC	34 ms 11,136 KB
testcase_58	AC	32 ms 11,136 KB
testcase_59	AC	32 ms 11,136 KB
testcase_60	AC	32 ms 11,008 KB
testcase_61	AC	33 ms 11,136 KB
testcase_62	AC	31 ms 11,136 KB
testcase_63	AC	33 ms 11,008 KB
testcase_64	AC	32 ms 11,008 KB
testcase_65	AC	32 ms 11,136 KB
testcase_66	AC	33 ms 11,136 KB
testcase_67	AC	32 ms 11,008 KB
testcase_68	AC	32 ms 11,136 KB
testcase_69	AC	32 ms 11,264 KB
testcase_70	AC	32 ms 11,008 KB
testcase_71	AC	32 ms 11,136 KB
testcase_72	AC	33 ms 11,136 KB
testcase_73	AC	32 ms 11,008 KB
testcase_74	AC	31 ms 11,136 KB
testcase_75	AC	32 ms 11,136 KB
testcase_76	AC	36 ms 11,136 KB
testcase_77	AC	31 ms 11,136 KB
testcase_78	AC	32 ms 11,136 KB
testcase_79	AC	32 ms 11,008 KB
testcase_80	AC	37 ms 11,264 KB
testcase_81	AC	44 ms 11,008 KB
testcase_82	AC	33 ms 11,008 KB
testcase_83	AC	32 ms 11,008 KB
testcase_84	AC	49 ms 11,136 KB
testcase_85	AC	32 ms 11,136 KB
testcase_86	AC	51 ms 11,008 KB
testcase_87	AC	54 ms 11,136 KB
testcase_88	AC	33 ms 11,136 KB
testcase_89	AC	33 ms 11,008 KB
testcase_90	AC	61 ms 11,008 KB
testcase_91	AC	33 ms 11,008 KB
testcase_92	AC	33 ms 11,136 KB
testcase_93	AC	33 ms 11,008 KB
testcase_94	AC	32 ms 11,136 KB
testcase_95	AC	33 ms 11,136 KB
testcase_96	AC	32 ms 11,136 KB
testcase_97	AC	32 ms 11,264 KB
testcase_98	AC	35 ms 11,008 KB
testcase_99	AC	32 ms 11,008 KB
testcase_100	AC	33 ms 11,008 KB
testcase_101	AC	64 ms 11,008 KB
testcase_102	AC	33 ms 11,136 KB
testcase_103	AC	32 ms 11,008 KB
testcase_104	AC	32 ms 11,136 KB
testcase_105	AC	33 ms 11,136 KB
testcase_106	AC	33 ms 11,136 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

結果

テストケース

ソースコード