結果

問題 No.1760 Setwise Coprime
ユーザー PCTprobabilityPCTprobability
提出日時 2021-10-16 06:33:40
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 267 ms / 2,000 ms
コード長 1,302 bytes
コンパイル時間 365 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 95,552 KB
最終ジャッジ日時 2024-06-11 05:51:05
合計ジャッジ時間 6,460 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 35 ms
52,352 KB
testcase_01 AC 37 ms
52,224 KB
testcase_02 AC 90 ms
76,800 KB
testcase_03 AC 35 ms
52,352 KB
testcase_04 AC 34 ms
52,224 KB
testcase_05 AC 35 ms
52,224 KB
testcase_06 AC 94 ms
76,544 KB
testcase_07 AC 83 ms
74,112 KB
testcase_08 AC 86 ms
76,672 KB
testcase_09 AC 91 ms
76,800 KB
testcase_10 AC 54 ms
65,280 KB
testcase_11 AC 89 ms
76,672 KB
testcase_12 AC 86 ms
76,672 KB
testcase_13 AC 46 ms
60,544 KB
testcase_14 AC 85 ms
76,672 KB
testcase_15 AC 83 ms
76,672 KB
testcase_16 AC 88 ms
76,800 KB
testcase_17 AC 84 ms
76,672 KB
testcase_18 AC 53 ms
64,512 KB
testcase_19 AC 84 ms
76,672 KB
testcase_20 AC 83 ms
76,544 KB
testcase_21 AC 217 ms
89,380 KB
testcase_22 AC 204 ms
87,592 KB
testcase_23 AC 238 ms
93,424 KB
testcase_24 AC 159 ms
81,904 KB
testcase_25 AC 215 ms
89,372 KB
testcase_26 AC 179 ms
84,904 KB
testcase_27 AC 139 ms
79,960 KB
testcase_28 AC 122 ms
78,212 KB
testcase_29 AC 243 ms
93,896 KB
testcase_30 AC 236 ms
94,012 KB
testcase_31 AC 149 ms
80,716 KB
testcase_32 AC 225 ms
92,484 KB
testcase_33 AC 182 ms
84,104 KB
testcase_34 AC 242 ms
93,880 KB
testcase_35 AC 209 ms
90,648 KB
testcase_36 AC 263 ms
95,172 KB
testcase_37 AC 267 ms
95,392 KB
testcase_38 AC 266 ms
95,552 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod=998244353
def modpow(a,b):
  res=1
  while b:
    if b%2:
      res*=a
      res%=mod
    a*=a
    a%=mod
    b//=2
  return res
def lcmconv(x):
  n=len(x)
  n-=1
  for i in range(n,0,-1):
    for j in range(2*i,n+1,i):
      x[j]+=x[i]
  for i in range(1,n+1):
    x[i]*=x[i]
    x[i]%=mod
  for i in range(1,n+1):
    for j in range(2*i,n+1,i):
      x[j]-=x[i]
  for i in range(1,n+1):
    x[i]%=mod
  return x
n=int(input())
mebius=[1 for i in range(n+1)]
prime=[True for i in range(n+1)]
prime[1]=False
for i in range(1,n+1):
  if prime[i]==False:
    continue
  for j in range(2*i,n+1,i):
    prime[j]=False
for i in range(1,n+1):
  if prime[i]==False:
    continue
  for j in range(i*i,n+1,i*i):
    mebius[j]=0
  for j in range(i,n+1,i):
    mebius[j]*=-1
mebiussum=0
for i in range(1,n+1):
  mebius[i]%=mod
  mebiussum+=mebius[i]
mebiussum%=mod
ans=0
ans+=mebiussum*mebiussum
ans%=mod
twosum=0
for i in range(1,n+1):
  twosum+=modpow(2,n//i)*mebius[i]
  twosum%=mod
ans-=2*twosum*mebiussum
ans%=mod
x=[0 for i in range(n+1)]
for i in range(1,n+1):
  x[i]=modpow(2,n//i)*mebius[i]
  x[i]%=mod
x=lcmconv(x)
twosum*=twosum
twosum%=mod
for i in range(1,n+1):
  twosum-=x[i]
twosum%=mod
ans+=twosum
z=3*modpow(4,mod-2)%mod
for i in range(1,n+1):
  ans+=x[i]*modpow(z,n//i)
  ans%=mod
print(ans)
0