知识小结

一般求有关素数的问题的题目，如果 数据量很大，就必须先进行预处理，对一定范围内的数进行线性素数筛选，再对其他相关问题进行处理。

比如对于Primes on Interval这一题，就必须先筛选素数，再预处理每一个区间有几个素数然后保存起来供后面使用。

而 Prime Matrix 这一题先筛选素数，在预处理每一个不是素数的数要加上几个1变成素数，把加1的个数保存起来。。真的好神奇

下面是两题的题目链接

Primes on Interval   http://acm.cugb.edu.cn/contest_showproblem.php?contestid=1036&problemid=J

Prime Matrix   http://acm.cugb.edu.cn/contest_showproblem.php?contestid=1035&problemid=J

Primes on Interval  代码：

#include <iostream>//C
#include <stdio.h>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
#include <map>
#include <stdlib.h>
#include <queue>
using namespace std;
int prim[1000005]= {0};
int pnum[1000005]= {0};
int a,b,k;
int main()
{
    int i,j;
    prim[1]=1;
    for(i=2; i<=1000005; i++)//帅选素数
    {
        for(j=2*i; j<=1000005; j+=i)
            prim[j]=1;
    }
    scanf("%d%d%d",&a,&b,&k);
    for(i=2; i<=b; i++) //先算出每一个数之前有几个素数包括它本身
    {
        pnum[i]=pnum[i-1];
        if(!prim[i])
            pnum[i]++;
    }
    if(pnum[b]-pnum[a-1]<k)
    {
        printf("-1\n");
    }
    else
    {
        int low=1,high=b-a+1,lmid,x;
        while(low<high)
        {
            int flag=0;
            lmid=(low+high)/2;
            for(x=a;x<=b-lmid+1;x++)
            {
                if(pnum[x+lmid-1]-pnum[x-1]<k)
                flag=1;
            }
            if(flag==0)
            high=lmid;
            else
            low=lmid+1;
        }
        printf("%d\n",low);

       /* int l,x;
        for( l=1; l<=b-a+1; l++)//导致超时的代码
        {
            int flag=0;
            for(x=a; x<=b-l+1; x++)
            {
                int sum=pnum[x+l-1]-pnum[x-1];
                if(sum<k)
                {
                    flag=1;
                    break;
                }
            }
            if(x==b-l+2 && flag==0)
            {
                printf("%d\n",l);
                break;
            }
        }*/
    }
    return 0;
}

Prime Matrix 代码：

#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using  namespace std;
const int M=1000010;
int p[M];
int sum[M];
int a[505][505];
int n,m;
int main()
{
    memset(sum,0,sizeof(sum));
    memset(p,0,sizeof(p));
    int i,j;
    p[1]=1;
    sum[1]=1;
    sum[0]=2;
    for(i=2; i<=M; i++)
    {
        for(j=2*i; j<=M; j+=i)
            p[j]=1;
    }
    for(i=2; i<=M; i++)
    {
        if(p[i])
        {
            for(j=i+1; j<=M; j++)
            {
                if(!p[j])
                {
                    sum[i]=j-i;//每一个数变成素数需要的步数
                    break;
                }
            }
        }
    }
    scanf("%d%d",&n,&m);
    int mmin=1000000001;
    for(i=0; i<n; i++)
    {
        for(j=0; j<m; j++)
        {
            scanf("%d",&a[i][j]);
        }
    }
    for(i=0; i<n; i++)
    {
        int s=0;
        for(j=0; j<m; j++)
        {
            s+=sum[a[i][j]];
        }
        if(s<mmin)
            mmin=s;
    }

    for(i=0; i<m; i++)
    {
        int s=0;
        for(j=0; j<n; j++)
        {
            s+=sum[a[j][i]];
        }
        if(s<mmin)
            mmin=s;
    }
    printf("%d\n",mmin);
    return 0;
}