Prime Palindromes

The number 151 is a prime palindrome because it is both a prime number and a palindrome (it is the same number when read forward as backward). Write a program that finds all prime palindromes in the range of two supplied numbers a and b (5 <= a < b <= 100,000,000); both a and b are considered to be within the range .

PROGRAM NAME: pprime

INPUT FORMAT

Line 1: Two integers, a and b

SAMPLE INPUT (file pprime.in)

5 500

OUTPUT FORMAT

The list of palindromic primes in numerical order, one per line.

SAMPLE OUTPUT (file pprime.out)

5
7
11
101
131
151
181
191
313
353
373
383 ——————————————————————题解
就是枚举回文数然后筛,写了一个logn的筛法(Miller Rabin),结果乘法没开longlong导致wa了一次
枚举素数我手敲了八个循环……后来ac看analysis……只要10-9999然后翻转就行orz我感觉自己真是个辣鸡
而且这位小哥以神一般的数学直觉便判断出任何二倍长的回文数一定是11的倍数,那么我们就不用考虑偶数个的了。小学数学老师你还要我吗要我吗……
这轻描淡写的一句话顶上了我30++行代码orz
那个版权不敢盗,所以不粘别人代码,贴上我易读的代码……以及写的很丑的Miller Rabin(check函数)
 /*
PROB: pprime
LANG: C++
ID: jiaqi si
*/
#include <iostream>
#include <string.h>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
#include <ctime>
#define ivory
#define mo 1000000007
#define siji(i,x,y) for(int i=(x);i<=(y);i++)
#define gongzi(j,x,y) for(int j=(x);j>=(y);j--)
#define xiaosiji(i,x,y) for(int i=(x);i<(y);i++)
#define sigongzi(j,x,y) for(int j=(x);j>(y);j--)
#define pii pair<int,int>
#define fi first
#define se second
#define mo 1000000007
using namespace std;
int ans[],cnt;
int a,b;
int mpow(int c,int d,int k) {
if(d==) {
return c%k;
}
else {
int tmp=mpow(c,d/,k)%k;
if(d&) {
return 1LL*tmp*tmp%k*c%k;
}
else {
return 1LL*tmp*tmp%k;
}
}
}
bool _check(int v,int s,int k) {
if(v==) return true;
siji(i,,s) {
if(1LL*v*v%k==) {
if(v%k==k-){return true;}
else return false;
}
v=1LL*v*v%k;
}
return false;
}
bool check(int k) {
int tmp=k-;
int s=;
while(tmp%(<<s)==) ++s;
--s;
int d=tmp/(<<s);
int tmp1=mpow(,d,k);
int tmp2=mpow(,d,k);
int tmp3=mpow(,d,k);
bool f=;
if(!_check(tmp1,s,k)) {f=;}
else if(k> && (!_check(tmp2,s,k))) {f=;}
else if(k> && (!_check(tmp3,s,k))) {f=;}
return f;
}
int binary(int k){
int l=,r=cnt;
while(l<r) {
int mid=(l+r+)>>;
if(ans[mid]>=k) r=mid-;
else l=mid;
}
return l;
}
int main() {
#ifdef ivory
freopen("pprime.in","r",stdin);
freopen("pprime.out","w",stdout);
#else
//freopen("f1.in","r",stdin);
#endif
ans[++cnt]=;
siji(i,,) {
if(i%!= && check(i)) ans[++cnt]=i;
}
siji(i,,) {
if(i%==)continue;
int tmp=i*+i;
if(check(tmp)) ans[++cnt]=tmp;
}
siji(i,,) {
if(i%==) continue;
siji(j,,) {
int tmp=i*+j*+i;
if(check(tmp)) ans[++cnt]=tmp;
}
}
siji(i,,) {
if(i%==) continue;
siji(j,,) {
int tmp=i*+j*+j*+i;
if(check(tmp)) ans[++cnt]=tmp;
}
}
siji(i,,) {
if(i%==) continue;
siji(j,,) {
siji(k,,) {
int tmp=i*+j*+k*+j*+i;
if(check(tmp)) ans[++cnt]=tmp;
}
}
}
siji(i,,) {
if(i%==) continue;
siji(j,,) {
siji(k,,) {
int tmp=i*+j*+k*+k*+j*+i;
if(check(tmp)) ans[++cnt]=tmp;
}
}
}
siji(i,,) {
if(i%==) continue;
siji(j,,) {
siji(k,,) {
siji(h,,) {
int tmp=i*+j*+k*+h*+k*+j*+i;
if(check(tmp)) ans[++cnt]=tmp;
} }
}
}
siji(i,,) {
if(i%==) continue;
siji(j,,) {
siji(k,,) {
siji(h,,) {
int tmp=i*+j*+k*+h*+h*+k*+j*+i;
if(check(tmp)) ans[++cnt]=tmp;
}
}
}
}
scanf("%d%d",&a,&b);
int il=binary(a);
int ir=binary(b);
if(ans[ir+]==b) ir++;
siji(i,il+,ir) {
printf("%d\n",ans[i]);
}
}

 

