Description   Give you a lot of positive integers, just to find out how many prime numbers there are. Input   There are a lot of cases. In each case, there is an integer N representing the number of integers to find. Each integer won’t exceed 32-bit…

HDU2138 给定N个32位大于等于2的正整数 输出其中素数的个数 用Miller Rabin 素数判定法 效率很高 数学证明比较复杂,略过, 会使用这个接口即可. #include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> using namespace std; typedef long long int LL; LL get…

hdu2138 How many prime numbers #include <bits/stdc++.h> using namespace std; typedef long long ll; ll prime[] = {,,,,}; ll qmul(ll a, ll b, ll mod) { ll res = ; while (b) { ) res = (res+a)%mod; a = (a+a)%mod; b >>= ; } return res; } ll qpow(ll…

#include<iostream> #include<cmath> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; long long mul(long long a,long long n,long long mo){ ; while (n){ ) ans=(ans+a)%mo; a=(a+a)%mo; n/=; } return ans;…

题意: 给出n个数,判断它是不是素数. SOL: 米勒拉宾裸题,思想方法略懂,并不能完全理解,所以实现只能靠背模板.... 好在不是很长... Code: /*========================================================================== # Last modified: 2016-03-21 10:09 # Filename: miller-rabin.cpp # Description: =================…

Miller_Rabin素数测试     Miller_Rabin判断单个素数的方法运用了费马小定理,可以说非常之快了.     Miller_Rabin曾经被称作“黑科技”,但是根据费马小定理其实完全可以自己写出来大半. 其算法的运行过程如下: (1)对于奇数M,使得N=(2^r)*M+1 (2)选取随机数使得A<N (3)对于任意i(i<r),若(A^(2^i)) Mod N=N - 1,则N为素数 (4)或者,若(A^M) Mod N=1,则N通过随机数A的测试 若对素数N进行T次测试,…

Prime Test Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 29046   Accepted: 7342 Case Time Limit: 4000MS Description Given a big integer number, you are required to find out whether it's a prime number. Input The first line contains the…

ll random(ll n) { return (ll)((double)rand()/RAND_MAX*n + 0.5); } ll pow_mod(ll a,ll p,ll n) { ) ; ll ans = pow_mod(a,p/,n); ans = ans*ans%n; ) ans = ans*a%n; return ans; } bool Witness(ll a,ll n) { ll m = n-; ; )) { j++; m >>= ; } ll x = pow_mod(a,…

之前一直对于这个神奇的素性判定方法感到痴迷而又没有时间去了解.借着学习<信息安全数学基础>将素性这一判定方法学习一遍. 首先证明一下费马小定理. 若p为素数,且gcd(a, p)=1, 则有 a^(p-1) = 1 (mod p) 基于以下定理 若(a, p)=1,{x| (x, p)=1}为模p下的一个完全剩余系,则{ax| (x, p)=1}也为模p下的一个完全剩余系. 又{0, 1, 2, ... p-1}为模p下一个剩余系   因此有, {a*0, a*1, a*2, ... a*(p…

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3864 题意:给定一个数n,求n的因子只有四个的情况. Miller_Rabin和Pollard_rho模板题,复杂度O(n^(1/4)),注意m^3=n的情况. //STATUS:C++_AC_62MS_232KB #include <functional> #include <algorithm> #include <iostream> //#include <ex…