计蒜客 18492.Upside down primes-米勒拉宾判大素数 (German Collegiate Programming Contest 2015 ACM-ICPC Asia Training League 暑假第一阶段第三场 K)

K. Upside down primes

传送门

这个题就是把大数按字符串输进去，判断一下是不是素数，然后反转180度，先判断反转之后的东西是不是一个数，如果是的话，再把这个数判一下是不是素数，如果都满足条件就yes。

直接调用两次米勒拉宾判大素数就可以了。

代码:

 //K-米勒拉宾判大素数
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<iomanip>
 #include<stdio.h>
 #include<stdlib.h>
 #include<math.h>
 #include<cstdlib>
 #include<set>
 #include<map>
 #include<ctime>
 #include<stack>
 #include<queue>
 #include<vector>
 #include<set>
 using namespace std;
 typedef long long ll;
 #define INF 0x7fffffff
 #define LIT 0x3f3f3f3f
 #define me(a,b) memset(a,b,sizeof(a))
 #define PI acos(-1.0)
 #define ios ios::sync_with_stdio(0),cin.tie(0);
 const int N=1e5+;
 const int maxn=1e9;
 const int S=;
 ll sw(ll a)
 {
     if(a==||a==||a==||a==||a==)
         return a;
     if(a==)
         return ;
     if(a==)
         return ;
     else
         return -;
 }
 ll mult_mod(ll a,ll b,ll mod)
 {
     a%=mod;b%=mod;
     ll ans=;
     while(b)
     {
         if(b&)
         {
             ans=ans+a;
             if(ans>=mod)
             ans=ans-mod;
         }
         a=a<<;
         if(a>=mod) a=a-mod;
         b=b>>;
     }
     return ans;
 }
 ll pow_mod(ll a,ll b,ll mod)
 {
     ll ans=;
     a=a%mod;
     while(b)
     {
         if(b&)
         {
             ans=mult_mod(ans,a,mod);
         }
         a=mult_mod(a,a,mod);
         b=b>>;
     }
     return ans;
 }
 bool check(ll a,ll n,ll x,ll t)
 {
     ll ret=pow_mod(a,x,n);
     ll last=ret;
     for(int i=;i<=t;i++)
     {
         ret=mult_mod(ret,ret,n);
         if(ret==&&last!=&&last!=n-)
             return true;
         last=ret;
     }
     if(ret!=) return true;
     else return false;
 }
 bool Miller_Rabin(ll n)
 {
     if(n<)return false;
     if(n==) return true;
     if((n&)==) return false;
     ll x=n-,t=;
     while((x&)==){x>>=;t++;}
     for(int i=;i<S;i++)
     {
         ll a=rand()%(n-)+;
         if(check(a,n,x,t))
             return false;
     }
     return true;
 }
 char ch[];
 int main()
 {
     cin>>ch;
     ll n=,bit=;
     int len=strlen(ch);
     bool flag=false;
     for(int i=len-;i>=;i--)
     {
         n+=(ch[i]-'')*bit;
         bit*=;
     }
     if(!Miller_Rabin(n))
     {
         cout<<"no"<<endl;
         return ;
     }
     n=;bit=;
     for(int i=;i<len;i++){
         ll cnt=sw(ch[i]-'');
         if(cnt==-){
             cout<<"no"<<endl;
             return ;
         }
         n+=cnt*bit;
         bit*=;
     }
     if(Miller_Rabin(n))
         cout<<"yes"<<endl;
     else
         cout<<"no"<<endl;
 }

就这些，没了。

因为这场没找到数据结构的题，就和一个队友一起刚了数论题，这场主要是图论和数论。难过，图论选手不理我，是真的回家休息去了，上一场比赛的A题和这场的图论题都比较好，如果，你能看到的话，可以去看看这些题，正好和你最近学的算法有关系。

唉，算了，都讨厌死我了。

就这些，备忘一下，米勒拉宾有点厉害。

最近在慢慢改代码风格，因为我发现队友们写代码都是把大括号对齐了写，现在再改回来。

继续去爱我的线段树去了，我滚了。