昨天在USACO做了一道判断素数的题,就想着学习一下Miller_Rabin素数测试算法,在网上找到两种模版,第一种十分简洁,运行速度也很快,但是会判错极少的几个非素数;第二种比较麻烦,运行速度很慢,所以我便想找到第一种模版不能判断的非素数特判一下,结果用了一天,电脑只找到10^8以下的,10^9内还有2个没找到,但正确的模版运行速度太慢,我的电脑又太渣,耗不起时间了,姑且先这样,等以后有深入理解有更好的方法再更新一下。
第一种:源自吉林大学ACM模版
刚开始用的是随机数测试,我想到以前了解过只要用2,7,61便可能判断unsigned内的数,便做了修改,交了USACO,由于特判了5,7,11,竟然水过了,然后就愉快地去完了。。。
后来做另一道素数题,终于暴露了它的缺点,但是感觉它判断的很快,而且只有极少的不能判断,便先找出它的特例特判一下,经过一天的无聊枚举终于找到了。
文本对比还是自己写来的快,网上下的用了结果内存爆了而且运行慢,自己写的10s出结果(50mb的文本)。
目前只针对2,3,5,找到的特例是:4097,1048577,16777217,25326001,(960946321是一个10^9内的特例,这个算法枚举一次10^9,我的渣电脑都要8min左右,别说更严密的那个了)
</