【BZOJ3629】【JLOI2014】聪明的燕姿 dfs 素数筛

链接:

#include <stdio.h>
int main()
{
    puts("转载请注明出处[vmurder]谢谢");
    puts("网址:blog.csdn.net/vmurder/article/details/44698555");
}

题解:

我们发现把一个数分解质因数以后然后可以根据每种质因数的个数算出这个数的约数和。所以我们可以暴力拆解每个数,根号时间复杂度分解。

就是枚举每种质数它用了多少,然后这个数除一下再往下一层深搜。

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 50100
using namespace std;
int prime[N],cnt;
bool vis[N];
void shake(int n)
{
    int i,j,k;
    for(i=2;i<=n;i++)
    {
        if(!vis[i])prime[++cnt]=i;
        for(j=1;j<=cnt&&i*prime[j]<=n;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}
bool check(int x)
{
    bool flag=0;
    for(int i=1;(long long)prime[i]*prime[i]<=x;i++)
        if(x%prime[i]==0)return 0;
    return 1;
}
int ans[N],num;
void dfs(long long now,int remain,int t)
{
    if(remain==1){ans[++num]=now;return ;}
    if(remain==2)return ;
    if(remain-1>=prime[t]&&check(remain-1))ans[++num]=now*(remain-1);
    for(int i=t;(long long)prime[i]*prime[i]<=remain;i++)
    {
        long long sum=prime[i]+1,x=prime[i];
        while(sum<=remain)
        {
            if(remain%sum==0)
                dfs(now*x,remain/sum,i+1);
            x*=prime[i],sum+=x;
        }
    }
}
int n;
int main()
{
//  freopen("test.in","r",stdin);

    shake(50000);
    while(scanf("%d",&n)!=EOF)
    {
        num=0;
        dfs(1,n,1);
        sort(ans+1,ans+num+1),n=0;
        for(int i=1;i<=num;i++)if(ans[i]!=ans[i-1])ans[++n]=ans[i];
        printf("%d\n",n);
        if(n>=1)printf("%d",ans[1]);
        for(int i=2;i<=n;i++)printf(" %d",ans[i]);
        if(num)puts("");
    }
    return 0;
}
时间: 2024-08-05 07:48:15

【BZOJ3629】【JLOI2014】聪明的燕姿 dfs 素数筛的相关文章

bzoj 3629: [JLOI2014]聪明的燕姿【线性筛+dfs】

数论+爆搜 详见这位大佬https://blog.csdn.net/eolv99/article/details/39644419 #include<iostream> #include<cstdio> #include<algorithm> using namespace std; const int N=100000; int s,p[N+5],tot,ans[N+5],con; bool v[N+5]; bool ok(int x) { if(x==1) retu

bzoj3629 / P4397 [JLOI2014]聪明的燕姿

P4397 [JLOI2014]聪明的燕姿 根据唯一分解定理 $n=q_{1}^{p_{1}}*q_{2}^{p_{2}}*q_{3}^{p_{3}}*......*q_{m}^{p_{m}}$ 而$n$的约数和为$\prod_{i=1}^{m} \sum_{j=0}^{p_{i}}q_{i}^j$ 于是我们可以暴搜枚举每个约数的个数,而且只要枚举到$\sqrt{S}$ tips:注意最后一个数字后不带空格 1 #include<iostream> 2 #include<cstdio&g

P4397 [JLOI2014]聪明的燕姿

P4397 [JLOI2014]聪明的燕姿 题目背景 阴天傍晚车窗外 未来有一个人在等待 向左向右向前看 爱要拐几个弯才来 我遇见谁会有怎样的对白 我等的人他在多远的未来 我听见风来自地铁和人海 我排着队拿着爱的号码牌 题目描述 城市中人们总是拿着号码牌,不停寻找,不断匹配,可是谁也不知道自己等的那个人是谁. 可是燕姿不一样,燕姿知道自己等的人是谁,因为燕姿数学学得好!燕姿发现了一个神奇的算法:假设自己的号码牌上写着数字 \(S\),那么自己等的人手上的号码牌数字的所有正约数之和必定等于 \(S

BZOJ 3629 JLOI2014 聪明的燕姿 约数和+DFS

根据约数和公式来拆s,最后再把答案乘出来,我们发先这样的话递归层数不会太大每层枚举次数也不会太多,然而我们再来个剪枝就好了 #include<cstdio> #include<iostream> #include<cmath> #include<algorithm> #include<vector> using namespace std; inline int read() { int sum=0; char ch=getchar(); whi

bzoj 3629 [JLOI2014]聪明的燕姿——约数和定理+dfs

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3629 如果要搜索,肯定得质因数分解吧:就应该朝这个方向想. **约数和定理: 对于任意一个大于1的正整数N可以分解正整数:N=P?^a? P?^a?-Pn^an,则由约数个数定理可知N的正约数有(a?+1)(a?+1)(a?+1)-(an+1)个,那么N的(a?+1)(a?+1)(a?+1)-(an+1)个正约数的和为f(N)=(P?^0+P?^1+P?^2+-P?^a?)(P?^0+P?

bzoj 3629 [JLOI2014]聪明的燕姿(约数和,搜索)

[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3629 [题意] 给定S,找出所有约数和为S的数. [思路] 若n=p1^a1*p2^a2*...*pk^ak 则约数和f(n)为(p1^0+p1+p1^2+...+p1^a1)*(p2^0+p2+p2^2+...+p2^a2)*...*(pk^0+pk+pk^2+...+pk^ak) 考虑搜索,使得和为S.至于这个搜索怎么写的,我能说我看不懂吗=_= [代码] 1 #include<

icodelab 找朋友(P4397 [JLOI2014]聪明的燕姿)

描述 老师给每个同学一个号码牌,假设小明的号码牌上写着数字 S,那么其他那些手上的号码牌数字的所有正约数之和等于 S的同学就是小明的朋友. 输入 输入包含 k 组数据. 对于每组数据,输入包含一个数字S. 输出 对于每组数据,输出有两行,第一行包含一个整数 m,表示有 m 个小明的朋友. 第二行包含相应的 m 个数,表示小明朋友的手中的数字. 注意:你输出的数字必须按照升序排列. 输入样例 1 42 输出样例 1 3 20 26 41 提示 对于 100%的数据,k≤100, S≤2×10^9

JLOI 2014--聪明的燕姿(DFS&amp;约数和定理)

写完这道题感觉人生都是灰暗的... 不存在的!!!数学学得好就知道自己等的人是谁? 这是不不不不不可能的!!!我到现在都没有女朋友!!! 题意 阴天傍晚车窗外 未来有一个人在等待 向左向右向前看 爱要拐几个弯才来 我遇见谁会有怎样的对白 我等的人他在多远的未来 我听见风来自地铁和人海 我排着队拿着爱的号码牌 城市中人们总是拿着号码牌,不停寻找,不断匹配,可是谁也不知道自己等的那个人是谁.可是燕姿不一样,燕姿知道自己等的人是谁,因为燕姿数学学得好!燕姿发现了一个神奇的算法:假设自己的号码牌上写着数

「JLOI2014」聪明的燕姿

传送门 Luogu 解题思路 很容易想到直接构造合法的数,但是这显然是会T飞的. 我们需要考虑这样一件事: 对于一个数 \(n\),对其进行质因数分解: \[n=\sum_{i=1}^x p_i^{c_i}\] 那么就会有: \[\sigma(n)=\prod_{i=1}^x \sum_{j=1}^{c_i}p^j\] 可以证明 \(\sigma(n)\) 和 \(n\) 同级,所以这个质因子 \(p_i\le \sqrt{S}\),所以我们可以直接爆搜出来所有小于 \(\sqrt{S}\) 的