Codeforces Round #315 (Div. 2C) 568A Primes or Palindromes? 素数打表+暴力

题目：Click here

题意：π(n)表示不大于n的素数个数，rub(n)表示不大于n的回文数个数，求最大n，满足π(n) ≤ A·rub(n)。A=p/q;

分析：由于这个题A是给定范围的，所以可以先暴力求下最大的n满足上式，可以想象下随着n的增大A也在增大(总体正相关，并不是严格递增的)，所以二分查找时不行的，所以对给定的A，n是一定存在的。这个题的关键就是快速得到素数表最好在O(n)的时间以内。（杭电15多校的一个题也用到了这个算法点这里查看）

#include <bits/stdc++.h>
using namespace std;
const int M = 1179858+3;

double p, q;
bool prime[M];
int ans, pri, pal;
void getprime( )    {   // 类似筛选法求素数 打表
    prime[1] = false;
    for( int i=2; i<M; i++ )    {
        if( !prime[i] )  continue;
        prime[i] = true;
        for( int j=2; j*i<M; j++ )
            prime[i*j] = false;
    }
}

bool palindromic( int n )   {   // 判断回文数
    int fn = 0;
    int aun = n;
    while( aun )    {
        fn *= 10;
        fn += aun%10;
        aun /= 10;
    }
    return fn == n ;
}

int main()  {
    memset( prime, true, sizeof(prime) );
    getprime( );
    while( ~scanf("%lf%lf", &p, &q ) )  {
        double A = p/q;
        pri = pal = 0;
        for( int i=1; i<M; i++ )    {
            pri += prime[i];
            pal += palindromic(i);
            if( pal*p/q >= pri )
                ans = i;
        }
        printf("%d\n", ans );
    }
    return 0;
}