【二分+思维】 HRBUST 1660 素数矩阵

素数矩阵

Time Limit: 2000 MS Memory Limit: 32768 K

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

题意

给你一个矩阵，规定若任意一行或者一列全为质数即为满足题目要求的素数矩阵。
对于矩阵内的每个数，你都可以加i直到使其变成素数。每次的花费即是你加的i。
现在问你需要最少需要多少花费使这个矩阵变成素数矩阵。

思路

其实就是找哪一行或者哪一列的所有数变成离其自身的素数最近的花费和最小，那么我们可以处理一下，把这个矩阵的全部数都变成离自己最近的素数所需的花费，因为素数是按顺序排列的，而且这个素数还一定大于等于当前这个数，我们可以用二分查找查找到这个离当前数最近的素数。

当前花费 = 离当前数最近的素数 - 当前数。

最后统计各行各列花费最小，输出即可。

坑点

无

AC代码

#include<bits/stdc++.h>
using namespace std;

#define maxn 100100
#define inf 0x3f3f3f3f
bool is_prime[maxn];
int prime[maxn];
int Map[505][505];
int n,m;
/*预处理，找到足够大足够多的质数*/
void ini(void)
{
    memset(prime,inf,sizeof prime);
    int p = 0;
    fill(is_prime,is_prime+maxn,true);
    is_prime[0] = is_prime[1] = false;
    for(int i = 2 ; i < maxn ; i++)
    {
        if(is_prime[i]){
            prime[p++] = i;
            for(int j = 2*i ; j <= maxn ; j+=i) is_prime[j] = false;
        }
    }
}

void solve(void)
{
    ini();
    while(~scanf("%d%d",&n,&m))
    {
        int hengmin = inf;
        int hengtemp = 0;
        for(int i = 0 ; i < n ; i++)
        {
            for(int j = 0 ; j < m ; j++)
            {
                int temp;
                scanf("%d",&temp);
                ///找到大于等于这个数最近的质数。直接计算出所需要的步数。
                Map[i][j] = *lower_bound(prime,prime+maxn,temp) - temp;

                hengtemp += Map[i][j];///横行花费之和。
            }
            hengmin = min(hengmin,hengtemp);///更新横行花费最小值。
            hengtemp = 0;
        }
        /*统计竖行花费最小值*/
        int zongmin = inf;
        int zongtemp = 0;
        for(int i = 0 ; i < m ; i++)
        {
            for(int j = 0 ; j < n ; j++)
            {
                zongtemp += Map[j][i];
            }
            zongmin = min(zongmin,zongtemp);
            zongtemp = 0;
        }

        printf("%d\n",min(hengmin,zongmin));
    }
}

int main(void)
{
    solve();
    return 0;
}