USACO 1.5 Prime Palindromes的更多相关文章

  1. USACO Section1.5 Prime Palindromes 解题报告

    pprime解题报告 —— icedream61 博客园(转载请注明出处)--------------------------------------------------------------- ...

  2. P1217 [USACO1.5]回文质数 Prime Palindromes(求100000000内的回文素数)

    P1217 [USACO1.5]回文质数 Prime Palindromes 题目描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的),所以 151 是回文质数. 写一个程序来找 ...

  3. 洛谷 P1217 [USACO1.5]回文质数 Prime Palindromes

    P1217 [USACO1.5]回文质数 Prime Palindromes 题目描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的),所以 151 是回文质数. 写一个程序来找 ...

  4. luogu P1217 [USACO1.5]回文质数 Prime Palindromes x

    P1217 [USACO1.5]回文质数 Prime Palindromes 题目描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的),所以 151 是回文质数. 写一个程序来找 ...

  5. 4190. Prime Palindromes 一亿以内的质数回文数

    Description The number 151 is a prime palindrome because it is both a prime number and a palindrome ...

  6. <Sicily>Prime Palindromes

    一.题目描述 The number 151 is a prime palindrome because it is both a prime number and a palindrome (it i ...

  7. USACO Prime Palindromes 构造回文数

    这道题目一点也不卡素数的判断 就是朴素的sqrt(n) 也不卡 所以~放心的用吧. 构造回文的时候看了HINT 其中是这么写的: Generate palindromes by combining d ...

  8. 【USACO 1.5】Prime Palindromes

    /* TASK: pprime LANG: C++ SOLVE: 枚举数的长度,dfs出对称的数,判断是否在范围内,是否是素数 原来想着枚举每个范围里的数,但是显然超时,范围最大是10^9. 对称的数 ...

  9. USACO Section 1.5 Prime Palindromes 解题报告

    题目 题目描述 题目就是给定一个区间[a,b]((5 <= a < b <= 100,000,000)),我们需要找到这个区间内所有既是回文串又是素数的数字. 输入样例 5 500 ...

随机推荐

  1. openstack私有云布署实践【7.1 keystone + memcache (科兴环境)】

    其实登录数据库集群中任意1台都可以创建库,它们会实时自动同步数据库和对应的数据库权限.   首先登录kxcontroller1创建kx_keystone数据库,并赋于远程和本地访问的权限.   mys ...

  2. pickle.dump()

    封装是一个将Python数据对象转化为字节流的过程,拆封是封装的逆操作,将字节文件或字节对象中的字节流转化为Python数据对象,不要从不收信任的数据源中拆封数据.可以封装和拆封几乎任何Python数 ...

  3. 【SQL】T-SQL基本语法复习

    数据库基本的几个对象 数据表.视图.存储过程.索引.触发器.函数 增删改查 Insert into test(name,sex,ago) values ('陈三','男',20) Update tes ...

  4. python——爬虫

    网络爬虫(又被称为网页蜘蛛,网络机器人,在FOAF社区中间,更经常的称为网页追逐者),是一种按照一定的规则,自动地抓取万维网信息的程序或者脚本.另外一些不常使用的名字还有蚂蚁.自动索引.模拟程序或者蠕 ...

  5. OSI模型第四层传输层--TCP协议

    1.传输层2个协议tcp和udp 2.tcp的可靠性(挂号信). 面向链接的:类似寄挂号信,对方收到了并且能够确认.所以也是可靠的传输. 最大报文传输:两端可以协商传输报文大小.(协商一个报文的大小) ...

  6. PRML 第二章mindmap

    PRML第二章的Mindmap,这一章读的比较快,主要是很多计算和证明的过程都跳过了,感觉不是特别需要认真去看每一个公式,能够记住每个小节的结论.公式就可以了.当然有能力.有时间的人还是可以认真读的, ...

  7. Listview右侧 IndexBar

    qq 好友聊天界面,右侧 IndexBar  A B C D ,点击跳转到相应的联系人名字 import android.content.Context; import android.graphic ...

  8. javaee 规范技术

    J2EE的13种核心技术 一.JDBC(Java Database Connectivity) JDBC API为访问不同的数据库提供了一种统一的途径,象ODBC一样,JDBC对开发者屏蔽了一些细节问 ...

  9. android studio 工程设置项

    1.在<工程根目录>\build\generated文件夹中 存在 xxx.jar 这个文件是用来做单元测试用的,但该功能目前还处于实验阶段,不想用可以关了. 去掉勾后,clean一下工程 ...

  10. js便利json 数组的方法

    这篇文章主要介绍了JQuery遍历json数组的3种方法,本文分别给出了使用each.for遍历json的方法,其中for又分成两种形式,需要的朋友可以参考下 $(function () { var ...