这道题数据范围较大,显然一般的素数筛法是不行的。那么如果要筛大质数,肯定要用到随机算法Miller_Rabin,结合费马小定理a^(p-1)≡1(mod p)

但网上的Miller_Rabin都非常繁琐,在豪哥的指导下,用他简单易懂的Miller_Rabin算法解决该题,orz。

不过感觉这个算法真的挺玄乎的,遇上一些奇怪的数据就真没法了。

#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int maxn=1e6+7;
int l,r,ans;
int Test[]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};
int Times=13;
long long ksm(long long x,long long n)
{long long base=1;long long mod=n+1;while(n){if(n&1) base=base*x%mod;n>>=1,x=(x*x)%mod;}return base%mod;
}
bool miller(int n)//简单易懂的Miller_Rabin
{if(n<=50){for(int i=0;i<=Times;i++){if(n<=Test[i])    break;if(n%Test[i]==0)    return false;}return true;}for(int i=0;i<=Times;i++){long long k=ksm(Test[i],n-1);//费马小定理判断if(k!=1) return false;}return true;
}
int main()
{//freopen("pcount.in","r",stdin);//freopen("pcount.out","w",stdout);scanf("%d%d",&l,&r);l=max(2,l);for(int i=l;i<=r;i++){ans+=miller(i);}printf("%d",ans);return 0;
}

转载于:https://www.cnblogs.com/LJB666/p/10922518.html

2019.5.25 Noip模拟测试2 T1题解相关推荐

  1. 2019.8.12 NOIP模拟测试18 反思总结

    写个博客总是符合要求的对吧 回来以后第一次悄悄参加考试,昨天全程围观- 然后喜提爆炸120分wwwwwwwww T1用了全机房最慢的写法,导致改掉死循环T掉的一个点以后还是死活过不了最后一个点.T2全 ...

  2. 2019.8.10 NOIP模拟测试16 反思总结【基本更新完毕忽视咕咕咕】

    一如既往先放代码,我还没开始改- 改完T1滚过来了,先把T1T2的题解写了[颓博客啊] 今天下午就要走了,没想到还有送行的饯别礼,真是欣喜万分[并没有] 早上刚码完前面的总结,带着不怎么有希望的心情开 ...

  3. [NOIP模拟测试34]反思+题解

    不要陷入思维定势,如果长时间没有突破就要考虑更改大方向. 不要把简单问题复杂化. 做完的题就先放下,不管能拿多少分.不能过一段时间就回来调一下. $Solutions:$ A.次芝麻 因为$n+m$始 ...

  4. [NOIP模拟测试37]反思+题解

    一定要分析清楚复杂度再打!!!窝再也不要花2h20min用暴力对拍暴力啦!!! 雨露均沾(滑稽),尽量避免孤注一掷.先把暴力分拿全再回来刚正解. 即使剩下的时间不多了也优先考虑认真读题+打暴力而非乱搞 ...

  5. NOIP 模拟赛 长寿花 题解

    NOIP 模拟赛 长寿花 题解 要放 \(n\) 层物品,第 \(i\) 层有 \(a_i\) 个位置放物品,物品有 \(m\) 中颜色,有约束条件: 同一层两个相邻物品颜色不能相同. 相邻两层颜色集 ...

  6. NOIP模拟测试19「count·dinner·chess」

    反思: 我考得最炸的一次 怎么说呢?简单的两个题0分,稍难(我还不敢说难,肯定又有人喷我)42分 前10分钟看T1,不会,觉得不可做,完全不可做,把它跳了 最后10分钟看T1,发现一个有点用的性质,仍 ...

  7. NOIP模拟测试25

    这次考试后面心态爆炸了...发现刚了2h的T2是假的之后就扔掉了,草率地打了个骗分 T1只会搜索和m=0 最先做的T3,主要是发现部分分很多,当时第一眼看上去有87分(眼瞎了). 后来想了想,感觉一条 ...

  8. 8.18 NOIP模拟测试25(B) 字符串+乌鸦喝水+所驼门王的宝藏

    T1 字符串 卡特兰数 设1为向(1,1)走,0为向(1,-1)走,限制就是不能超过$y=0$这条线,题意转化为从(0,0)出发,走到(n+m,n-m)且不越过$y=0$,然后就裸的卡特兰数,$ans ...

  9. 【2019.7.16 NOIP模拟赛 T1】洗牌(shuffle)(找环)

    找环 考虑每次洗牌其实是一次置换的过程,而这样必然就会有循环出现. 因此我们直接通过枚举找出每一个循环,询问时只要找到环上对应的位置就可以了. 貌似比我比赛时被卡成\(30\)分的倍增简单多了? 代码 ...

  10. NOIP模拟测试25「字符串·乌鸦喝水·所陀门王的宝藏(陀螺王)」

    字符串 题解 没看出catalan怎么办 dp打表啊! 考虑大力dp拿到30分好成绩!顺便收获一张表 打表发现$C_{n+m}^{m}-C_{n+m}^{m-1}$ 仔细观察然后发现其实就是之前的网格 ...

最新文章

  1. c语言间接级别不同_例9:C语言实现成绩的等级判别
  2. 网络基础知识----------iptables
  3. 报名开启!AI+科学计算专场直播:大咖齐聚×独家礼物,限量200份
  4. java二叉树合并_Java(树的前中后序遍历构造二叉树题型整合)前序和中序、中序和后序、前序和后序遍历序列构造二叉树算法整合归纳...
  5. ubuntu下安装oracle
  6. 什么是Vue.nextTick()
  7. centos mysql单向同步,虚拟机下实现Centos5.6下Mysql双向同步配置
  8. Python——杂记
  9. Java 3 年,25K 多吗?
  10. 直接拿来用!Google 新推 AI 模型打破了现有 CNN 技术壁垒
  11. js 自定义 $ 选择器
  12. 【VirtualBox】设置NAT端口映射-SSH登录
  13. 管理感悟:先做难度大的工作
  14. css 绝对定位底部居中,css – 在另一个元素的中心下方水平居中绝对定位元素
  15. xp php环境搭建 iis,Windows XP环境下快速配置IIS+PHP详细过程
  16. JPEG压缩原理详解
  17. 手机怎么将视频制作成GIF?手机端动态图片制作的方法
  18. android 色彩搭配,设计学堂:关于APP配色的一些常用色彩搭配技巧
  19. word每行的末尾箭头怎么去掉_word文档后面的箭头怎么去掉
  20. 淘宝怎么选品技巧秘籍,做淘宝7分靠选货品

热门文章

  1. Looper中的消息队列处理机制
  2. 面向对象之自动生成助记码
  3. 03K个数或第k个数算法
  4. 07堆排序 python
  5. 关于信息熵的简单理解
  6. Docker教程:dokcer machine的概念和安装
  7. 人脸检测-Haar分类器方法
  8. base cap 分布式_分布式系统一致性问题、CAP定律以及 BASE 理论
  9. java调用python脚本_python脚本估算Java线程池参数配置
  10. 力扣-102. 二叉树的层序遍历