Codeforces Round #166 (Div. 2) B.Prime Matrix【二分】【欧拉筛】

JA_yichao

已于 2022-01-21 16:44:44 修改

阅读量246

点赞数

CC 4.0 BY-SA版权

分类专栏：题解 nowcoder&&CodeForces 前缀和文章标签： c++ 题解前缀和暴力 Codeforces

于 2022-01-21 16:44:17 首次发布

本文链接：https://blog.csdn.net/Jackma_mayichao/article/details/122624781

题解同时被 3 个专栏收录

414 篇文章

订阅专栏

前缀和

14 篇文章

订阅专栏

nowcoder&&CodeForces

6 篇文章

订阅专栏

本文介绍了一种算法，利用欧拉筛法找到质数，并通过二分查找法计算每个元素到下一个较大质数的距离，最终求得每行每列距离之和的最小值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目

在这里插入图片描述

思路

直接用欧拉筛筛出质数
然后二分把每个位置到离第一个比它大的质数的距离处理出来
累加求min即可

代码

#include<iostream>
#include<cstdio>
using namespace std;
int a[510][510],js[510][510];
int prime[1000010],v[1000010];
int n,m,c,sum,ans=1e9;
void priming()
{
	for(int i=2; i<=200000; i++)
	 {
	 	if(v[i]==0)
	 	  prime[++c]=i;
	 	for(int j=1; i*prime[j]<=200000&&j<=c; j++)
	 	 {
	 	 	v[i*prime[j]]=1;
	 	 	if(i%prime[j]==0)
	 	 	  break;
		 }
	 }
}
int middle(int x)
{
	int mid,l=1,r=c,news=1e9;
	while(l<=r)
	 {
	 	mid=(l+r)>>1;
	 	if(prime[mid]>=x)
	 	  r=mid-1,news=min(news,prime[mid]);
	 	else
	 	  l=mid+1;
	 }
	return news-x;
}
int main()
{
	priming();
	scanf("%d%d",&n,&m);
	for(int i=1; i<=n; i++)
	 for(int j=1; j<=m; j++)
	  {
	    scanf("%d",&a[i][j]);
	    js[i][j]=middle(a[i][j]);
	  }
	for(int i=1; i<=n; i++)
	 {
	 	sum=0;
		for(int j=1; j<=m; j++)
		   sum+=js[i][j];
		ans=min(ans,sum);
	 }
	for(int j=1; j<=m; j++)
	 {
	 	sum=0;
		for(int i=1; i<=n; i++)
		   sum+=js[i][j];
		ans=min(ans,sum);
	 }
	printf("%d",ans);
	return 0;
}