8月月赛A 洛谷p3861

最新推荐文章于 2023-08-07 15:10:39 发布

MrTinTin

最新推荐文章于 2023-08-07 15:10:39 发布

阅读量454

点赞数

分类专栏：动态规划

本文链接：https://blog.csdn.net/MrTinTin/article/details/78075120

版权

动态规划专栏收录该内容

54 篇文章 0 订阅

订阅专栏

题目描述

给定一个整数 $n$ ，求将 $n$ 分解为互不相同的不小于 $2$ 的数的乘积的方案数。答案模 $998244353$ 。

输入输出格式

输入格式：

第一行一个整数 $T$ ，表示数据组数。

接下来 $T$ 行，每行一个整数 $n$ ，意义如描述所述。

输出格式：

一共 $T$ 行，每行一个整数，表示答案。

输入输出样例

输入样例#1：

1
688

输出样例#1：

说明

样例中，因为

$\times 4 \times 86= 2 \times 8 \times 43= 2 \times 344= 4 \times 172= 8 \times 86= 16 \times 43$

所以答案为 $6$

对于 $10%10\%$ 的数据，保证 $n$ 为质数

对于 $20%20\%$ 的数据，保证 $\leq n \leq 10^4$

对于 $50%50\%$ 的数据，保证 $\leq n \leq 10^7$

对于 $100%100\%$ 的数据，保证 $\leq n \leq 10^{12}$

所有数据满足 $\leq T \leq 5$

对于每一个N，我们先求出它的第 $i$ 小因子 $fac_{i}$ , 然后记录下因子 $j$ 所对应的位置 $pos_{j}$ 。这里有一个小技巧：

+ $i f$ $\le \sqrt{N}$ , 则令 $pos1_{j}=pos_{j}$

+ $i f$ $\sqrt{N}$ , 则令 $pos2Nj=posjpos2_{\frac{N}{j}}=pos_{j}$

这样就可以把空间压缩为 $O(N)O(\sqrt{N})$

然后令 $dp_{ij}$ 为把 $fac_{i}$ 分解成若干个小于等于 $fac_{j}$ 的数（包括 $1$ ）的积的方案总数，则：

$dp_{ij}=dp_{i j-1}$

$i f$ $i = j$ , $dp_{ij}++$

$i f$ $faci≡0,mod(facj)fac_{i} \equiv 0,mod( fac_{j})$ , $dpij+=dpposfacifacjj−1dp_{ij}+=dp_{pos_{\frac{fac{i}}{fac_{j}}} j-1}$

#include<iostream>
#include<cstring>
#include<map>
#include<cmath>
#include<algorithm>
#define f(i,l,r) for(i=(l);i<=(r);i++)
using namespace std;
const long long MOD=998244353,MAXN=11000;
long long n,fac[MAXN],q;
int f[MAXN][MAXN],pos1[1000011],pos2[1000010],num;
map<long long,int> ans;
int main()
{
	ios::sync_with_stdio(false);
	int T,i,j;
	cin>>T;
	while(T--){
		num=0;
		cin>>n;
		if(ans[n]){
			cout<<ans[n]<<endl;
			continue;
		}
		q=sqrt(n);
		if(n%q==0){
			fac[++num]=q;
			if(q*q!=n){
				fac[++num]=n/q;
			}
		}
		f(i,1,q-1){
			if(n%i==0){
				fac[++num]=i;
				fac[++num]=n/i;
			}
		}
		f(i,1,num){
			f(j,1,num){
				f[i][j]=0;
			}
		}
		sort(fac+1,fac+1+num);
		f(i,1,num){
			f[i][i]=1;
			if(fac[i]<=q){
				pos1[fac[i]]=i;
			}
			else{
				pos2[n/fac[i]]=i;
			}
		}
		f(i,1,num){
			f(j,1,num){
				f[i][j]+=f[i][j-1];
				if(i<=j||fac[i]%fac[j]){
				//	cout<<i<<' '<<j<<' '<<f[i][j]<<endl;
					continue;
				}
				long long tmp=fac[i]/fac[j];
				if(tmp<=q){
					f[i][j]+=f[pos1[tmp]][j-1];
				}
				else f[i][j]+=f[pos2[n/tmp]][j-1];
			//	cout<<i<<' '<<j<<' '<<f[i][j]<<endl;
				f[i][j]%=MOD;
			}
		}
		ans[n]=f[num][num]-1;
		cout<<f[num][num]-1<<endl;
	}
	return 0;
}