結果

問題 No.1955 Not Prime
ユーザー titiatitia
提出日時 2022-05-22 18:13:22
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 2,361 bytes
コンパイル時間 152 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 210,196 KB
最終ジャッジ日時 2024-09-20 12:46:36
合計ジャッジ時間 6,476 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 83 ms
85,632 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 AC 78 ms
85,632 KB
testcase_06 AC 136 ms
100,480 KB
testcase_07 AC 324 ms
153,692 KB
testcase_08 AC 284 ms
135,696 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 AC 447 ms
210,192 KB
testcase_14 AC 193 ms
100,972 KB
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 AC 442 ms
208,500 KB
testcase_20 AC 83 ms
87,296 KB
testcase_21 AC 264 ms
111,248 KB
testcase_22 AC 81 ms
85,760 KB
testcase_23 WA -
testcase_24 AC 79 ms
85,504 KB
testcase_25 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
input = sys.stdin.readline
import math 
N=int(input())

AB=[tuple(map(int,input().split())) for i in range(N)]

x=10**6+10

L=math.floor(math.sqrt(x)) # 平方根を求める

Primelist=[i for i in range(x+1)]
Primelist[1]=0 # 1は素数でないので0にする.
 
for i in Primelist:
    if i>L:
        break
    if i==0:
        continue
    for j in range(2*i,x+1,i):
        Primelist[j]=0

Primes=[Primelist[j] for j in range(x+1) if Primelist[j]!=0]
Primes=set(Primes)

EDGE=[[] for i in range(2*N)]
EDGE_INV=[[] for i in range(2*N)]

# 0...N-1: x
# N...2N-1: ¬x

for i in range(N):
    a,b=AB[i]

    for j in range(N):
        c,d=AB[j]

        if not(int(str(a)+str(c)) in Primes):
            # ¬(x ⋀ y)
            EDGE[i].append(j+N)
            EDGE[j].append(i+N)

            EDGE_INV[j+N].append(i)
            EDGE_INV[i+N].append(j)

        if not(int(str(a)+str(d)) in Primes):
            # ¬(x ⋀ ¬y)
            EDGE[i].append(j)
            EDGE[j+N].append(i+N)

            EDGE_INV[j].append(i)
            EDGE_INV[i+N].append(j+N)


        if not(int(str(b)+str(c)) in Primes):
            # ¬(¬x ⋀ y)
            EDGE[i+N].append(j+N)
            EDGE[j].append(i)

            EDGE_INV[j+N].append(i+N)
            EDGE_INV[i].append(j)

        if not(int(str(b)+str(d)) in Primes):
            # ¬(¬x ⋀ ¬y)
            EDGE[i+N].append(j)
            EDGE[j+N].append(i)

            EDGE_INV[j].append(i+N)
            EDGE_INV[i].append(j+N)
            
            
QUE = list(range(2*N))
check=[0]*(2*N)
TOP_SORT=[]
 
def dfs(x):
    if check[x]==1:
        return
    check[x]=1
    
    for to in EDGE[x]:
        if check[to]==0:
            dfs(to)
 
    TOP_SORT.append(x) # 全ての点からDFSを行い, 帰りがけに点を答えに入れる
    check[x]=1
 
while QUE:
    x=QUE.pop()
    dfs(x)
 
USE=[0]*(2*N)
SCC=[]
 
def dfs2(x):
    Q=[x]
    USE[x]=1
    ANS=[]
 
    while Q:
        x=Q.pop()
        ANS.append(x)
        for to in EDGE_INV[x]:
            if USE[to]==0:
                USE[to]=1
                Q.append(to)
    return ANS
 
for x in TOP_SORT[::-1]:
    if USE[x]==0:
        SCC.append(dfs2(x))


for scc in SCC:
    SET=set(scc)

    for x in scc:
        if x+N in SET or x-N in SET:
            print("No")
            exit()

print("Yes")
0