BZOJ3136: [Baltic2013]brunhilda

最新推荐文章于 2020-01-13 20:23:43 发布

原创最新推荐文章于 2020-01-13 20:23:43 发布 · 1.4k 阅读

1 ·

CC 4.0 BY-SA版权

BZOJ 同时被 2 个专栏收录

374 篇文章

订阅专栏

138 篇文章

订阅专栏

本文介绍了一种通过递归寻找最优路径来最大程度减少操作次数的算法。该算法采用贪心策略，每次尽可能多地去除人员，并利用质数特性进行优化。通过对每个节点的状态进行预处理，实现了时间复杂度为O(n)的高效求解。

首先每次都去掉尽量多的人，可以证明这样贪心是最优的
然后递推，f[i]表示i最少操作多少次到0，j表示当前离i最远的i-i%p[x]的点，f[i]=f[j]+1，如果j不满足条件了就要向后找j
因为f[i]是单调不降的，所以我们需要找最小的j使得j能对i产生贡献
j能对当前的i产生贡献，易知j一定是题目给定某一个的质数p[x]的倍数，需要知道j的质因子里面含题目给定的质数的最大的那个，预处理出每个数最小的质因数minp，对于i含的最大的p[x]的只要看minp和i/minp

时间复杂度O(n)

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 10001000;
const int maxm = 101000;

int minp[maxn],pr[maxn/10],pri;
int v[maxn];
int pi[maxm],p[maxm],f[maxn],n,m,up;
int q[maxm];

void pre()
{
    pri=0; minp[1]=1;
    for(int i=2;i<up;i++)
    {
        if(!v[i])
        {
            pr[++pri]=i;
            minp[i]=i;
        }
        for(int j=1;j<=pri;j++)
        {
            int k=pr[j]*i;
            if(k>=up) break;
            v[k]=1;
            minp[k]=pr[j];
            if(i%pr[j]==0) break;
        }
    }
}

int main()
{
    up=0;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&pi[i]);
    for(int i=1;i<=m;i++) scanf("%d",&q[i]),up=max(up,q[i]+1);

    pre();
    memset(v,0,sizeof v);

    sort(pi+1,pi+n+1);
    int nn=0;
    for(int i=1;i<=n;i++) 
        if(pi[i]!=pi[i-1]) 
            p[++nn]=pi[i],v[p[nn]]=p[nn];
    n=nn;
    v[0]=p[n];

    f[0]=0;
    int j=0;
    for(int i=1;i<up;i++)
    {
        while(i-j>=v[j]||!v[j])
        {
            j++;
            if(i==j) break;
        }
        if(i==j) {up=i; break;}
        f[i]=f[j]+1;
        v[i]=max(v[minp[i]],v[i/minp[i]]);
    }
    for(int i=1;i<=m;i++)
    {
        int x=q[i];
        if(x>=up) printf("oo\n");
        else printf("%d\n",f[x]);
    }

    return 0;
}