vijos 1042 捕风捉影——回文素数

捕风捉影

背景

《风云》里有详细记载，自己去看吧^_…

描述

有一天，雄霸传授本人风神腿法第一式：捕风捉影…的步法（弟子一：堂主，你大喘气呀。风：你给我闭嘴。）捕风捉影的关键是换气（换不好就会大喘气…）。

使用捕风捉影这一招时并不是每一步都喘气，而是在特定的步数喘气。一般来说功力越高，喘气越稀疏。喘气的步数符合特定规律：第一要是SUSHU（弟子二：哇塞！堂主，你还会鸟语，我好好崇拜你呦！可是SUSHU是什么意思呢？风：笨蛋，那是汉语拼音！）第二要是一个回文数，回文数就是正反念一样的数，如：123321，121，5211314（弟子三：堂主，最后一个好象不是…风：废话，当然不是了，我是考察一下你们的纠错能力！）现在给出两个数M，N(5<=M<N<=100,000,000)，你要算出M，N之间需要换气的都有哪几步。（包括M，N）。算出来的可以提升为本堂一级弟子，月薪（1000000000000000000000000000000000000000000 MOD 10 ）元。

格式

输入格式

两个整数M，N。用空格隔开。

输出格式

在M，N之间的换气点，每个一行。

样例1

样例输入1

100 500

样例输出1

101
131
151
181
191
313
353
373
383

题解： 见代码

/*
    乍一看数据范围很大
    其实只需要注意两点就行了
    这个回文素数必定是2n-1位 2n位回文数必定被11整除
    回文数第一位只能是1 3 7 9（如果是2、4、6、8就能被2整除不符合素数的要求）
    整个数据量不会超过10^7，再往上就没得回文数了
    数字的位数为4、6、8的回文数都能被11 整除，所有不用管（例如2992、438834)
    把回文数表打出来在判断就好了 
 */
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
typedef long long ll;
#define MAXN 10010
using namespace std;

ll n, m, cnt = 1;

ll arr[MAXN];

void read(ll &x)	// 单个字符的读取速度较快（即getchar速度比scanf快）
{
	x = 0;
	int f = 1;
	char c = getchar();
	while (c > '9' || c < '0')
	{
		if (c == '-')
			f = -1;
		c = getchar();
	}
	while ('0' <= c && c <= '9')
	{
		x = (x << 1) + (x << 3) + c - 48;
		c = getchar();
	}
	x = x * f;
}

bool sushu(ll x)
{
	for (ll i = 2; i <= sqrt(x); i++)
		if (!(x % i))
			return false;
	return true;
}

void palindromic()
{
	for (int i = 1; i <= 9; i += 2)
	{
		if (sushu(i * 10 + i))
			arr[cnt++] = i * 10 + i;
		if (i != 5)
		{
			for (int j = 0; j <= 9; j++)
			{
				if (sushu(i * 100 + j * 10 + i))
					arr[cnt++] = i * 100 + j * 10 + i;
				for (int k = 0; k <= 9; k++)
				{
					if (sushu(i * 10000 + j * 1000 + k * 100 + j * 10 + i))
						arr[cnt++] = i * 10000 + j * 1000 + k * 100 + j * 10 + i;
					for (int l = 0; l <= 9; l++)
						if (sushu(i * 1000000 + j * 100000 + k * 10000 + l * 1000 + k * 100 + j * 10 + i))
							arr[cnt++] = i * 1000000 + j * 100000 + k * 10000 + l * 1000 + k * 100 + j * 10 + i;
				}
			}
		}
	}
}

int main()
{
	read(n);
	read(m);
	// cout << n << ' ' << m;
	arr[cnt++] = 5;
	arr[cnt++] = 7;

	palindromic();
	sort(arr + 1, arr + 1 + cnt);
	for (ll i = 1; i <= cnt; i++)
	{
		if (n <= arr[i] && arr[i] <= m)
			printf("%lld\n", arr[i]);
	}
	// system("pause");
	return 0;
}

文章参考自：https://blog.csdn.net/WhiStLenA/article/details/78240180?utm_source=app

刚开始写这个题的时候傻乎乎的用质数筛选再进行判断，然后结果就是爆内存。选拔的时候还把欧拉筛法忘记了。。。幸亏队友带了埃氏筛法的板子。。。