【USACO1.6.3】Prime Palindromes【数论，数学】【模拟】

版权声明：若希望转载，在评论里直接说明即可，谢谢！ https://blog.csdn.net/SSL_ZYC/article/details/83035135

题目大意：

题目链接：http://train.usaco.org/usacoprob2?a=iLZIJL4lyhX&S=pprime
求 $l$到 $r$之间的回文质数。

---

思路：

由于 $100000000$之内有 $5761455$个质数，而且我懒到不想打线性筛，所以就只能先筛回文再筛质数。
回文数很好筛。先把它的每一位储存到数组里，然后再用两个指针判断它是否回文即可。
质数的话就一边根号过去就可以了。
但是这样会T三个点。
其实，除了11以外的偶数位回文数都是11的倍数。

证明：
先设出一般形式： $an...a2a1a1a2...an$
然后可将其改写(首尾两两依次配对)：
$an...a2a1a1a2...an=an\times(10^{(2n-1)}+1)+...+a2\times(10^{(2\times2-1)}+1)\times10^{(n-2)}+a1\times(10^{(2\times1-1)}+1)\times10^{(n-1)}$
可以看到求和的每一项均有因式 $10^{(2k-1)}+1$,而该因式又含有因式 $10+1=11$,故和是 $11$的倍数.

那么就剪个枝就好了。

---

代码：

/*
ID:ssl_zyc2
TASK:pprime
LANG:C++
*/

#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

int l,r;

bool check(int x)
{
    int a[101],len=0,y=x;
    memset(a,0,sizeof(a));
    while (y)  //求出每一位
    {
        a[++len]=y%10;
        y/=10;
    }
    for (int i=1;i<=len/2;i++)
     if (a[i]!=a[len-i+1]) return 0;  //判回文
    for (int i=2;i<=sqrt(x);i++)
     if (!(x%i)) return 0;
    return 1;
}

int main()
{
    scanf("%d%d",&l,&r);
    for (register int i=l;i<=r;i++)
    {
        if (i==1000||i==100000||i==10000000) i*=10;  //不能为偶数位（两位除外）
        if (i>r) break;
        if (i>1&&check(i)) printf("%d\n",i); 
    } 
    return 0;
}