結果

問題 No.732 3PrimeCounting
ユーザー vwxyzvwxyz
提出日時 2022-10-03 02:39:15
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 3,290 bytes
コンパイル時間 268 ms
コンパイル使用メモリ 82,104 KB
実行使用メモリ 96,468 KB
最終ジャッジ日時 2024-06-07 19:18:22
合計ジャッジ時間 16,081 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 150 ms
96,468 KB
testcase_01 AC 151 ms
89,464 KB
testcase_02 AC 154 ms
89,168 KB
testcase_03 AC 150 ms
89,364 KB
testcase_04 AC 148 ms
89,152 KB
testcase_05 AC 149 ms
89,576 KB
testcase_06 AC 150 ms
89,120 KB
testcase_07 AC 154 ms
89,108 KB
testcase_08 AC 150 ms
89,120 KB
testcase_09 AC 153 ms
89,260 KB
testcase_10 AC 152 ms
89,004 KB
testcase_11 AC 152 ms
89,528 KB
testcase_12 AC 153 ms
89,388 KB
testcase_13 AC 152 ms
89,088 KB
testcase_14 AC 153 ms
89,204 KB
testcase_15 AC 152 ms
89,364 KB
testcase_16 AC 147 ms
89,544 KB
testcase_17 AC 150 ms
89,320 KB
testcase_18 AC 153 ms
89,252 KB
testcase_19 AC 148 ms
89,196 KB
testcase_20 AC 159 ms
90,124 KB
testcase_21 AC 183 ms
90,284 KB
testcase_22 AC 180 ms
90,332 KB
testcase_23 AC 152 ms
89,952 KB
testcase_24 AC 154 ms
89,756 KB
testcase_25 AC 217 ms
90,648 KB
testcase_26 AC 196 ms
91,000 KB
testcase_27 AC 166 ms
89,904 KB
testcase_28 AC 169 ms
90,048 KB
testcase_29 AC 178 ms
90,464 KB
testcase_30 AC 174 ms
90,388 KB
testcase_31 AC 185 ms
90,040 KB
testcase_32 AC 201 ms
90,792 KB
testcase_33 AC 205 ms
90,940 KB
testcase_34 AC 205 ms
91,068 KB
testcase_35 AC 183 ms
90,324 KB
testcase_36 AC 186 ms
90,008 KB
testcase_37 AC 159 ms
90,136 KB
testcase_38 AC 158 ms
90,044 KB
testcase_39 AC 191 ms
90,320 KB
testcase_40 AC 187 ms
90,176 KB
testcase_41 AC 180 ms
90,280 KB
testcase_42 AC 184 ms
90,332 KB
testcase_43 AC 181 ms
90,296 KB
testcase_44 AC 180 ms
90,308 KB
testcase_45 AC 172 ms
89,968 KB
testcase_46 AC 172 ms
89,900 KB
testcase_47 AC 175 ms
90,240 KB
testcase_48 AC 228 ms
91,092 KB
testcase_49 AC 226 ms
90,868 KB
testcase_50 AC 182 ms
90,232 KB
testcase_51 AC 180 ms
89,924 KB
testcase_52 AC 166 ms
90,128 KB
testcase_53 AC 364 ms
93,804 KB
testcase_54 TLE -
testcase_55 -- -
testcase_56 -- -
testcase_57 -- -
testcase_58 -- -
testcase_59 -- -
testcase_60 -- -
testcase_61 -- -
testcase_62 -- -
testcase_63 -- -
testcase_64 -- -
testcase_65 -- -
testcase_66 -- -
testcase_67 -- -
testcase_68 -- -
testcase_69 -- -
testcase_70 -- -
testcase_71 -- -
testcase_72 -- -
testcase_73 -- -
testcase_74 -- -
testcase_75 -- -
testcase_76 -- -
testcase_77 -- -
testcase_78 -- -
testcase_79 -- -
testcase_80 -- -
testcase_81 -- -
testcase_82 -- -
testcase_83 -- -
testcase_84 -- -
testcase_85 -- -
testcase_86 -- -
testcase_87 -- -
testcase_88 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

from ast import Mod
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

class Prime:
    def __init__(self,N):
        assert N<=10**8
        self.smallest_prime_factor=[None]*(N+1)
        for i in range(2,N+1,2):
            self.smallest_prime_factor[i]=2
        n=int(N**.5)+1
        for p in range(3,n,2):
            if self.smallest_prime_factor[p]==None:
                self.smallest_prime_factor[p]=p
                for i in range(p**2,N+1,2*p):
                    if self.smallest_prime_factor[i]==None:
                        self.smallest_prime_factor[i]=p
        for p in range(n,N+1):
            if self.smallest_prime_factor[p]==None:
                self.smallest_prime_factor[p]=p
        self.primes=[p for p in range(N+1) if p==self.smallest_prime_factor[p]]

    def Factorize(self,N):
        assert N>=1
        factors=defaultdict(int)
        if N<=len(self.smallest_prime_factor)-1:
            while N!=1:
                factors[self.smallest_prime_factor[N]]+=1
                N//=self.smallest_prime_factor[N]
        else:
            for p in self.primes:
                while N%p==0:
                    N//=p
                    factors[p]+=1
                if N<p*p:
                    if N!=1:
                        factors[N]+=1
                    break
                if N<=len(self.smallest_prime_factor)-1:
                    while N!=1:
                        factors[self.smallest_prime_factor[N]]+=1
                        N//=self.smallest_prime_factor[N]
                    break
            else:
                if N!=1:
                    factors[N]+=1
        return factors

    def Divisors(self,N):
        assert N>0
        divisors=[1]
        for p,e in self.Factorize(N).items():
            pow_p=[1]
            for _ in range(e):
                pow_p.append(pow_p[-1]*p)
            divisors=[i*j for i in divisors for j in pow_p]
        return divisors

    def Is_Prime(self,N):
        return N==self.smallest_prime_factor[N]

    def Totient(self,N):
        for p in self.Factorize(N).keys():
            N*=p-1
            N//=p
        return N

    def Mebius(self,N):
        fact=self.Factorize(N)
        for e in fact.values():
            if e>=2:
                return 0
        else:
            if len(fact)%2==0:
                return 1
            else:
                return -1

N=int(readline())
P=Prime(3*N)
primes=[p for p in P.primes if p<=N]
cnt=defaultdict(int)
le=len(primes)
ans=0
for b in range(le-1,-1,-1):
    for a in range(b-1,-1,-1):
        ans+=cnt[primes[a]+primes[b]]
    for p in P.primes:
        cnt[p-primes[b]]+=1
print(ans)
0