QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#753187 | #9553. The Hermit | ucup-team5296# | AC ✓ | 179ms | 28224kb | C++20 | 2.9kb | 2024-11-16 11:41:40 | 2024-11-18 19:45:05 |
Judging History
answer
#include<bits/stdc++.h>
#define ll long long
#define eb emplace_back
#define ep emplace
#define pii pair<int,int>
#define fi first
#define se second
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define mems(arr,x) memset(arr,x,sizeof(arr))
#define memc(arr1,arr2) memcpy(arr1,arr2,sizeof(arr2))
using namespace std;
const int maxn=4e5+10,mod=998244353;
int n,m;
int dp[22][maxn]; // dp[i][j] 表示前 i 个数最后一个是 j 的方案数
int f[22][maxn];
namespace FastMod{
inline void madd(int &x,int y){x+=y;(x>=mod)&&(x-=mod);}
inline void mdel(int &x,int y){x-=y;(x<0)&&(x+=mod);}
inline void mmul(int &x,int y){x=1ll*x*y%mod;}
inline int imadd(int x,int y){madd(x,y);return x;}
inline int imdel(int x,int y){mdel(x,y);return x;}
inline int immul(int x,int y){mmul(x,y);return x;}
inline int qpow(int x,int y){int res=1;while(y){if(y&1) res=1ll*res*x%mod;x=1ll*x*x%mod;y>>=1;}return res;}
}
using namespace FastMod;
int mu[maxn];bool ispri[maxn];
int fac[maxn],inv[maxn];
vector<int> pri;
void init(){
int N=max(n,m);
fac[0]=1;
for(int i=1;i<=N;i++) fac[i]=immul(fac[i-1],i);
inv[N]=qpow(fac[N],mod-2);
for(int i=N-1;~i;i--) inv[i]=immul(inv[i+1],i+1);
mems(ispri,1);
mu[1]=1;
for(int i=2;i<=N;i++){
if(ispri[i]){pri.eb(i);mu[i]=-1;}
for(int j:pri){
if(i*j>n) break;
ispri[i*j]=false;
if(i%j==0) mu[i*j]=0;
else mu[i*j]=mu[i]*mu[j];
if(i%j==0) break;
}
}
}
inline int C(int x,int y){return x<y?0:1ll*fac[x]*inv[y]%mod*inv[x-y]%mod;}
int ans=0;
int main(){
scanf("%d%d",&n,&m);
init();
for(int i=1;i<=n;i++) dp[1][i]=1;
for(int i=2;i<=min(20,m-1);i++)
for(int j=1;j<=n;j++)
for(int k=j+j;k<=n;k+=j) madd(dp[i][k],dp[i-1][j]);
for(int l=m;l>=m-20-1&&l;l--){
for(int i=1;i<=n;i++){ // l 个数 gcd=i 的方案数
for(int j=i;j<=n;j+=i){
if(!mu[j/i]) continue;
if(mu[j/i]==1) madd(f[m-l][i],C(n/j,l));
else mdel(f[m-l][i],C(n/j,l));
}
// printf("%d %d %d\n",l,i,f[m-l][i]);
}
}
{
int res=C(n,m);
for(int y=1;y<=n;y++){
// int cnt=1ll*dp[i][x];
// madd(ans,immul(cnt,m-i));
mdel(res,C(n/y-1,m-1));
}
// printf("0 %d\n",res);
madd(ans,immul(res,m));
}
for(int i=1;i<=min(20,m-1);i++){
for(int x=1;x<=n;x++){
int res=C(n/x-1,m-i);
for(int y=x+x;y<=n;y+=x){
// int cnt=1ll*dp[i][x];
// madd(ans,immul(cnt,m-i));
mdel(res,C(n/y-1,m-i-1));
}
// printf("%d %d : %d %d\n",i,x,dp[i][x],res);
madd(ans,immul(dp[i][x],immul(res,m-i)));
}
}
printf("%d\n",ans);
}
这程序好像有点Bug,我给组数据试试?
詳細信息
Test #1:
score: 100
Accepted
time: 0ms
memory: 10008kb
input:
4 3
output:
7
result:
ok 1 number(s): "7"
Test #2:
score: 0
Accepted
time: 0ms
memory: 9892kb
input:
11 4
output:
1187
result:
ok 1 number(s): "1187"
Test #3:
score: 0
Accepted
time: 177ms
memory: 26536kb
input:
100000 99999
output:
17356471
result:
ok 1 number(s): "17356471"
Test #4:
score: 0
Accepted
time: 12ms
memory: 9948kb
input:
11451 1919
output:
845616153
result:
ok 1 number(s): "845616153"
Test #5:
score: 0
Accepted
time: 179ms
memory: 26572kb
input:
99998 12345
output:
936396560
result:
ok 1 number(s): "936396560"
Test #6:
score: 0
Accepted
time: 7ms
memory: 12540kb
input:
100000 1
output:
0
result:
ok 1 number(s): "0"
Test #7:
score: 0
Accepted
time: 133ms
memory: 21148kb
input:
100000 15
output:
190067060
result:
ok 1 number(s): "190067060"
Test #8:
score: 0
Accepted
time: 2ms
memory: 10004kb
input:
10 3
output:
299
result:
ok 1 number(s): "299"
Test #9:
score: 0
Accepted
time: 0ms
memory: 10024kb
input:
10 4
output:
743
result:
ok 1 number(s): "743"
Test #10:
score: 0
Accepted
time: 1ms
memory: 7988kb
input:
10 5
output:
1129
result:
ok 1 number(s): "1129"
Test #11:
score: 0
Accepted
time: 1ms
memory: 8036kb
input:
15 6
output:
28006
result:
ok 1 number(s): "28006"
Test #12:
score: 0
Accepted
time: 0ms
memory: 7924kb
input:
15 7
output:
42035
result:
ok 1 number(s): "42035"
Test #13:
score: 0
Accepted
time: 1ms
memory: 8128kb
input:
123 45
output:
214851327
result:
ok 1 number(s): "214851327"
Test #14:
score: 0
Accepted
time: 3ms
memory: 10252kb
input:
998 244
output:
964050559
result:
ok 1 number(s): "964050559"
Test #15:
score: 0
Accepted
time: 4ms
memory: 10400kb
input:
1919 810
output:
379720338
result:
ok 1 number(s): "379720338"
Test #16:
score: 0
Accepted
time: 3ms
memory: 10324kb
input:
1048 576
output:
216543264
result:
ok 1 number(s): "216543264"
Test #17:
score: 0
Accepted
time: 2ms
memory: 8264kb
input:
999 777
output:
635548531
result:
ok 1 number(s): "635548531"
Test #18:
score: 0
Accepted
time: 168ms
memory: 26360kb
input:
99999 77777
output:
448144614
result:
ok 1 number(s): "448144614"
Test #19:
score: 0
Accepted
time: 46ms
memory: 15620kb
input:
34527 6545
output:
748108997
result:
ok 1 number(s): "748108997"
Test #20:
score: 0
Accepted
time: 4ms
memory: 11144kb
input:
12345 12
output:
777496209
result:
ok 1 number(s): "777496209"
Test #21:
score: 0
Accepted
time: 0ms
memory: 10056kb
input:
1 1
output:
0
result:
ok 1 number(s): "0"
Test #22:
score: 0
Accepted
time: 168ms
memory: 26240kb
input:
100000 10101
output:
855985819
result:
ok 1 number(s): "855985819"
Test #23:
score: 0
Accepted
time: 176ms
memory: 26480kb
input:
100000 91919
output:
92446940
result:
ok 1 number(s): "92446940"
Test #24:
score: 0
Accepted
time: 171ms
memory: 26528kb
input:
100000 77979
output:
106899398
result:
ok 1 number(s): "106899398"
Test #25:
score: 0
Accepted
time: 9ms
memory: 10908kb
input:
10000 11
output:
326411649
result:
ok 1 number(s): "326411649"
Test #26:
score: 0
Accepted
time: 14ms
memory: 11272kb
input:
100000 2
output:
15322970
result:
ok 1 number(s): "15322970"
Test #27:
score: 0
Accepted
time: 26ms
memory: 13464kb
input:
100000 3
output:
93355797
result:
ok 1 number(s): "93355797"
Test #28:
score: 0
Accepted
time: 168ms
memory: 26272kb
input:
100000 99998
output:
331850772
result:
ok 1 number(s): "331850772"
Test #29:
score: 0
Accepted
time: 164ms
memory: 26320kb
input:
100000 99996
output:
885066226
result:
ok 1 number(s): "885066226"
Test #30:
score: 0
Accepted
time: 11ms
memory: 12268kb
input:
13115 2964
output:
0
result:
ok 1 number(s): "0"
Test #31:
score: 0
Accepted
time: 147ms
memory: 23112kb
input:
100000 17
output:
425792977
result:
ok 1 number(s): "425792977"
Test #32:
score: 0
Accepted
time: 136ms
memory: 22364kb
input:
99991 16
output:
667323936
result:
ok 1 number(s): "667323936"
Test #33:
score: 0
Accepted
time: 151ms
memory: 22988kb
input:
99991 17
output:
627396741
result:
ok 1 number(s): "627396741"
Test #34:
score: 0
Accepted
time: 155ms
memory: 23928kb
input:
99991 18
output:
874158501
result:
ok 1 number(s): "874158501"
Test #35:
score: 0
Accepted
time: 172ms
memory: 28224kb
input:
100000 100000
output:
99999
result:
ok 1 number(s): "99999"
Test #36:
score: 0
Accepted
time: 160ms
memory: 25704kb
input:
94229 94229
output:
94228
result:
ok 1 number(s): "94228"
Test #37:
score: 0
Accepted
time: 165ms
memory: 25644kb
input:
94229 94223
output:
476599876
result:
ok 1 number(s): "476599876"
Test #38:
score: 0
Accepted
time: 2ms
memory: 9896kb
input:
2 1
output:
0
result:
ok 1 number(s): "0"
Test #39:
score: 0
Accepted
time: 0ms
memory: 9988kb
input:
2 2
output:
0
result:
ok 1 number(s): "0"
Test #40:
score: 0
Accepted
time: 0ms
memory: 10068kb
input:
3 1
output:
0
result:
ok 1 number(s): "0"
Test #41:
score: 0
Accepted
time: 2ms
memory: 9992kb
input:
3 2
output:
2
result:
ok 1 number(s): "2"
Test #42:
score: 0
Accepted
time: 2ms
memory: 9880kb
input:
3 3
output:
2
result:
ok 1 number(s): "2"
Test #43:
score: 0
Accepted
time: 2ms
memory: 10004kb
input:
9 2
output:
44
result:
ok 1 number(s): "44"
Test #44:
score: 0
Accepted
time: 0ms
memory: 10060kb
input:
9 3
output:
206
result:
ok 1 number(s): "206"
Test #45:
score: 0
Accepted
time: 1ms
memory: 7980kb
input:
9 4
output:
441
result:
ok 1 number(s): "441"
Test #46:
score: 0
Accepted
time: 1ms
memory: 7996kb
input:
9 7
output:
224
result:
ok 1 number(s): "224"
Test #47:
score: 0
Accepted
time: 119ms
memory: 19696kb
input:
70839 22229
output:
0
result:
ok 1 number(s): "0"
Test #48:
score: 0
Accepted
time: 88ms
memory: 20264kb
input:
65536 17
output:
698801006
result:
ok 1 number(s): "698801006"
Test #49:
score: 0
Accepted
time: 93ms
memory: 18684kb
input:
65535 17
output:
433312902
result:
ok 1 number(s): "433312902"
Test #50:
score: 0
Accepted
time: 163ms
memory: 26188kb
input:
99856 317
output:
932131332
result:
ok 1 number(s): "932131332"
Test #51:
score: 0
Accepted
time: 174ms
memory: 26524kb
input:
99856 318
output:
398997854
result:
ok 1 number(s): "398997854"
Test #52:
score: 0
Accepted
time: 21ms
memory: 11204kb
input:
99856 2
output:
984791559
result:
ok 1 number(s): "984791559"
Test #53:
score: 0
Accepted
time: 167ms
memory: 24544kb
input:
100000 50000
output:
309108799
result:
ok 1 number(s): "309108799"
Extra Test:
score: 0
Extra Test Passed