博客
关于我
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共享锁与排他锁
查看>>
MySQL内存表使用技巧
查看>>
MySQL再叙(体系结构、存储引擎、索引、SQL执行过程)
查看>>
mysql出现错误的解决办法
查看>>
MySQL函数
查看>>
mysql函数汇总之字符串函数
查看>>
mysql函数汇总之数学函数
查看>>
mysql函数汇总之日期和时间函数
查看>>
mysql函数汇总之条件判断函数
查看>>
mysql函数汇总之系统信息函数
查看>>
MySQL函数简介
查看>>
mysql函数遍历json数组
查看>>
MySQL函数(转发)
查看>>
mysql分区表
查看>>
MySQL分层架构与运行机制详解
查看>>
mysql分库分表中间件简书_MySQL分库分表
查看>>
MySQL分库分表会带来哪些问题?分库分表问题
查看>>
MySQL分组函数
查看>>
MySQL分组查询
查看>>
Mysql分表后同结构不同名称表之间复制数据以及Update语句只更新日期加减不更改时间
查看>>