标签:prime 燕姿 ch int long JLOI2014 mod P4397 define
题目
分析
根据因数和公式,显然一个数的质因数和质数都不会太多,类似反素数的套路。
所以直接搜索即可。
代码
#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
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。