【Codeforces Round#166】B. Prime Matrix【质数】

分析

题意：每次操作可以将一个数+1,问至少多少次操作能将某行/列变成全是质数。

欧拉筛存100000之内质数，然后对于每个数，二分查找比他大的第一个质数，求差值，然后看哪一行/列操作步数最少。

上代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdio>
using namespace std;

int n,m,a[510][510],v[100001],b[510][510],cnt;
int p[10001];


void pre()
{
	for(int i=2;i<=100010;i++)
	{
		if(v[i]==0)
		{
			p[++cnt]=i;
		}
		for(int j=1;j<=cnt&&i*p[j]<=100010;j++)//筛去质数倍 
		{
			v[i*p[j]]=1;
			if(i%p[j]==0) break;
		}
	}
}

int main()
{
	pre();
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&a[i][j]);
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			b[i][j]=p[lower_bound(p+1,p+cnt+1,a[i][j])-p]-a[i][j];
		}
	}
	int mnrow,mncolomn[501]={0},mn=999999999;
	for(int i=1;i<=n;i++)
	{
		mnrow=0;
		for(int j=1;j<=m;j++)
		{
			mnrow+=b[i][j];
			mncolomn[j]+=b[i][j];
		}
		mn=min(mnrow,mn);
	}
	for(int i=1;i<=m;i++)
	{
		mn=min(mn,mncolomn[i]);
	}
	cout<<mn;
	return 0;
}