bzoj3643 Phi的反函数搜索

最新推荐文章于 2018-08-15 17:36:08 发布

原创最新推荐文章于 2018-08-15 17:36:08 发布 · 982 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#函数 #phi #dfs #素数筛 #线性筛

bzoj 专栏收录该内容

346 篇文章

订阅专栏

本文介绍了一种优化算法，通过枚举质因数并搜索满足特定条件的整数，解决数学问题。利用质因数分解的性质，算法有效减少了计算复杂度，特别适用于大规模数据集。详细步骤包括定义质因数筛选、递归搜索以及特殊情况处理，最后通过代码实现验证算法的有效性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

由于不同的质因数只有11个，所以就可以搜索了。枚举x的一个质因数u，u符合条件当且仅当(u-1)|n。然后特判一下x的>n^0.5的质因数就好了。

AC代码如下：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#define ll long long
using namespace std;

int n,m,cnt,c[100005]; bool vis[100005]; ll ans;
bool prm(int x){
	int sz=(int)sqrt(x),i;
	for (i=1; c[i]<=sz; i++) if (!(x%c[i])) return 0;
	return 1;
}
void dfs(int k,int now,ll sum){
	if (sum>=ans) return; int i;
	if (now==1){ ans=sum; return; }
	if (now>m && prm(now+1)) ans=min(ans,sum*(now+1));
	for (i=k+1; c[i]-1<=m; i++){
		if (c[i]-1>now) break;
		if (!(now%(c[i]-1))){
			int x=now/(c[i]-1); ll y=sum*c[i]; dfs(i,x,y);
			while (!(x%c[i])){
				x/=c[i]; y*=c[i];
				dfs(i,x,y);
			}
		}
	}
}
int main(){
	scanf("%d",&n); m=(int)sqrt(n); int i,j;
	for (i=2; i<=50000; i++){
		if (!vis[i]) c[++cnt]=i;
		for (j=1; j<=cnt && i*c[j]<=50000; j++){
			vis[i*c[j]]=1;
			if (!(i%c[j])) break;
		}
	}
	ans=(ll)1073741824*2; dfs(1,n,1);
	if (ans<=2147483647) printf("%lld\n",ans);
		else puts("-1");
	return 0;
}

by lych

2016.2.22