[USACO1.5] 回文质数
Posted liylho
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[USACO1.5] 回文质数相关的知识,希望对你有一定的参考价值。
题目描述
因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的),所以 151 是回文质数。
写一个程序来找出范围[a,b](5 <= a < b <= 100,000,000)( 一亿)间的所有回文质数;
输入格式:
第 1 行: 二个整数 a 和 b .
输出格式:
输出一个回文质数的列表,一行一个。
思路:枚举一半长度的整数即可。
代码:
1 #include"bits/stdc++.h" 2 #define ll long long 3 #define ci(x) scanf("%d",&x) 4 #define cd(x) scanf("%lf",&x) 5 #define cl(x) scanf("%lld",&x) 6 #define pi(x) printf("%d\n",x) 7 #define pd(x) printf("%f\n",x) 8 #define pl(x) printf("%lld\n",x) 9 using namespace std; 10 const int N = 1e6 + 5; 11 const int mod = 1e9 + 7; 12 ll a[N],b[N]; 13 14 ll mod_mul(ll a, ll b, ll mod) { 15 ll ret = 0; 16 while(b) { 17 if(b & 1) ret = (ret + a) % mod; 18 a = (a + a) % mod; 19 b >>= 1; 20 } 21 return ret; 22 } 23 24 ll mod_exp(ll a, ll b, ll mod) { 25 ll ret = 1; 26 while(b) { 27 if(b & 1) ret = mod_mul(ret, a, mod); 28 a = mod_mul(a, a, mod); 29 b >>= 1; 30 } 31 return ret; 32 } 33 34 bool check(ll a, ll n) { 35 ll x = n - 1, y; 36 int t = 0; 37 while((x & 1) == 0) { 38 x >>= 1; 39 t++; 40 } 41 x = mod_exp(a, x, n); 42 43 for(int i = 1; i <= t; i++) { 44 y = mod_exp(x, 2, n); 45 if(y == 1 && x != 1 && x != n - 1) return true; 46 x = y; 47 } 48 if(y != 1) return true; 49 return false; 50 } 51 52 bool mill(ll n, int times = 10) { 53 if(n == 2) return true; 54 if(n == 1 || !(n & 1)) return false; 55 56 for(int i = 1; i <= times; i++) { 57 ll a = rand() % (n - 1) + 1; 58 if(check(a, n)) return false; 59 } 60 return true; 61 } 62 int p=1,q=0; 63 int n; 64 int l,r; 65 void init()//预处理出所有的回文质数 66 { 67 string f,e; 68 a[0]=11; 69 for(ll i=5;i<=10000;i++){ 70 e=to_string(i); 71 f=e; 72 reverse(f.begin(),f.end()); 73 e=e+f; 74 ll x=(ll)atoi(e.data()); 75 if(mill(x)==1) b[q++]=x; 76 77 e=to_string(i); 78 f=e.substr(0,e.size()-1); 79 reverse(f.begin(),f.end()); 80 e=e+f; 81 ll y=(ll)atoi(e.data()); 82 if(mill(y)==1) a[p++]=y; 83 } 84 int k=0; 85 for(int i=p;i<p+q;i++) a[i]=b[k++]; 86 sort(a,a+p+q); 87 n=p+q; 88 } 89 int main () 90 { 91 init(); 92 ios::sync_with_stdio(0); 93 cin>>l>>r; 94 for(int i=0;i<n;i++) if(l<=a[i]&&a[i]<=r) pl(a[i]); 95 return 0; 96 }
以上是关于[USACO1.5] 回文质数的主要内容,如果未能解决你的问题,请参考以下文章
洛谷 P1217 [USACO1.5]回文质数 Prime Palindromes
luogu P1217 [USACO1.5]回文质数 Prime Palindromes x
P1217 [USACO1.5]回文质数 Prime Palindromes
题解 P1217 [USACO1.5]回文质数 Prime Palindromes