博客
关于我
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事务原理以及MVCC详解
查看>>
MySQL事务及其特性与锁机制
查看>>
mysql事务理解
查看>>
MySQL事务详解结合MVCC机制的理解
查看>>
MySQL事务隔离级别:读未提交、读已提交、可重复读和串行
查看>>
MySQL事务隔离级别:读未提交、读已提交、可重复读和串行
查看>>
webpack css文件处理
查看>>
mysql二进制包安装和遇到的问题
查看>>
MySql二进制日志的应用及恢復
查看>>
mysql互换表中两列数据方法
查看>>
mysql五补充部分:SQL逻辑查询语句执行顺序
查看>>
mysql交互式连接&非交互式连接
查看>>
MySQL什么情况下会导致索引失效
查看>>
Mysql什么时候建索引
查看>>
MySql从入门到精通
查看>>
MYSQL从入门到精通(一)
查看>>
MYSQL从入门到精通(二)
查看>>
mysql以下日期函数正确的_mysql 日期函数
查看>>