Prime Matrix

　　【B.Prime Matrix】

　　http://www.codeforces.com/contest/271/problem/B

　　题目大意：n×m的格子，每次操作可以将一个格子中的数字加1，问最少操作次数使得存在一行或一列全为质数。

　　预处理出每个格子需要多少次操作成为质数，然后求最小行、列和即可。素数筛表即可，注意多筛一些。

　　一开始没打素数表,直接+1 +1的判断,导致超时...然后就有了上一篇文章...

#include <iostream>
#include <cmath>
#define min(a,b) (((a)<(b))? (a):(b))

const  int MAXN=101000;
bool flag[MAXN];
int prime[MAXN/3],pi,num[501][501];
using namespace std;
bool isPrime(int x)
{//不打素数表的时候用来判断一个数是否为素数
    if(x==1)
        return false;
    if (x==2)
        return true;
    for (int i=2;i<=(int)sqrt((long double)x);i++)
    {
        if (x%i==0)
            return false;
    }
    return true;
}
void GetPrime_Improved()
{    //打素数表.任何一个数与其它所有素数的乘积必然也不是素数（这是因为每个合数必有一个最小素因子）。
    //初始置为false,素数为false,非素数为true
    int i,j;
    pi=0;
    memset(flag,false,sizeof(flag));
    flag[1]=true;
    for (i=2;i<MAXN;i++)
    {
        if (!flag[i])
            prime[pi++]=i;
        for ( j = 0 ; ( j<pi )&&( i * prime[j]<MAXN ) ; j++)
        {
            flag[i*prime[j]]=true;
            if ( i % prime[j] == 0) //保证每个非素数只被筛一次
                break;
        }
    }
}
int main()
{
    int i,j,k,n,m,temp,rowCount,columnCount,countMin=0x7FFFFFFF;
    cin>>n>>m;
    GetPrime_Improved();    
    for (i=0;i<n;i++)
        for (j=0;j<m;j++)
        {
            cin>>temp;
            k=0;
            while (flag[temp+k])
                k++;
            num[i][j]=k;    
        }
    for (i=0;i<n;i++)
        {
            rowCount=0;
            for (j=0;j<m;j++)
                rowCount+=num[i][j];
            countMin=min(rowCount,countMin);
        }
    for (j=0;j<m;j++)
        {
            columnCount=0;
            for (i=0;i<n;i++)
                columnCount+=num[i][j];
            countMin=min(columnCount,countMin);
        }
        cout<<countMin<<endl;
    return 0;
}