码迷,mamicode.com
首页 > 其他好文 > 详细

题解 P1036 【选数】

时间:2019-03-17 23:41:43      阅读:226      评论:0      收藏:0      [点我收藏+]

标签:分析   loss   情况   scan   赋值   new   scanf   奇迹   thunder   

关于 P1036 【选数】

嗯,新手试炼场的,错了两次,对,我是蒟蒻。

因为这道题对我有帮助,所以,它是好题。

错啦两次,好尬的。

49——17——100;

不费话了,过程函数与递推。

当然要递推:

49分的不说了,从未先编译一下试试。

跟着题目走,判断质数。

来一段辣鸡代码


#include<bits/stdc++.h>

using namespace std;

int n,k;
int x[25];
int ans;
bool judge_prime(int x)
{
    for(register int i;i*i<=x;i++)
    {
        if(x%i==0)
        {
            return false;
        }
    }
    return true;
}
int rule(int choose_left_num,int already_sum,int start,int end)
{
    if(choose_left_num==0)
    {
        return judge_prime(already_sum);
    }
    int sum=0;
    for(register int i=start;i<=end;i++)
    {
        sum+=rule(choose_left_num-1,already_sum+x[i],i+1,end);
    }
    return sum;
}
int main()
{
    scanf("%d%d",&n,&k);
    for(register int i=1;i<=n;i++)
    {
        scanf("%d",&x[i]);
    }
    ans=rule(k,0,0,n-1);
    printf("%d",ans);
    return 0;
}

递推函数:

choose_left_num为剩余的k;

already_sum为前面累加的和;

start和end为全组合(常见)剩下数字的选取范围;

生成全组合,在过程中逐渐把k个数相加,当选取的数个数为0时,直接返回前面的累加和是否为质数即可;

想了好久才懂得

谢谢题解

nonetheless,

however,

whereas,

but,

就你这种辣鸡水平还想递推,一边呆着!

所谓“骗分过样例,暴力出奇迹” 这段话搜狗拼音第一个(^-^)

暴搜啊,想什么哪?

所以:


#include<bits/stdc++.h>

using namespace std;

int n,k;
bool notp[10000010];
int a[25];
void caesar(void)
{
    for(int i=2;i<10000010;i++)
    {
        if(!notp[i])
        {
            for(int j=i+i;j<10000010;j+=i)
            {
                notp[j]=1;
            }
        }
    }
    return;
}
int cae(int x,int sum,int y)
{
    if(x==0)
    {
        if(notp[sum]) return 0;
        else return 1;
    }
    if(n-y<x)return 0;
    int ans=0;
    for(register int i=y+1;i<=n;i++)
    {
        ans+=cae(x-1,sum+a[i],i);
    }
    return ans;
}
int main()
{
    caesar();
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",a+i);
    }
    printf("%d",cae(k,0,0));
    return 0;
}

老老实实地暴搜。

考虑可能需要进行多次判断,即使如此优化仍可能超时;

考虑一个更优的方法:

注意到相加得到的和最大不超过107,可以构造一个质数表;

每次检查n是否为质数只需要查表即可.如何构造一个这样的质数表呢?

由于一个质数的任意整倍数(除其自身外)都不是质数,任意一个整数都可以被表示为某(几)个质数的乘积,也就是说任意的整数都是某(几)个质数的整倍数;

因此只需将任意一个质数的二倍及以上的整倍数标记为"非质数",就可以获得这样的一张质数表了;

是这段:


int n,k;
bool notp[10000010];//我叫质数婊,呸,质数表
int a[25];
void caesar(void)//我要创造质数,额(⊙o⊙)…质数表
{
    //当前的i是质数,因此将其所有二倍及以上的整倍数都设为非质数,千万不要越界
    for(int i=2;i<10000010;i++)
    {
        if(!notp[i])
        {
            for(int j=i+i;j<10000010;j+=i)
            {
                notp[j]=1;
            }
        }
    }
    return;
}

