社区讨论

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 条回复,欢迎继续交流。

正在加载回复...