QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#69652#5031. 核zhouhuanyi100 ✓1986ms355916kbC++112.5kb2022-12-29 14:17:022022-12-29 14:17:05

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2022-12-29 14:17:05]
  • 评测
  • 测评结果:100
  • 用时:1986ms
  • 内存:355916kb
  • [2022-12-29 14:17:02]
  • 提交

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;
}

Details

Tip: Click on the bar to expand more detailed information

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"