2023集dfs 1 反素数

最新推荐文章于 2024-09-27 14:28:55 发布

renxvqing0522

最新推荐文章于 2024-09-27 14:28:55 发布

阅读量32

点赞数

文章标签：深度优先算法图论

本文链接：https://blog.csdn.net/2201_75785857/article/details/132038213

版权

[POI2001] [HAOI2007] 反素数

题目描述

对于任何正整数 $x$ ，其约数的个数记作 $g (x)$ 。例如 $g (1) = 1$ ， $g (6) = 4$ 。

如果某个正整数 $x$ 满足： $\forall 0 \lt i \lt x$ ，都有 $\gt g(i)$ ，则称 $x$ 为反质数。例如，整数 $1, 2, 4, 6$ 等都是反质数。

现在给定一个数 $N$ ，你能求出不超过 $N$ 的最大的反质数么？

输入格式

一个数 $N$ 。

输出格式

不超过 $N$ 的最大的反质数。

样例 #1

样例输入 #1

样例输出 #1

提示

$\leq N \leq 2 \times 10^9$

思路

约数个数

假设x= p1^k1* p2^k2* …pm^km,dx为x的约数，dx= p1^a1 p2^a2* …* pm^km
那么一定能得a1<=k1,a2<=k2,…,am<=km
所以dx有p1^0, p1^1 ,…,p1^k1共k1+1种可能
由此得x= p1^k1* p2^k2* …pi^ki,g(x)=(k1+1)(k2+1)…(ki+1)
例如12=2^2*3，g(12)=(2+1)(1+1)=6

剪枝

性质：pi保持单调递增，ki保持不上升
证明：即x= p1^k1* p2^k2* …pm^km中，p1<p2<…<pm ,k1>=k2>=…>=km
设k1<k2,xx=p1^k2 * p2^k1 … *pm^km
那么xx=(p1p2)^k1 * p1^(k2-k1) …
又因为 x=(p1p2)^k1 * p2^(k2-k1) …
可知xx<x且约数个数相等，所以x必然不是反素数
例：18=2 * 3 * 3 g(18)=2 * 3=6 12=2 * 2 * 3,g(12)=2 * 3=6
所以18不是反素数

完整代码

#include<bits/stdc++.h>
using namespace std;
int n,tot;//tot是可能成为反素数的数的个数
pair<int,int> a[10005];//第一个存储的是数，第二个存储的是第一个数所对应的约数
int p[]={0,2,3,5,7,11,13,17,19,23};//这里第一个必须要是0，用来站位 ，
//因为数据范围小，把这几个数乘起来，就超范围了
void dfs(int i,int lst,int x,int g)//i表示当前是第几个质数，
//lst表示上一个质数的指数，x表示当前的数，g表示当前的数有多少个约数
{
	a[++tot]=make_pair(x,g);//读入这个数，待会枚举
	if(i>9)	return;//一共九个质数
	for(int j=1;j<=lst;j++)// 剪枝
	{
		if(1ll*p[i]*x>n)break;
		x*=p[i];
		dfs(i+1,j,x,g*(j+1));
	}
}
int main()
{
	cin>>n;
	dfs(1,100,1,1);
	int ans=0,mx=0;
	sort(a+1,a+1+tot); //一定要排序 
	for(int i=1;i<=tot;i++)//筛选
	{
		if(a[i].second>mx)
		{
			mx=a[i].second;
			ans=a[i].first;
		}
	}
	cout<<ans;
	return 0;
}