QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#69652 | #5031. 核 | zhouhuanyi | 100 ✓ | 1986ms | 355916kb | C++11 | 2.5kb | 2022-12-29 14:17:02 | 2022-12-29 14:17:05 |
Judging History
answer
#include<iostream>
#include<cstdio>
#define N 10000001
#define Base 100000
using namespace std;
int read()
{
char c=0;
int sum=0;
while (c<'0'||c>'9') c=getchar();
while ('0'<=c&&c<='9') sum=sum*10+c-'0',c=getchar();
return sum;
}
int n,q,invq,mod,pwa[N+1],pwb[N+1],pw[N+1],s1[N+1],s2[N+1],invs[N+1],invpw[N+1],qnum[N+1],qfac[N+1],invqfac[N+1],dp[N+1],ans;
void Adder(int &x,int d)
{
x+=d;
if (x>=mod) x-=mod;
return;
}
void Adder2(int &x,int d)
{
x+=d;
if (x<0) x+=mod;
return;
}
int MD(int x)
{
return (x>=mod)?x-mod:x;
}
int MD2(int x)
{
return (x<0)?x+mod:x;
}
int fast_pow(int a,int b)
{
int res=1,mul=a;
while (b)
{
if (b&1) res=1ll*res*mul%mod;
mul=1ll*mul*mul%mod,b>>=1;
}
return res;
}
int qC(int x,int y)
{
if (x<y) return 0;
return 1ll*qfac[x]*invqfac[y]%mod*invqfac[x-y]%mod;
}
int get_fast_pow(int x)
{
return 1ll*pwa[x%Base]*pwb[x/Base]%mod;
}
int main()
{
int cst,rst,dst=1,rest,delta,ds,d1,d2;
n=read(),q=read(),mod=read(),qfac[0]=pw[0]=invpw[0]=pwa[0]=pwb[0]=s1[0]=s2[n+2]=1,cst=MD2(q-1),rst=fast_pow(cst,mod-2),invq=fast_pow(q,mod-2),rest=1ll*cst*invq%mod;
for (int i=1;i<=Base;++i) pwa[i]=3ll*pwa[i-1]%mod;
for (int i=1;i<=Base;++i) pwb[i]=1ll*pwb[i-1]*pwa[Base]%mod;
for (int i=1;i<=n+1;++i) pw[i]=1ll*pw[i-1]*q%mod,invpw[i]=1ll*invpw[i-1]*invq%mod;
delta=1ll*MD2(pw[n+1]-pw[n])*rest%mod,ds=MD(MD2(pw[n+1]-pw[n])+rest);
for (int i=1;i<=n+1;++i) s1[i]=1ll*s1[i-1]*MD2(1-pw[i])%mod;
s2[n+2]=fast_pow(s1[n+1],mod-2);
for (int i=n+1;i>=1;--i) s2[i]=1ll*s2[i+1]*MD2(1-pw[i])%mod;
for (int i=1;i<=n+1;++i) invs[i]=1ll*s1[i-1]*s2[i+1]%mod;
for (int i=1;i<=n+1;++i) qnum[i]=1ll*MD2(pw[i]-1)*rst%mod;
for (int i=1;i<=n+1;++i) qfac[i]=1ll*qfac[i-1]*qnum[i]%mod;
invqfac[n+1]=fast_pow(qfac[n+1],mod-2);
for (int i=n;i>=0;--i) invqfac[i]=1ll*invqfac[i+1]*qnum[i+1]%mod;
dp[0]=1,dp[1]=MD2(MD2(pw[n]-pw[n-1])-1);
for (int i=2;i<=n;++i)
{
d1=MD2(MD(1ll*MD(cst+rest)*pw[i-1]%mod+1ll*MD2(pw[n]-pw[n-1])*invpw[i-1]%mod)-ds);
d2=MD2(1ll*cst*rest%mod*pw[i-2]%mod-delta);
dp[i]=1ll*MD(1ll*d1*dp[i-1]%mod+1ll*d2*dp[i-2]%mod)*invs[i]%mod;
}
for (int i=0;i<=n;++i) dp[i]=1ll*dp[i]*invqfac[n-i]%mod;
for (int i=1;i<=n;++i) dst=1ll*dst*q%(mod-1);
for (int i=0,res=1;i<=n;++i,res=1ll*res*dst%(mod-1)) Adder(ans,1ll*dp[n-i]*get_fast_pow(res)%mod);
printf("%lld\n",1ll*ans*qfac[n]%mod);
return 0;
}
詳細信息
Subtask #1:
score: 10
Accepted
Test #1:
score: 10
Accepted
time: 6ms
memory: 6348kb
input:
1 2 540053233
output:
9
result:
ok 1 number(s): "9"
Test #2:
score: 0
Accepted
time: 6ms
memory: 6500kb
input:
2 2 156542707
output:
43046970
result:
ok 1 number(s): "43046970"
Test #3:
score: 0
Accepted
time: 2ms
memory: 6408kb
input:
1 2 186225229
output:
9
result:
ok 1 number(s): "9"
Test #4:
score: 0
Accepted
time: 2ms
memory: 6636kb
input:
3 3 109884329
output:
100602209
result:
ok 1 number(s): "100602209"
Test #5:
score: 0
Accepted
time: 5ms
memory: 6388kb
input:
1 2 144802297
output:
9
result:
ok 1 number(s): "9"
Subtask #2:
score: 20
Accepted
Dependency #1:
100%
Accepted
Test #6:
score: 20
Accepted
time: 5ms
memory: 6384kb
input:
20 21992843 328859143
output:
110137213
result:
ok 1 number(s): "110137213"
Test #7:
score: 0
Accepted
time: 5ms
memory: 6632kb
input:
22 332524739 654888401
output:
410922781
result:
ok 1 number(s): "410922781"
Test #8:
score: 0
Accepted
time: 3ms
memory: 6528kb
input:
26 302215049 566649113
output:
221720840
result:
ok 1 number(s): "221720840"
Test #9:
score: 0
Accepted
time: 3ms
memory: 6452kb
input:
15 111009527 722130737
output:
648834664
result:
ok 1 number(s): "648834664"
Test #10:
score: 0
Accepted
time: 2ms
memory: 6472kb
input:
82 110032063 394529383
output:
111730592
result:
ok 1 number(s): "111730592"
Test #11:
score: 0
Accepted
time: 2ms
memory: 6396kb
input:
9 11172911 259650437
output:
68381774
result:
ok 1 number(s): "68381774"
Test #12:
score: 0
Accepted
time: 5ms
memory: 6536kb
input:
86 12016027 354886243
output:
263687778
result:
ok 1 number(s): "263687778"
Test #13:
score: 0
Accepted
time: 5ms
memory: 6500kb
input:
91 273689959 454097881
output:
114436127
result:
ok 1 number(s): "114436127"
Test #14:
score: 0
Accepted
time: 3ms
memory: 6396kb
input:
73 148878967 694206977
output:
176215101
result:
ok 1 number(s): "176215101"
Test #15:
score: 0
Accepted
time: 5ms
memory: 6328kb
input:
45 205982233 227598247
output:
156769598
result:
ok 1 number(s): "156769598"
Subtask #3:
score: 10
Accepted
Dependency #2:
100%
Accepted
Test #16:
score: 10
Accepted
time: 6ms
memory: 6436kb
input:
2778 122825869 147297463
output:
43419574
result:
ok 1 number(s): "43419574"
Test #17:
score: 0
Accepted
time: 6ms
memory: 6360kb
input:
289 7729669 589652893
output:
552952137
result:
ok 1 number(s): "552952137"
Test #18:
score: 0
Accepted
time: 6ms
memory: 6544kb
input:
2281 35651417 203950963
output:
21659018
result:
ok 1 number(s): "21659018"
Test #19:
score: 0
Accepted
time: 6ms
memory: 6584kb
input:
1684 258745639 373223677
output:
355596229
result:
ok 1 number(s): "355596229"
Test #20:
score: 0
Accepted
time: 6ms
memory: 6464kb
input:
2107 86850989 455823859
output:
245960059
result:
ok 1 number(s): "245960059"
Test #21:
score: 0
Accepted
time: 5ms
memory: 6572kb
input:
1323 43290799 791120419
output:
509649562
result:
ok 1 number(s): "509649562"
Test #22:
score: 0
Accepted
time: 3ms
memory: 6408kb
input:
2401 34064903 185314627
output:
70571452
result:
ok 1 number(s): "70571452"
Test #23:
score: 0
Accepted
time: 2ms
memory: 6500kb
input:
1073 82288187 564447959
output:
168200843
result:
ok 1 number(s): "168200843"
Test #24:
score: 0
Accepted
time: 6ms
memory: 6412kb
input:
1926 29995039 129122281
output:
60921463
result:
ok 1 number(s): "60921463"
Test #25:
score: 0
Accepted
time: 6ms
memory: 6432kb
input:
3000 66915659 765705179
output:
222619979
result:
ok 1 number(s): "222619979"
Subtask #4:
score: 20
Accepted
Test #26:
score: 20
Accepted
time: 204ms
memory: 41560kb
input:
998818 198334853 998244353
output:
153251445
result:
ok 1 number(s): "153251445"
Test #27:
score: 0
Accepted
time: 170ms
memory: 38416kb
input:
914379 128814383 998244353
output:
477606145
result:
ok 1 number(s): "477606145"
Test #28:
score: 0
Accepted
time: 185ms
memory: 37676kb
input:
944474 478445339 998244353
output:
174204073
result:
ok 1 number(s): "174204073"
Test #29:
score: 0
Accepted
time: 189ms
memory: 38052kb
input:
948637 711592127 998244353
output:
178256114
result:
ok 1 number(s): "178256114"
Test #30:
score: 0
Accepted
time: 199ms
memory: 37836kb
input:
927564 14465663 998244353
output:
315244613
result:
ok 1 number(s): "315244613"
Test #31:
score: 0
Accepted
time: 188ms
memory: 37660kb
input:
934615 392799073 998244353
output:
892700270
result:
ok 1 number(s): "892700270"
Test #32:
score: 0
Accepted
time: 182ms
memory: 36600kb
input:
917196 124972031 998244353
output:
782017412
result:
ok 1 number(s): "782017412"
Test #33:
score: 0
Accepted
time: 192ms
memory: 38104kb
input:
957149 392606173 998244353
output:
159348443
result:
ok 1 number(s): "159348443"
Test #34:
score: 0
Accepted
time: 204ms
memory: 39660kb
input:
997042 184649453 998244353
output:
464643024
result:
ok 1 number(s): "464643024"
Test #35:
score: 0
Accepted
time: 171ms
memory: 39152kb
input:
953353 14071961 998244353
output:
391688875
result:
ok 1 number(s): "391688875"
Subtask #5:
score: 10
Accepted
Test #36:
score: 10
Accepted
time: 1949ms
memory: 353196kb
input:
9909956 720431399 720431401
output:
86883659
result:
ok 1 number(s): "86883659"
Test #37:
score: 0
Accepted
time: 1886ms
memory: 353544kb
input:
9924163 267052829 267052831
output:
75754681
result:
ok 1 number(s): "75754681"
Test #38:
score: 0
Accepted
time: 1945ms
memory: 354904kb
input:
9967885 197873129 197873131
output:
16653739
result:
ok 1 number(s): "16653739"
Test #39:
score: 0
Accepted
time: 1912ms
memory: 354436kb
input:
9952642 101872151 101872153
output:
0
result:
ok 1 number(s): "0"
Test #40:
score: 0
Accepted
time: 1930ms
memory: 354528kb
input:
9955909 167874431 167874433
output:
130012020
result:
ok 1 number(s): "130012020"
Test #41:
score: 0
Accepted
time: 1888ms
memory: 355740kb
input:
9994785 399509567 399509569
output:
153324498
result:
ok 1 number(s): "153324498"
Test #42:
score: 0
Accepted
time: 1874ms
memory: 354612kb
input:
9954011 108819131 108819133
output:
101671540
result:
ok 1 number(s): "101671540"
Test #43:
score: 0
Accepted
time: 1896ms
memory: 355824kb
input:
9997570 213315827 213315829
output:
57441081
result:
ok 1 number(s): "57441081"
Test #44:
score: 0
Accepted
time: 1914ms
memory: 355916kb
input:
9995867 113028299 113028301
output:
67837072
result:
ok 1 number(s): "67837072"
Test #45:
score: 0
Accepted
time: 1908ms
memory: 353212kb
input:
9909335 247275617 247275619
output:
202966817
result:
ok 1 number(s): "202966817"
Subtask #6:
score: 30
Accepted
Dependency #1:
100%
Accepted
Dependency #2:
100%
Accepted
Dependency #3:
100%
Accepted
Dependency #4:
100%
Accepted
Dependency #5:
100%
Accepted
Test #46:
score: 30
Accepted
time: 1954ms
memory: 353596kb
input:
9921815 38466881 725310841
output:
601117286
result:
ok 1 number(s): "601117286"
Test #47:
score: 0
Accepted
time: 1928ms
memory: 353400kb
input:
9919464 4830599 747345523
output:
168521454
result:
ok 1 number(s): "168521454"
Test #48:
score: 0
Accepted
time: 1913ms
memory: 355456kb
input:
9981374 3616373 154722097
output:
2696288
result:
ok 1 number(s): "2696288"
Test #49:
score: 0
Accepted
time: 1888ms
memory: 352948kb
input:
9906664 12433457 558159149
output:
538699014
result:
ok 1 number(s): "538699014"
Test #50:
score: 0
Accepted
time: 1986ms
memory: 355372kb
input:
9985736 46853 410275823
output:
258567756
result:
ok 1 number(s): "258567756"
Test #51:
score: 0
Accepted
time: 1963ms
memory: 354744kb
input:
9962926 33790087 203505083
output:
40932778
result:
ok 1 number(s): "40932778"
Test #52:
score: 0
Accepted
time: 1866ms
memory: 352900kb
input:
9903735 146658401 157137433
output:
154493145
result:
ok 1 number(s): "154493145"
Test #53:
score: 0
Accepted
time: 1894ms
memory: 353204kb
input:
9913516 105010771 110717611
output:
67979325
result:
ok 1 number(s): "67979325"
Test #54:
score: 0
Accepted
time: 1974ms
memory: 354464kb
input:
9953517 268142489 675913921
output:
523115756
result:
ok 1 number(s): "523115756"
Test #55:
score: 0
Accepted
time: 1897ms
memory: 355372kb
input:
9981005 11993207 114120883
output:
7261617
result:
ok 1 number(s): "7261617"
Test #56:
score: 0
Accepted
time: 1920ms
memory: 354064kb
input:
9945956 36522077 168104303
output:
82398556
result:
ok 1 number(s): "82398556"
Test #57:
score: 0
Accepted
time: 1900ms
memory: 354908kb
input:
9967933 15301477 352827883
output:
242773007
result:
ok 1 number(s): "242773007"
Test #58:
score: 0
Accepted
time: 1882ms
memory: 352804kb
input:
9911781 83845891 360130933
output:
158254305
result:
ok 1 number(s): "158254305"
Test #59:
score: 0
Accepted
time: 1879ms
memory: 353108kb
input:
9916390 100404191 108138473
output:
103346432
result:
ok 1 number(s): "103346432"
Test #60:
score: 0
Accepted
time: 1939ms
memory: 355084kb
input:
9974438 7828049 430399297
output:
76675277
result:
ok 1 number(s): "76675277"