結果

問題 No.843 Triple Primes
ユーザー Navier_Boltzmann
提出日時 2024-01-10 22:29:45
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 163 ms / 2,000 ms
コード長 1,336 bytes
コンパイル時間 152 ms
コンパイル使用メモリ 82,300 KB
実行使用メモリ 111,080 KB
最終ジャッジ日時 2024-09-27 20:22:05
合計ジャッジ時間 8,275 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import *
from itertools import *
from functools import *
from heapq import *
import sys,math
input = sys.stdin.readline


class prime_factorize():
    
    def __init__(self,M=10**6):
        self.sieve = [-1]*(M+1)
        self.sieve[1] = 1
        self.p = [False]*(M+1)
        self.mu = [1]*(M+1)
        
        for i in range(2,M+1):
            if self.sieve[i] == -1:
                self.p[i] = True
                
                i2 = i**2
                for j in range(i2,M+1,i2):
                    
                    self.mu[j] = 0
                
                
                for j in range(i,M+1,i):
                    self.sieve[j] = i
                    
                    self.mu[j] *= -1
                    
    def factors(self,x):
        tmp = []
        while self.sieve[x] != x:
            tmp.append(self.sieve[x])
            x //= self.sieve[x]
        tmp.append(self.sieve[x])
        return tmp
        
    def is_prime(self,x):
        return self.p[x]
        
    def mobius(self,x):
        return self.mu[x]

N = int(input())
pf = prime_factorize()
ans = 0
S = set()
P = set()
for i in range(2,N+1):
    if pf.is_prime(i):
        P.add(i)
        S.add(i**2)

if N==1:
    print(0)
    exit()
p = 2
for q in P:
    if (q+2) in S:
        ans += 1
print(ans*2 - 1)
0