题目

P4397 [JLOI2014]聪明的燕姿

分析

根据因数和公式，显然一个数的质因数和质数都不会太多，类似反素数的套路。

所以直接搜索即可。

代码

#include<bits/stdc++.h>
using namespace std;
template <typename T>
inline void read(T &x){
	x=0;char ch=getchar();bool f=false;
	while(!isdigit(ch)){if(ch=='-'){f=true;}ch=getchar();}
	while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	x=f?-x:x;
	return ;
}
template <typename T>
inline void write(T x){
	if(x<0) putchar('-'),x=-x;
	if(x>9) write(x/10);
	putchar(x%10^48);
	return ;
}
#define ll long long
#define ull unsigned long long
#define inc(x,y,mod) (((x)+(y))>=(mod)?(x)+(y)-(mod):(x)+(y))
#define dec(x,y,mod) ((x)-(y)<0?(x)-(y)+(mod):(x)-(y))
#define rep(i,x,y) for(int i=(x);i<=(y);i++)
#define dep(i,y,x) for(int i=(y);i>=(x);i--)
const int N=1e7+5,M=5e4+5e3,INF=1e9+7;
int S,prime[M],tot,ans[N];
bool vis[M+5];
inline void GetPrimes(){
	vis[0]=vis[1]=true;
	for(int i=2;i<=M;i++){
		if(!vis[i]) prime[++prime[0]]=i;
		for(int j=1;j<=prime[0]&&prime[j]*i<=M;j++){vis[i*prime[j]]=true;if(i%prime[j]==0) break;}
	}
	return ;
}
inline bool Check(int x){
	if(x==2) return true;
	if(x%2==0) return false;
	for(int i=3;i*i<=x;i+=2) if(x%i==0)return false;
	return true;
}
void dfs(int st,int ed,int now){
	if(ed==1){ans[++tot]=now;return;}
	if(ed>prime[st]&&Check(ed-1)) ans[++tot]=now*(ed-1);
	for(int i=st;prime[i]*prime[i]<=ed;i++){
		ll t=prime[i],sum=prime[i]+1;
		for(;sum<=ed;t*=prime[i],sum+=t) if(ed%sum==0) dfs(i+1,ed/sum,now*t);
	}
	return ;
}
int main(){
	GetPrimes();
	while(scanf("%d",&S)==1){
		tot=0;dfs(1,S,1);
		write(tot);puts("");
		sort(ans+1,ans+1+tot);
		for(int i=1;i<=tot;i++) write(ans[i]),putchar(' ');
		if(tot) puts("");
	}
	return 0;
}

标签：prime,燕姿,ch,int,long,JLOI2014,mod,P4397,define
来源： https://www.cnblogs.com/Akmaey/p/15168829.html