博客
关于我
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五补充部分:SQL逻辑查询语句执行顺序
查看>>
mysql交互式连接&非交互式连接
查看>>
MySQL什么情况下会导致索引失效
查看>>
Mysql什么时候建索引
查看>>
MySql从入门到精通
查看>>
MYSQL从入门到精通(一)
查看>>
MYSQL从入门到精通(二)
查看>>
mysql以下日期函数正确的_mysql 日期函数
查看>>
mysql以服务方式运行
查看>>
mysql优化--索引原理
查看>>
MySQL优化之BTree索引使用规则
查看>>
MySQL优化之推荐使用规范
查看>>
Webpack Critical CSS 提取与内联教程
查看>>
mysql优化概述(范式.索引.定位慢查询)
查看>>
MySQL优化的一些需要注意的地方
查看>>
mysql优化相关
查看>>
MySql优化系列-优化版造数据(存储过程+函数+修改存储引擎)-2
查看>>
MySql优化系列-进阶版造数据(load data statment)-3
查看>>
MySql优化系列-造数据(存储过程+函数)-1
查看>>
MySQL优化配置详解
查看>>