这样一来就更优了.简单胡乱分析一波:

外层循环将重复约107次,若i不是质数则内层将仅进行1次判断,否则在判断后还将进行约(107-i)/i次赋值操作;

实验证明隐藏自己不打算推式子的事实,该函数的循环仅会运行3.95x107次,反复改变MAX的值可以发现运行次数与MAX的一次方成正比,即该方法可以在线性时间内构造质数表.

到目前为止,我们获得了初始化时间复杂度为O(n),查询时间复杂度为O(1)的质数表.下面,我们尝试回到原来的问题,看看我们还缺些什么?

我们需要从n个数中选出k个并计算它们的和.为了暴搜,我们需要寻找一种方法唯一地定义当前状态?

很容易注意到,每次只向后看是不会漏选的;

举例,如当前选到的数中下标最大的是low,那么选下一个数的时候只从下标为[low+1,n]的数中选不会导致缺失情况也不会导致重复依然不证明

因此可以重新定义当前状态:

还须选i个数,当前的和为sum,当前已经看完了下标为k及以前的数.这样就确保了每次不需要检查重复的情况,也不需要管具体选择了哪些数;

只要一一枚举k以后的数并递归调用就可以了,极大的简化了问题.

递归边界也不难得出:

当还须选0个数的时候,已经选够了,判断当前的和是否为质数,如果是返回1,否则返回0;

而非边界的时候则进行枚举,计算所有情况得到的返回值的和并返回,我们的函数返回值即使在当前限定条件下(从下标大于k的数中选i个数,在sum的基础上加上选的数的和结果是质数)的情况数;

显然在主函数中调用cae(k,0,0)得到的返回值即为结果

最后还有一个小小的可行性剪枝:若cal函数中发现n-k<i则返回0

它:


#include<bits/stdc++.h>

using namespace std;

int n,k;
bool notp[10000010];
int a[25];
void caesar(void)
{
    for(int i=2;i<10000010;i++)
    {
        if(!notp[i])
        {
            for(int j=i+i;j<10000010;j+=i)
            {
                notp[j]=1;
            }
        }
    }
    return;
}
int cae(int x,int sum,int y)
{
    if(x==0)
    {
        if(notp[sum]) return 0;
        else return 1;
    }
    if(n-y<x)return 0;
    int ans=0;
    for(register int i=y+1;i<=n;i++)
    {
        ans+=cae(x-1,sum+a[i],i);
    }
    return ans;
}
int main()
{
    caesar();
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",a+i);
    }
    printf("%d",cae(k,0,0));
    return 0;
}

差不多就这样。

谢谢!


人生总有那么一段大片大片空白的时光。你在等待,你在坚忍,你在静默。你在等一场春华秋实,你在等新一轮的春暖花开,你在等从未有过的雷霆万钧。这静默的日子有些长,有些闷,但是我也会等下去。我相信人的青春不止有一次,有时候,时光会给你额外的惊喜。

There is always a big blank space in life. You are waiting, you are persevering, you are silent. You are waiting for a spring and autumn harvest. You are waiting for a new round of spring blossoms. You are waiting for a thunderbolt that has never happened before. The silent days are long and stuffy, but I will wait. I believe that people have more than one youth. Sometimes, time gives you extra surprises.

人生にはいつも大きな空白の時間がある。あなたは待っていて、あなたは我慢して、あなたは静かにしています。あなたは春の美しさを待っていて、あなたは新しい1ラウンドの春の暖かい花が咲いて、あなたはかつてない雷の危機を待っています。この静かな日は少し長くて、少し退屈ですが、私も待っています。私は人の青春が一度だけではないと信じています。時には、時間はあなたに余分なサプライズを与えます。

题解 P1036 【选数】

标签:分析   loss   情况   scan   赋值   new   scanf   奇迹   thunder   

原文地址:https://www.cnblogs.com/XSZCaesar/p/10549502.html

(0)
(1)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!