hdu4910 找规律+大素数检测

题意：给定一个数字，求出1 - n之间与他互质的数的乘积mod n

           规律：对于 1，2，4，P^n, 2*P^n,答案为 N-1，其余情况都是1。也就是说，1，2，4，以及只有一个质因子(奇数)或者它的1/2只有一个

                       质因子(偶数)，答案是N-1.其余情况答案均是 1

然后由于这题的n很大，也没法直接判定，可以这样搞，先去试10^6以内的素数，判断可不可以，如果不行，再利用米勒拉宾判下是否是素数，如果不是的话，把这个数字开根在平方，判断是不是完全平方数，这样做的原因是数字最大10^18,如果没有10^6以内的质因子，又不是质数的话，那么他最多只能包含2个质因子了，那么如果他不是一个完全平方数的话，那么就肯定不是了

小数版

int p[MAXN],cont[50];
bool flag[MAXN];
int k,num;
void get_prime(){
memset(flag,false,sizeof(flag));
k = 0;
for(int i = 2;i < MAXN; ++i){
if(!flag[i])
p[k++] = i;
for(int j = 0;j < k && i*p[j] < MAXN; ++j){
flag[i*p[j]] = true;
if(i % p[j] == 0) break;
}
}
}
bool cal(int n){
num = 0;
for(int i = 0;p[i] <= n && i < k; ++i){
if(n%p[i] == 0){
while(n%p[i] == 0)
n /= p[i];
num++;
if(num > 1) return false;
}
if(n == 1) break;
}
if(n > 1) num++;
if(num > 1) return false;
return true;
}


int main()
{
int n;
get_prime();
while(scanf("%d",&n) != EOF){
if(n != 4 && n % 4 == 0){
printf("1/n");
continue;
}
if(n%2 == 0){
if(cal(n/2)){
printf("%d/n",n-1);
continue;
}else{
printf("1/n");
continue;
}
}else{
if(cal(n)){
printf("%d/n",n-1);
continue;
}else{
printf("1/n");
continue;
}
}
}
}

AC代码

#include <iostream>

#include <cstdio>
#include <cstring>
#include <queue>
#include <cmath>
#include <algorithm>
#include <vector>
using namespace std;
#define PI 3.1415926535897932
#define E 2.718281828459045
#define INF 0x3f3f3f3f
#define mod 9973
#define N 70
typedef long long ll;
const int M=1005;
ll n,m;
int cnt;
int sx,sy,sz;
char mp[M];
int pa[M],rankk[M];
int head[M*6],vis[M*6];
double dis[M][10];
ll prime[M*1000];
bool isprime[M*1000];
int lowcost[M],closet[M];
char st1[5050],st2[5050];
int len[M*6];
//pair<int ,int> a[M];
vector<int> g[M];
int sum[2005];
int dp[2005][2005];
int add,rem;
struct node
{
    int st,s;
};
bool cmp(node a,node b)
{
    if(a.st!=b.st)
        return a.st<b.st;
    return a.s<b.s;
}
int month[13]= {0,31,59,90,120,151,181,212,243,273,304,334,0};
bool isrun(int year)
{
    if((year%4==0&&year%100!=0)||(year%400==0))
        return true;
    return false;
}
void getpri()
{
    ll i;
    int j;
    cnt=0;
    memset(isprime,false,sizeof(isprime));
    for(i=2; i<1000000LL; i++)
    {
        if(!isprime[i])prime[cnt++]=i;
        for(j=0; j<cnt&&prime[j]*i<1000000LL; j++)
        {
            isprime[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}
ll qk_mul(ll a,ll b,ll mo)
{
    ll t=0;
    while(b)
    {
        if(b&1)
            t=(t+a)%mo;
        a=(a<<1)%mo;
        b>>=1;
    }
    t%=mo;
    return t;
}
ll qk_mod(ll a,ll b,ll mo)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=qk_mul(ans,a,mo);
        a=qk_mul(a,a,mo);
        b>>=1;
    }
    ans%=mo;
    return ans;
}
bool miller_rabbin(ll u)
{
    if (u < 2) return false;
    if (u == 2) return true;
    if (u % 2 == 0) return false;
    for (int i = 0; i < 20; i++)
    {
        ll a = rand() % (u - 1) + 1;
        if (qk_mod(a, u - 1, u) != 1)
            return false;
    }
    return true;
    /**
        if(u==2)return true;
        if(n<2||!(n&1))return false;
        int t=0;
        ll a,x,y,v=u-1;
        while((v&1)==0)t++,v>>=1;
        for(int i=0; i<20; i++) //检测次数
        {
            a=rand()%(u-1)+1;
            x=qk_mod(a,v,u);
            for(int j=0; j<t; j++)
            {
                y=qk_mul(x,x,u);
                if(y==1&&x!=1&&x!=n-1)return false;
                ///其中用到定理，如果对模n存在1的非平凡平方根，则n是合数。
                ///如果一个数x满足方程x^2≡1 (mod n),但x不等于对模n来说1的两个‘平凡’平方根：1或-1，则x是对模n来说1的非平凡平方根
                x=y;
            }
            if(x!=1)return false;
        }
        return true;
        **/
}
bool cal(ll u)
{
    ll tmp=u;
    int num=0;
    for(int i=0; prime[i]<=tmp&&i<cnt; i++)
    {
        if(tmp%prime[i]==0)
        {
            while(tmp%prime[i]==0)
                tmp/=prime[i];
            num++;
            if(num>1)return false;
            if(tmp==1)
                return true;
            else return false;
        }


    }
    //if(tmp>1)num++;
    //if(num>1)return false;
    if(miller_rabbin(u))return true;
    ll x=(ll)sqrt(u*1.0);
    if(x*x==u)return true;
    return false;
}
int gcd(int a,int b)
{
    return b == 0 ? a : gcd(b, a%b);
}
int main()
{
    int i,j;
    char c;
    //n=50;
    int ans;
    getpri();
    while(~scanf("%I64d",&n)&&n>0)
    {
        if(n!=4&&n%4==0)
        {
            printf("1\n");
            continue;
        }
        if(n%2==0)
        {
            if(cal(n/2))
            {
                printf("%I64d\n",n-1);
                continue;
            }
            else
            {
                printf("1\n");
                continue;
            }
        }
        else
        {
            if(cal(n))
            {
                printf("%I64d\n",n-1);
                continue;
            }
            else
            {
                printf("1\n");
                continue;
            }
        }
        //printf("%d %d\n",i,ans%i);
    }
    return 0;
}