[BZOJ]4320: ShangHai2006 Homework 分块

Description

1：在人物集合 S 中加入一个新的程序员，其代号为 X,保证 X 在当前集合中不存在。
2：在当前的人物集合中询问程序员的mod Y 最小的值。 （为什么统计这个？因为拯救
过世界的人太多了，只能取模）

题解：

考虑对值进行分块。对于每个询问的$Y$，若$Y<=\sqrt{300000}$，那么用$a[i]$表示$mod\ i$最小的余数是多少，那么询问的时候可以$O(1)$，插入的时候可以$O(\sqrt{300000})$更新；若$Y>\sqrt{300000}$，考虑枚举每个$Y$的倍数$X$，那么只要找大于等于$X$的最小值，这个最小值减去$X$就是一个可能的答案，找这个最小值也可以分块，维护两个东西，$mn[i]$表示第$i+1$一块到最后一块的最小值，$mn1[i]$表示$i$所在块内，大于等于$i$的最小值，这些都是可以在插入的时候用$O(\sqrt{300000})$的时间维护。

代码：

#include<bits/stdc++.h>
using namespace std;
#define pa pair<int,int>
#define LL long long
const int N=300000;
const int Maxn=100010;
const int inf=2147483647;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*f;
}
int bel[N+10],size;
int n,a[600],mn[600],mn1[N+10];
void ins(int x)
{
    for(int i=1;i<=size;i++)a[i]=min(a[i],x%i);
    for(int i=bel[x]-1;i;i--)mn[i]=min(mn[i],x);
    for(int i=x;;i--)
    {
        if(i==-1||(bel[i]!=bel[i+1]&&i!=x))break;
        mn1[i]=min(mn1[i],x);
    }
}
int query(int x)
{
    if(x<=size)return a[x];
    int re=inf;
    for(int i=0;i<=N;i+=x)re=min(re,min(mn[bel[i]],mn1[i])-i);
    return re;
}
int main()
{
    memset(a,63,sizeof(a));
    memset(mn,63,sizeof(mn));
    memset(mn1,63,sizeof(mn1));
    size=(int)(sqrt(N));
    bel[0]=1;for(int i=1;i<=N;i++)bel[i]=(i-1)/size+1;
    n=read();
    while(n--)
    {
        char op[3];
        scanf("%s",op);
        if(op[0]=='A')ins(read());
        else printf("%d\n",query(read()));
    }
}