結果
| 問題 | No.3691 Calculate Mu Sum |
| コンテスト | |
| ユーザー |
titia
|
| 提出日時 | 2026-09-09 03:39:54 |
| 言語 | PyPy3 (7.3.23 + ACL) |
| 結果 |
AC
不安定
|
| 実行時間 | 1,805 ms / 2,000 ms |
| + 318µs | |
| コード長 | 1,232 bytes |
| 記録 | |
| コンパイル時間 | 78 ms |
| コンパイル使用メモリ | 81,408 KB |
| 実行使用メモリ | 240,020 KB |
| 最終ジャッジ日時 | 2026-09-09 03:42:32 |
| 合計ジャッジ時間 | 12,935 ms |
|
ジャッジサーバーID (参考情報) |
judge2_0 / judge3_0 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 11 |
ソースコード
import sys
input = sys.stdin.readline
# エラトステネスの篩を用いた素因数分解・約数列挙
MAX=10**7+10 # 使いたい最大値を指定
# Sieve[i]で、iの最も小さい1でない約数を返す。
Sieve=[i for i in range(MAX)]
for i in range(2,MAX):
if Sieve[i]!=i:
continue
for j in range(i,MAX,i):
if Sieve[j]==j:
Sieve[j]=i
# 素因数分解
def fact(x):
D=dict()
while x!=1:
k=Sieve[x]
if k in D:
D[k]+=1
else:
D[k]=1
x//=k
return D
# 約数列挙
def faclist(x):
LIST=[1]
while x!=1:
k=Sieve[x]
count=0
while x%k==0:
count+=1
x//=k
LIST2=[]
for l in LIST:
for i in range(1,count+1):
LIST2.append(l*k**i)
LIST+=LIST2
return LIST
#print("!")
N=int(input())
ANS=[-10]*(N+1)
ANS[0]=0
ANS[1]=1
for i in range(2,N):
if i*i<=N:
ANS[i*i]=0
if ANS[i]==0:
for j in range(i,N+1,i):
ANS[j]=0
for i in range(1,N+1):
if ANS[i]==-10:
if len(fact(i))%2==0:
ANS[i]=1
else:
ANS[i]=-1
print(sum(ANS))
titia