博客
关于我
bzoj 3629: [JLOI2014]聪明的燕姿
阅读量:258 次
发布时间:2019-03-01

本文共 1486 字,大约阅读时间需要 4 分钟。

求所有约数和等于S的数。

强大的dfs。

先要知道d是怎么求的。

d(n)=ki(1+p1+p2+p3+pa)
所以……直接上代码

void cal(int last,int tot,int num){    if(tot==1){ans[++n]=num;return;}    if(tot-1>sqrts&&isprime(tot-1))        ans[++n]=(tot-1)*num;    for(int i=last+1;prime[i]<=sqrts;i++)    {        int tmp=1,t=prime[i];         for (int j=1;tmp+t<=tot;j++)        {            tmp+=t;t*=prime[i];            if(tot%tmp==0)                cal(i,tot/tmp,num*(t/prime[i]));        }    }}

tot-1>sqrts是为了使选出来的数不重复,后面依次枚举 1+p,1+p+p2

最后排序,注意这题卡输出格式。
code:

#include
#include
#include
#include
#include
#include
#define LL long longusing namespace std;int prime[100010],pr=0;int sqrts,ans[1000010],n;bool v[100010];void pre(){ memset(v,true,sizeof(v)); for(int i=2;i<=100000;i++) { if(v[i]) prime[++pr]=i; for(int j=1;j<=pr&&i*prime[j]<=100000;j++) { v[i*prime[j]]=false; if(i%prime[j]==0) break; } }}bool isprime(int x){ if(x<=100000) return v[x]; int tmp=sqrt(x); for(int i=1;prime[i]<=tmp;i++) if(x%prime[i]==0) return false; return true;}void cal(int last,int tot,int num){ if(tot==1){ans[++n]=num;return;} if(tot-1>sqrts&&isprime(tot-1)) ans[++n]=(tot-1)*num; for(int i=last+1;prime[i]<=sqrts;i++) { int tmp=1,t=prime[i]; for (int j=1;tmp+t<=tot;j++) { tmp+=t;t*=prime[i]; if(tot%tmp==0) cal(i,tot/tmp,num*(t/prime[i])); } }}int main(){ int S;pre(); while(scanf("%d",&S)!=EOF) { sqrts=sqrt(S); n=0;cal(0,S,1); printf("%d\n", n); sort(ans+1,ans+n+1); for (int i=1;i

转载地址:http://gdza.baihongyu.com/

你可能感兴趣的文章
MySQL Cluster 7.0.36 发布
查看>>
Multimodal Unsupervised Image-to-Image Translation多通道无监督图像翻译
查看>>
MySQL Cluster与MGR集群实战
查看>>
multipart/form-data与application/octet-stream的区别、application/x-www-form-urlencoded
查看>>
mysql cmake 报错,MySQL云服务器应用及cmake报错解决办法
查看>>
Multiple websites on single instance of IIS
查看>>
mysql CONCAT()函数拼接有NULL
查看>>
multiprocessing.Manager 嵌套共享对象不适用于队列
查看>>
multiprocessing.pool.map 和带有两个参数的函数
查看>>
MYSQL CONCAT函数
查看>>
multiprocessing.Pool:map_async 和 imap 有什么区别?
查看>>
MySQL Connector/Net 句柄泄露
查看>>
multiprocessor(中)
查看>>
mysql CPU使用率过高的一次处理经历
查看>>
Multisim中555定时器使用技巧
查看>>
MySQL CRUD 数据表基础操作实战
查看>>
multisim变压器反馈式_穿过隔离栅供电:认识隔离式直流/ 直流偏置电源
查看>>
mysql csv import meets charset
查看>>
multivariate_normal TypeError: ufunc ‘add‘ output (typecode ‘O‘) could not be coerced to provided……
查看>>
MySQL DBA 数据库优化策略
查看>>