社区讨论
50pts求调必关
P3861拆分参与者 1已保存回复 0
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 0 条
- 当前快照
- 1 份
- 快照标识符
- @mi8o2pwf
- 此快照首次捕获于
- 2025/11/21 17:37 3 个月前
- 此快照最后确认于
- 2025/11/21 19:12 3 个月前
如题
CPP#include<bits/stdc++.h>
#include<cstdio>
#define ll long long
using namespace std;
const int mod = 998244353;
const int N = 1000010;
int f[6721][6721],pos1[N],pos2[N],T,tot,pos;
ll x[N],n,sqrtn,tmp;
int read(){
int sum = 0 , f = 1;
char c = getchar();
while(c < '0' or c > '9')
{
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' and c <= '9')
{
sum = (sum << 1) + (sum << 3) + c - '0';
c = getchar();
}
return sum * f;
}
void uppp(){
while (T--){
n=read();
tot = 0,sqrtn = (ll) sqrt(n);
for (ll i=1;i*i<=n;i++){
if (n%i==0){
x[++tot]=i;
if (i*i!=n){
x[++tot]= n/i;
}
}
}
sort(x+1,x+tot+1);
for (int i=1;2*i<=tot+1;i++){
pos1[x[i]]=i;
pos2[x[i]]=tot-i+1;
}
for (int i=1;i<=tot;i++){
f[i][1]=i==1?1:0;
for (int j=2;j<=tot;j++){
f[i][j]=f[i][j-1];
if (j>i){
continue;
}
if (x[i]%x[j]==0){
tmp=x[i]/x[j];
pos=tmp<=sqrtn?pos1[tmp]:pos2[n/tmp];
f[i][j]=(f[i][j]+f[pos][j-1])%mod;
}
}
}
printf("%d\n",f[tot][tot]-1);
}
}
int main(){
T=read();
uppp();
return 0;
}
回复
共 0 条回复,欢迎继续交流。
正在加载回复...