結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 150 ms
228,732 KB
testcase_01 AC 149 ms
89,088 KB
testcase_02 AC 149 ms
89,404 KB
testcase_03 AC 150 ms
88,960 KB
testcase_04 AC 148 ms
89,344 KB
testcase_05 AC 151 ms
89,216 KB
testcase_06 AC 148 ms
89,344 KB
testcase_07 AC 150 ms
89,472 KB
testcase_08 AC 150 ms
89,472 KB
testcase_09 AC 152 ms
89,344 KB
testcase_10 AC 151 ms
89,088 KB
testcase_11 AC 149 ms
89,216 KB
testcase_12 AC 148 ms
89,344 KB
testcase_13 AC 152 ms
89,216 KB
testcase_14 AC 154 ms
89,088 KB
testcase_15 AC 153 ms
89,088 KB
testcase_16 AC 150 ms
89,472 KB
testcase_17 AC 150 ms
89,216 KB
testcase_18 AC 151 ms
89,344 KB
testcase_19 AC 154 ms
89,088 KB
testcase_20 AC 159 ms
89,900 KB
testcase_21 AC 182 ms
90,064 KB
testcase_22 AC 181 ms
89,984 KB
testcase_23 AC 154 ms
89,728 KB
testcase_24 AC 154 ms
89,856 KB
testcase_25 AC 215 ms
91,076 KB
testcase_26 AC 197 ms
91,008 KB
testcase_27 AC 166 ms
89,728 KB
testcase_28 AC 160 ms
89,856 KB
testcase_29 AC 175 ms
90,300 KB
testcase_30 AC 170 ms
90,240 KB
testcase_31 AC 183 ms
90,228 KB
testcase_32 AC 201 ms
90,752 KB
testcase_33 AC 200 ms
91,136 KB
testcase_34 AC 200 ms
90,880 KB
testcase_35 AC 181 ms
90,112 KB
testcase_36 AC 180 ms
90,112 KB
testcase_37 AC 159 ms
90,112 KB
testcase_38 AC 156 ms
89,728 KB
testcase_39 AC 186 ms
90,112 KB
testcase_40 AC 183 ms
90,496 KB
testcase_41 AC 176 ms
90,368 KB
testcase_42 AC 181 ms
90,112 KB
testcase_43 AC 177 ms
90,112 KB
testcase_44 AC 177 ms
90,112 KB
testcase_45 AC 167 ms
89,984 KB
testcase_46 AC 168 ms
89,856 KB
testcase_47 AC 171 ms
90,240 KB
testcase_48 AC 223 ms
90,880 KB
testcase_49 AC 224 ms
90,880 KB
testcase_50 AC 179 ms
90,240 KB
testcase_51 AC 178 ms
90,112 KB
testcase_52 AC 162 ms
89,896 KB
testcase_53 AC 360 ms
93,568 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