hrbust/哈理工oj 1660素数矩阵【预处理+思维】

原创 2016年05月31日 20:04:54

素数矩阵

Time Limit: 2000 MS

Memory Limit: 32768 K

 

Total Submit: 120(40 users)

Total Accepted: 57(39 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

 

思路:因为只要有一行或者是一列的全部数是素数就行的话,我们假设知道了每一行中的每一个元素到达距离它最近的素数需要多少次操作,那么我们直接贪心即可,假如:

样例中第二组:

4 8 8

9 2 9

其中每个数据变成距离自己最近的素数需要的操作次数为:

1 3 3

2 0 2

这样我们就可以暴力贪心直接找最小和的那一行或者一列即可。


那么我们要如何知道这个元素距离它最近的素数需要多少次操作呢?如果每一组数据对每一个数都进行处理的话,一个是看起来麻烦,另外一个是时限问题比较尴尬,所以我们采用预处理的方式来做,首先我们使用筛法求出10^5内的全部素数,然后暴力求出10^5内的全部数据到达距离它最近的素数需要多少次操作。这样我们就很容易搞定了这个问题了~


AC代码:


#include<stdio.h>
#include<math.h>
#include<string.h>
#include<iostream>
using namespace std;
int a[505][505];
int Is_or[100050];//0表示素数
int dis[100050];
void init()
{
    memset(dis,0,sizeof(dis));
    memset(Is_or,0,sizeof(Is_or));
    Is_or[0]=1;
    Is_or[1]=1;
    int a,b;
    for(int j=2;j<sqrt(100050);j++)
    {
        if(Is_or[j]==0)
        for(int k=j+j;k<=100050;k+=j)
        {
            Is_or[k]=1;
        }
    }
    for(int i=1;i<=100050;i++)
    {
        if(Is_or[i]==0)
        {
            dis[i]=0;
        }
        else
        {
            int j=i;
            while(1)
            {
                j++;
                if(Is_or[j]==0)
                {
                    dis[i]=j-i;
                    break;
                }
            }
        }
    }
}
int main()
{
    init();
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<m;j++)
            {
                scanf("%d",&a[i][j]);
                a[i][j]=dis[a[i][j]];
            }
        }
        int output=0x3f3f3f3f;
        for(int i=0;i<n;i++)
        {
            int sum=0;
            for(int j=0;j<m;j++)
            {
                sum+=a[i][j];
            }
            output=min(output,sum);
        }
        for(int i=0;i<m;i++)
        {
            int sum=0;
            for(int j=0;j<n;j++)
            {
                sum+=a[j][i];
            }
            output=min(output,sum);
        }
        printf("%d\n",output);
    }
}








 

版权声明:0.0小白有写的不好不对的地方请指出~

2.数据的预处理

数据导入MATLAB之后,通常需要对数据进行一些预处理,例如平滑处理(去噪)、标准化变换和极差归一化变换等。 1.数据的平滑处理    1.1 smooth函数平滑处理 1.2 smoothts...
  • MATLAB_matlab
  • MATLAB_matlab
  • 2017年02月04日 08:34
  • 5233

hrbustoj新生赛

A 棋盘村A.棋盘村 Time Limit: 1000 MS Memory Limit: 32768 K Total Submit: 868 (212 users) Total Accepted:...
  • RedPolya
  • RedPolya
  • 2016年11月28日 17:10
  • 1221

hrbustoj I.行编辑器(2016级新生程序设计全国邀请赛)

Description   这次我们要写一个简单的行编辑器,当按下‘#’时代表按下了一次退格符,当按下‘@’时代表一个退行符(使当前行的字符全部无效)。例如,假设从终端接收...
  • sinat_35121480
  • sinat_35121480
  • 2016年11月26日 11:04
  • 373

Scikit-learn:数据预处理Preprocessing data

http://blog.csdn.net/pipisorry/article/details/52247679本blog内容有标准化、数据最大最小缩放处理、正则化、特征二值化和数据缺失值处理。基础知识...
  • pipisorry
  • pipisorry
  • 2016年08月19日 10:25
  • 6108

hrbust/哈理工oj 1660素数矩阵【预处理+思维】

素数矩阵 Time Limit: 2000 MS Memory Limit: 32768 K   Total Submi...
  • mengxiang000000
  • mengxiang000000
  • 2016年05月31日 20:04
  • 360

哈理工OJ 1184 早起一水(水题)

题目链接: http://acm.hrbust.edu.cn/index.php?m=ProblemSet&a=showProblem&problem_id=1184 早起一水 Time Lim...
  • qq_32866009
  • qq_32866009
  • 2016年05月17日 09:51
  • 540

hrbust 哈理工oj 2179 组合【基础dfs过】

组合 Time Limit: 1000 MS Memory Limit: 32768 K Total Submit: 86(55 users) ...
  • mengxiang000000
  • mengxiang000000
  • 2016年03月09日 17:19
  • 394

矩阵

矩阵Time Limit: 1000 ms Memory Limit: 65536 KiBSubmit StatisticProblem DescriptionX终于开始玩一些比较高大上的东西了,比如...
  • qq_38074938
  • qq_38074938
  • 2018年02月15日 22:15
  • 2

Python机器学习库SKLearn:数据集转换之预处理数据

""" 数据集转换之预处理数据: 将输入的数据转化成机器学习算法可以使用的数据。包含特征提取和标准化。 原因:数据集的标准化(服从均值为0方差为1的标准正态分布(高斯分布))是...
  • cheng9981
  • cheng9981
  • 2017年03月12日 21:47
  • 3313

哈理工oj 1079 I can do it【贪心】

哈理工oj 1079,hrbust 1079
  • mengxiang000000
  • mengxiang000000
  • 2016年03月19日 10:14
  • 708
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:hrbust/哈理工oj 1660素数矩阵【预处理+思维】
举报原因:
原因补充:

(最多只允许输入30个字)