QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#753187 | #9553. The Hermit | ucup-team5296# | AC ✓ | 179ms | 28112kb | C++20 | 2.9kb | 2024-11-16 11:41:40 | 2024-11-16 11:41:41 |
Judging History
你现在查看的是测评时间为 2024-11-16 11:41:41 的历史记录
- [2024-11-18 19:43:48]
- hack成功,自动添加数据
- (/hack/1196)
- [2024-11-16 11:41:40]
- 提交
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: 1ms
memory: 9932kb
input:
4 3
output:
7
result:
ok 1 number(s): "7"
Test #2:
score: 0
Accepted
time: 1ms
memory: 8020kb
input:
11 4
output:
1187
result:
ok 1 number(s): "1187"
Test #3:
score: 0
Accepted
time: 175ms
memory: 26764kb
input:
100000 99999
output:
17356471
result:
ok 1 number(s): "17356471"
Test #4:
score: 0
Accepted
time: 16ms
memory: 12120kb
input:
11451 1919
output:
845616153
result:
ok 1 number(s): "845616153"
Test #5:
score: 0
Accepted
time: 171ms
memory: 25988kb
input:
99998 12345
output:
936396560
result:
ok 1 number(s): "936396560"
Test #6:
score: 0
Accepted
time: 6ms
memory: 10540kb
input:
100000 1
output:
0
result:
ok 1 number(s): "0"
Test #7:
score: 0
Accepted
time: 124ms
memory: 21192kb
input:
100000 15
output:
190067060
result:
ok 1 number(s): "190067060"
Test #8:
score: 0
Accepted
time: 1ms
memory: 9972kb
input:
10 3
output:
299
result:
ok 1 number(s): "299"
Test #9:
score: 0
Accepted
time: 1ms
memory: 7900kb
input:
10 4
output:
743
result:
ok 1 number(s): "743"
Test #10:
score: 0
Accepted
time: 0ms
memory: 10088kb
input:
10 5
output:
1129
result:
ok 1 number(s): "1129"
Test #11:
score: 0
Accepted
time: 1ms
memory: 8040kb
input:
15 6
output:
28006
result:
ok 1 number(s): "28006"
Test #12:
score: 0
Accepted
time: 0ms
memory: 10108kb
input:
15 7
output:
42035
result:
ok 1 number(s): "42035"
Test #13:
score: 0
Accepted
time: 2ms
memory: 10280kb
input:
123 45
output:
214851327
result:
ok 1 number(s): "214851327"
Test #14:
score: 0
Accepted
time: 2ms
memory: 8320kb
input:
998 244
output:
964050559
result:
ok 1 number(s): "964050559"
Test #15:
score: 0
Accepted
time: 3ms
memory: 8388kb
input:
1919 810
output:
379720338
result:
ok 1 number(s): "379720338"
Test #16:
score: 0
Accepted
time: 2ms
memory: 10240kb
input:
1048 576
output:
216543264
result:
ok 1 number(s): "216543264"
Test #17:
score: 0
Accepted
time: 2ms
memory: 10388kb
input:
999 777
output:
635548531
result:
ok 1 number(s): "635548531"
Test #18:
score: 0
Accepted
time: 178ms
memory: 26388kb
input:
99999 77777
output:
448144614
result:
ok 1 number(s): "448144614"
Test #19:
score: 0
Accepted
time: 54ms
memory: 15556kb
input:
34527 6545
output:
748108997
result:
ok 1 number(s): "748108997"
Test #20:
score: 0
Accepted
time: 8ms
memory: 11248kb
input:
12345 12
output:
777496209
result:
ok 1 number(s): "777496209"
Test #21:
score: 0
Accepted
time: 1ms
memory: 10004kb
input:
1 1
output:
0
result:
ok 1 number(s): "0"
Test #22:
score: 0
Accepted
time: 170ms
memory: 26264kb
input:
100000 10101
output:
855985819
result:
ok 1 number(s): "855985819"
Test #23:
score: 0
Accepted
time: 171ms
memory: 28108kb
input:
100000 91919
output:
92446940
result:
ok 1 number(s): "92446940"
Test #24:
score: 0
Accepted
time: 174ms
memory: 28112kb
input:
100000 77979
output:
106899398
result:
ok 1 number(s): "106899398"
Test #25:
score: 0
Accepted
time: 9ms
memory: 8864kb
input:
10000 11
output:
326411649
result:
ok 1 number(s): "326411649"
Test #26:
score: 0
Accepted
time: 20ms
memory: 11320kb
input:
100000 2
output:
15322970
result:
ok 1 number(s): "15322970"
Test #27:
score: 0
Accepted
time: 26ms
memory: 11988kb
input:
100000 3
output:
93355797
result:
ok 1 number(s): "93355797"
Test #28:
score: 0
Accepted
time: 170ms
memory: 26684kb
input:
100000 99998
output:
331850772
result:
ok 1 number(s): "331850772"
Test #29:
score: 0
Accepted
time: 174ms
memory: 26236kb
input:
100000 99996
output:
885066226
result:
ok 1 number(s): "885066226"
Test #30:
score: 0
Accepted
time: 19ms
memory: 12304kb
input:
13115 2964
output:
0
result:
ok 1 number(s): "0"
Test #31:
score: 0
Accepted
time: 152ms
memory: 23032kb
input:
100000 17
output:
425792977
result:
ok 1 number(s): "425792977"
Test #32:
score: 0
Accepted
time: 145ms
memory: 21956kb
input:
99991 16
output:
667323936
result:
ok 1 number(s): "667323936"
Test #33:
score: 0
Accepted
time: 145ms
memory: 22868kb
input:
99991 17
output:
627396741
result:
ok 1 number(s): "627396741"
Test #34:
score: 0
Accepted
time: 154ms
memory: 25560kb
input:
99991 18
output:
874158501
result:
ok 1 number(s): "874158501"
Test #35:
score: 0
Accepted
time: 175ms
memory: 26652kb
input:
100000 100000
output:
99999
result:
ok 1 number(s): "99999"
Test #36:
score: 0
Accepted
time: 163ms
memory: 25684kb
input:
94229 94229
output:
94228
result:
ok 1 number(s): "94228"
Test #37:
score: 0
Accepted
time: 158ms
memory: 25816kb
input:
94229 94223
output:
476599876
result:
ok 1 number(s): "476599876"
Test #38:
score: 0
Accepted
time: 1ms
memory: 10036kb
input:
2 1
output:
0
result:
ok 1 number(s): "0"
Test #39:
score: 0
Accepted
time: 0ms
memory: 8064kb
input:
2 2
output:
0
result:
ok 1 number(s): "0"
Test #40:
score: 0
Accepted
time: 0ms
memory: 7856kb
input:
3 1
output:
0
result:
ok 1 number(s): "0"
Test #41:
score: 0
Accepted
time: 1ms
memory: 9972kb
input:
3 2
output:
2
result:
ok 1 number(s): "2"
Test #42:
score: 0
Accepted
time: 1ms
memory: 7936kb
input:
3 3
output:
2
result:
ok 1 number(s): "2"
Test #43:
score: 0
Accepted
time: 0ms
memory: 10068kb
input:
9 2
output:
44
result:
ok 1 number(s): "44"
Test #44:
score: 0
Accepted
time: 0ms
memory: 10108kb
input:
9 3
output:
206
result:
ok 1 number(s): "206"
Test #45:
score: 0
Accepted
time: 1ms
memory: 8008kb
input:
9 4
output:
441
result:
ok 1 number(s): "441"
Test #46:
score: 0
Accepted
time: 1ms
memory: 10072kb
input:
9 7
output:
224
result:
ok 1 number(s): "224"
Test #47:
score: 0
Accepted
time: 117ms
memory: 21948kb
input:
70839 22229
output:
0
result:
ok 1 number(s): "0"
Test #48:
score: 0
Accepted
time: 95ms
memory: 16560kb
input:
65536 17
output:
698801006
result:
ok 1 number(s): "698801006"
Test #49:
score: 0
Accepted
time: 92ms
memory: 18604kb
input:
65535 17
output:
433312902
result:
ok 1 number(s): "433312902"
Test #50:
score: 0
Accepted
time: 167ms
memory: 27884kb
input:
99856 317
output:
932131332
result:
ok 1 number(s): "932131332"
Test #51:
score: 0
Accepted
time: 170ms
memory: 26312kb
input:
99856 318
output:
398997854
result:
ok 1 number(s): "398997854"
Test #52:
score: 0
Accepted
time: 16ms
memory: 10968kb
input:
99856 2
output:
984791559
result:
ok 1 number(s): "984791559"
Test #53:
score: 0
Accepted
time: 179ms
memory: 24744kb
input:
100000 50000
output:
309108799
result:
ok 1 number(s): "309108799"
Extra Test:
score: 0
Extra Test Passed