QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#409342#5040. Dirichlet $k$-th rootkangkunmaTL 122ms24532kbPython3889b2024-05-11 22:09:402024-05-11 22:09:40

Judging History

你现在查看的是最新测评结果

  • [2024-05-11 22:09:40]
  • 评测
  • 测评结果:TL
  • 用时:122ms
  • 内存:24532kb
  • [2024-05-11 22:09:40]
  • 提交

answer

def b(x):
    t=[]
    while x:t+=[x%2];x//=2
    return t
D=[[]for i in range(100001)]
for i in range(1,317):
    for j in range(i*i,100001,i):D[j]+=[i]
k,n=map(int,input().split());l=[*map(int,input().split())];M=998244353;B=b(n);L=len(B);div=pow(n,-1,M)
X=[[0,1]+[0]*(k-1) for i in range(len(B))]
Y=[[0,1]+[0]*(k-1) for i in range(sum(B)+1)]
for i in range(2,k+1): #구할 func.값
    c=0
    for j in range(L-1): 
        for p in D[i]:X[j+1][i]=(X[j+1][i]+2*X[j][p]*X[j][i//p]-pow(X[j][p],2,M)*(p**2==i))
    for j in range(L):
        if B[j]:
            for p in D[i]:Y[c+1][i]=(Y[c+1][i]+Y[c][p]*X[j][i//p]+Y[c][i//p]*X[j][p]-Y[c][p]*X[j][p]*(p**2==i))
            c+=1
    a=((l[i-1]-Y[-1][i])*div)%M
    t,d=0,0
    for j in range(L):
        w=pow(2,j,M);X[j][i]=(X[j][i]+(a*w))%M
        if B[j]:t=(t+w)%M;d+=1;Y[d][i]=(Y[d][i]+a*t)%M
print(*X[0][1::])

详细

Test #1:

score: 100
Accepted
time: 106ms
memory: 24532kb

input:

5 2
1 8 4 26 6

output:

1 4 2 5 3

result:

ok 5 number(s): "1 4 2 5 3"

Test #2:

score: 0
Accepted
time: 107ms
memory: 24336kb

input:

1000 10
1 516238479 765833625 408427111 593035929 984641872 149832309 743527500 533342020 315774476 81082172 461044503 235289865 539969249 564590101 736090497 933661679 394417077 377363949 717587561 143428748 321957671 413061277 166846744 532438287 802900949 101236134 182814904 328306207 83339054 87...

output:

1 750394895 575705539 383612918 758074640 870483127 713754278 172378328 328377975 148157627 607054829 769165678 522651163 979451027 321052733 428161561 792137215 595294435 736507442 784558508 35541967 129771597 141130563 667029666 976562087 614420517 890540347 922620141 132655056 898078185 586990958...

result:

ok 1000 numbers

Test #3:

score: 0
Accepted
time: 116ms
memory: 24296kb

input:

1000 10
1 836321215 611344545 176796302 574251412 137508398 761666720 357050290 801608050 863745155 811911311 829248030 265662561 683320811 920044294 338428495 840771823 558535163 926728514 375903056 522656042 909868778 153190814 812151247 960179729 591861004 169384211 295767075 841728125 81742549 5...

output:

1 582754298 560256631 598696700 656371753 428716095 76166672 762502364 234864311 233488066 380664437 883840099 326039562 115077358 562124625 491468731 982497100 88908472 292321722 938842515 680178128 840510694 214967952 569450139 579949075 97329455 173736151 491492784 583294989 144817858 358079034 9...

result:

ok 1000 numbers

Test #4:

score: 0
Accepted
time: 122ms
memory: 24292kb

input:

1000 10
1 158159598 758845348 641419963 949731479 892599046 605485292 272562963 71629727 807736065 544496098 895461672 992289725 360646418 47025620 104802564 747881967 584699437 779838610 638198320 671654821 729764045 195310235 459211397 159448304 78831175 237532287 408719246 520941761 546171999 897...

output:

1 415113701 475182276 672451493 793744195 950229534 659495141 726381508 631317964 412740072 453747351 272892516 598351149 977147905 739966693 418651375 174612632 596342698 77983861 481155841 428325132 962771540 518653200 864432859 639266798 599162983 114398480 12687680 351567482 742362461 888336134 ...

result:

ok 1000 numbers

Test #5:

score: -100
Time Limit Exceeded

input:

100000 483386906
1 925386004 817279817 436846237 363754513 36363197 132023702 507929492 314426859 89798592 832041984 299550534 283316534 540665653 105421174 607783420 916340928 533190624 313917727 118229052 732510050 447767998 823833186 132028675 284839462 164439702 201511962 983444172 268467355 346...

output:


result: