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;
}