結果
| 問題 | No.103 素因数ゲーム リターンズ | 
| コンテスト | |
| ユーザー |  chocorusk | 
| 提出日時 | 2020-09-13 09:18:15 | 
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 62 ms / 5,000 ms | 
| コード長 | 900 bytes | 
| コンパイル時間 | 211 ms | 
| コンパイル使用メモリ | 12,672 KB | 
| 実行使用メモリ | 12,032 KB | 
| 最終ジャッジ日時 | 2024-06-11 16:27:45 | 
| 合計ジャッジ時間 | 2,896 ms | 
| ジャッジサーバーID (参考情報) | judge3 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 5 | 
| other | AC * 20 | 
ソースコード
import sys
read=sys.stdin.buffer.read
readline=sys.stdin.buffer.readline
readlines=sys.stdin.buffer.readlines
n=int(readline())
m=list(map(int, readline().split()))
def prime_sieve(n):
    isprime=[True]*(n+1)
    isprime[0]=False
    isprime[1]=False
    for i in range(2, n+1):
        if isprime[i]:
            for j in range(2*i, n+1, i):
                isprime[j]=False
    return isprime
mx=10000
isprime=prime_sieve(mx)
f=[[] for _ in range(mx+1)]
for p, b in enumerate(isprime):
    if not b:
        continue
    for i in range(p, mx+1, p):
        f[i].append(p)
dp=[0]*(mx+1)
for x in range(2, mx+1):
    st=set()
    for p in f[x]:
        x1=x//p
        st.add(dp[x1])
        if x1%p==0:
            st.add(dp[x1//p])
    for i in range(mx+1):
        if i not in st:
            dp[x]=i
            break
s=0
for x in m:
    s^=dp[x]
if s:
    print("Alice")
else:
    print("Bob")
            
            
            
        