BZOJ 4070 [Apio2015]雅加达的摩天楼 ——分块 SPFA

挺有趣的分块的题目。

直接暴力建边SPFA貌似是$O(nm)$的。

然后考虑分块，$sqrt n$一下用虚拟节点辅助连边，

以上的直接暴力连边即可。

然后卡卡时间，卡卡空间。

终于在UOJ上T掉辣。

遂弃疗

#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define ll long long
#define mp make_pair
#define inf 0x3f3f3f3f
#define maxn 30005
 
int n,m,b[maxn],p[maxn],cnt=0;
int num[maxn][101],T,s,t;
int h[maxn*101],to[maxn*500],ne[maxn*500],en=0,dis[maxn*101];
int w[maxn*500],inq[maxn*101];
queue <int> q;
void add(int a,int b,int c)
{to[en]=b;ne[en]=h[a];w[en]=c;h[a]=en++;}
 
void SPFA()
{
    memset(dis,0x3f,sizeof dis);
    dis[s]=0;inq[s]=1;q.push(s);
    while (!q.empty())
    {
        int x=q.front();q.pop();inq[x]=0;
        for (int i=h[x];i>=0;i=ne[i])
            if (dis[to[i]]>dis[x]+w[i])
            {
                dis[to[i]]=dis[x]+w[i];
                if (!inq[to[i]])
                {
                    inq[to[i]]=1;
                    q.push(to[i]);
                }
            }
    }
    if (dis[t]!=inf) printf("%d
",dis[t]);
    else printf("-1
");
    return;
}
 
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*10+ch-'0'; ch=getchar();}
    return x*f;
}
 
int main()
{
    memset(h,-1,sizeof h);
    n=read();m=read();T=min((int)sqrt(n),100);cnt=n;
    F(i,1,n) F(j,1,T) num[i][j]=++cnt;
    F(i,1,m) b[i]=read(),p[i]=read(),b[i]++;
    s=b[1];t=b[2];
    F(i,1,n) F(j,1,T)
    {
        if (i+j<=n) add(num[i][j],num[i+j][j],1);
        if (i-j>0) add(num[i][j],num[i-j][j],1);
        add(num[i][j],i,0);
    }
    F(i,1,m)
    {
        if (p[i]<=T) add(b[i],num[b[i]][p[i]],0);
        else
        {
            for (int j=1;b[i]+j*p[i]<=n;j++) add(b[i],b[i]+j*p[i],j);
            for (int j=1;b[i]-j*p[i]>=1;j++) add(b[i],b[i]-j*p[i],j);
        }
    }
    SPFA();
    return 0;
}