結果
問題 | No.1747 Many Formulae 2 |
ユーザー |
![]() |
提出日時 | 2023-05-25 06:23:10 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 50 ms / 2,000 ms |
コード長 | 1,521 bytes |
コンパイル時間 | 106 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 11,904 KB |
最終ジャッジ日時 | 2024-12-24 03:29:08 |
合計ジャッジ時間 | 2,200 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 19 |
ソースコード
import bisect import copy import decimal import fractions import heapq import itertools import math import random import sys import time from collections import Counter,deque,defaultdict from functools import lru_cache,reduce from heapq import heappush,heappop,heapify,heappushpop,_heappop_max,_heapify_max def _heappush_max(heap,item): heap.append(item) heapq._siftdown_max(heap, 0, len(heap)-1) def _heappushpop_max(heap, item): if heap and item < heap[0]: item, heap[0] = heap[0], item heapq._siftup_max(heap, 0) return item from math import gcd as GCD read=sys.stdin.read readline=sys.stdin.readline readlines=sys.stdin.readlines write=sys.stdout.write def Miller_Rabin_Primality_Test(N): if N==1: return False NN=N-1 NN=NN//(NN&-NN) if N<4759123141: lst=[2,7,61] elif N<341550071728321: lst=[2,3,5,7,11,13,17] else: lst=[2,3,5,7,11,13,17,19,23,29,31,37] if N in lst: return True for a in lst: n=NN p=pow(a,n,N) if p==1: continue while p!=N-1: p=p*p%N if p==1 or n==N-1: return False n<<=1 return True S=readline().rstrip() le=len(S) ans=0 for bit in range(1<<le-1): bound=[0] for i in range(1,le): if bit&1<<i-1: bound.append(i) bound.append(le) s=0 for l,r in zip(bound,bound[1:]): s+=int(S[l:r]) if Miller_Rabin_Primality_Test(s): ans+=1 print(ans)