HRBUST OJ 1660 素数矩阵(素数打表，暴力)

题目传送门：点击打开链接

素数矩阵

Time Limit: 2000 MS Memory Limit: 32768 K Total Submit: 125(43 users) Total Accepted: 60(41 users) Rating: Special Judge: No

Description

有一个n x m 大小的矩阵，该矩阵由整数组成。对于每一次操作，可以任选择其中一个数并为这个数加 1。如果这个矩阵满足以下至少一个性质，这个矩阵就可以被称做素数矩阵： ● 矩阵中某一行元素全为质数. ● 矩阵中某一列元素全为质数. 你的任务是计算出最少的操作次数将给出的矩阵变成素数矩阵。

Input

多组测试数据。 每组测试数据第一行输入两个整数n，m(1<=n,m<=500)分别表示矩阵的行数和列数。 接下来的n行，每行有m个整数，表示矩阵的初始状态，所有的整数都是正整数，并且不超过10^5.

Output

每组数据输出一个整数表示操作的最少次数，每组输出占一行。

Sample Input

3 3 1 2 3 5 6 1 4 4 1 2 3 4 8 8 9 2 9 2 2 1 3 4 2

Sample Output

1 3 0

Author

`Wind @hrbust

   如果要节省时间的话。首先要进行素数打表，然后对矩阵进行暴力。对矩阵中的每个元素进行讨论。比如

1 2 3
5 6 1
4 4 1 

这个矩阵，处理之后变为。

1 0 0

0 1 1

5 5 1

所以这个矩阵的答案是 1

上代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
#define MM(a) memset(a,0,sizeof(a));

int prime[100005];
bool check[1000005];
int jz[505][505];
void getprime()///素数打表部分
{
int flag=0;
for(int i=2; i<=100005; i++)
{
if(!check[i])
{
prime[flag++]=i;
}
for(int j=i+i; j<=100005; j+=i)
{
check[j]=1;
}
}

}
int main()
{
MM(prime);
MM(check);

int n,m;

while(~scanf("%d %d",&n,&m))
{
MM(jz);
getprime();
for(int i=1; i<=n; i++)
{
for(int j=1; j<=m; j++)
{

scanf("%d",&jz[i][j]);///对每个元素讨论
for(int k=0;; k++)
{
if(prime[k]>=jz[i][j])
{

jz[i][j]=prime[k]-jz[i][j];

break;
}
}
}
}

int sum=0;
int sum1=88888888;
int sum2=88888888;

for(int i=1; i<=n; i++)
{
for(int j=1; j<=m; j++)
{
sum+=jz[i][j];
}
if(sum<sum1)
{
sum1=sum;

}
sum=0;
}
sum=0;
for(int i=1; i<=m; i++)
{
for(int j=1; j<=n; j++)
{
sum+=jz[j][i];
}
if(sum<sum2)
{
sum2=sum;

}
sum=0;
}
printf("%d\n",min(sum2,sum1));

}

return 0;
}