luoguP1036 选数 暴力AC题解(非正解)
俗话说得好:暴力出奇迹,打表拿省一。 对于一些暴力就能拿分的题,暴力就好啦QWQ
题目描述
输入格式
输出格式
输入输出样例
定义变量
我们令输入的第一行分别为 n , k ;
第二行的数由 a [ 25 ] 来存储。
题目分析
1)制作素数筛子
看完这个题之后,我们需要用到一个判断素数的筛子。可以定义一个函数,如果是素数就返回1,否则返回0.
判断一个数是不是素数的方法也有很多种。我用的属于直观判断法。
根据定义,因为质数除了1和本身之外没有其他约数。我们知道,一个数若可以进行因数分解,那么分解时得到的两个数一定是一个小于等于sqrt(n),一个大于等于sqrt(n),据此,代码中并不需要遍历到 n-1 ,遍历到 sqrt(n) 即可,因为若sqrt(n)左侧找不到约数,那么右侧也一定找不到约数。
bool prime(long long y)
{
if (y==1||!y) return 0;
//判断1和0的情况
for (int i=2;i<=sqrt(y);i++)
if (!(y%i)) return 0;
//判断 y%i 是不是=0,如果值为0说明能被整除,不是素数
return 1; //遍历完之后如果没有返回0,则返回1.
}
这个“素数筛子”就做好了。
2)暴力循环
因为是判断 k 个数的和是不是素数,k的范围也不是特别大( 1 ≤ n ≤ 20 ,k < n )
所以,我们可以用20个 if ,从k=1开始暴力,一直到k=20。在暴力的过程中用一个计数器( tt ) ,来计算是素数的个数。
暴力也要有方法,不能无脑暴力,不然喜提TLE……
k=1时:
只需遍历一遍所有的数,看看它本身是不是素数。
此时用到了我们刚才制作的素数筛子。
分析一下:如果 a [ i ] 为素数,那么prime ( a [ i ] ) 的值就为1,if 满足条件,执行下面的 tt++ 。
相反的,如果 a [ i ] 不是素数,那么prime ( a [ i ] ) 的值就为0,if 不满足条件,什么都不执行,继续 for 循环直到 i>n。
if(k==1)
for(int i=1;i<=n;++i)
{
if(prime(a[i]))
tt++;
}
k=2时:
这时 a数组 有2个数组成最终要进行判断的数,我们可以用2层循环,把所有可能的情况都遍历一遍(暴力枚举),如果这两个数的和为素数,计数器+1.
注意:此时第二层循环的变量为第一次循环变量值+1.( int b = i + 1 , ……那里)这样可以防止出现重复判断的情况,节省了一半的时间。
另外,判断素数时,prime 括号内的部分为 a [ i ] 与 a [ b ] 之和。
if(k==2)
for(int i=1;i<=n;++i)
for(int b=i+1;b<=n;++b)
{
if(prime(a[i]+a[b]))
tt++;
}
k=3时:
同理。3个数相加,遍历一遍,不要忘记下层循环为上层+1.
在判断素数的时候也不要忘记 prime ( a [ i ] + a [ b ] + a [ c ] )。
if(k==3)
for(int i=1;i<=n;++i)
for(int b=i+1;b<=n;++b)
for(int c=b+1;c<=n;++c)
{
if(prime(a[i]+a[b]+a[c]))
tt++;
}
对!就这样!一鼓作气!打出20个 if 吧!……
AC 代码
链接:https://www.luogu.com.cn/record/35531313
瞧瞧这速度!(我想大概 也许可能 是数据有水分)
1 /*---------------------------------
2 *Title number: luoguP1036 选数
3 *Creation date: 2020-07-22 afternoon
4 *Author: EdisonBa
5 *-------------------------------*/
6 #define fastcall __attribute__((optimize("-O3")))
7 #pragma GCC optimize(2)
8 #include
9 #include
10 #include<string>
11 #include
12 #include
13 #include
14 #include
15 #include
16 #include
17 #include
403 return 0;
114514 }
这是本蒟蒻发表的第二篇题解,继承了第一篇题解的暴力传统。
这是一道橙题,我觉得打这个暴力对付它来说有点小亏。
不过也顺便锻炼了一下自己的耐力和代码能力
既然您认真地看完了,点个关注,推荐一下不香嘛!~
谢谢您的支持!
2020.8.16
EdisonBa