Luogu2662 牛场围栏（最短路）

　　小凯的疑惑升级版的升级版。答案若存在不会超过3000²-3000，暴力dp似乎勉强可以过。当然这不优美。

　　注意到如果能拼出长度为l的围栏，就一定能拼出长度为l+kx的围栏，其中x为最短的（或任意一个）围栏长度。这样将值域范围缩小到了3000以内。于是将同余类间连长为木料长度的边，求出0为源点到每个点的最短路（据说spfa有奇效），答案即为max{dis}-x。

　　这个东西就叫同余最短路。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 110
#define M 3010
int n,m,a[N],q[M],d[M],f[M][M],p[M],t=,P=,ans;
bool flag[M];
struct data{int to,nxt,len;
}edge[M*M];
void addedge(int x,int y,int z){t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].len=z,p[x]=t;}
int inc(int &x){x++;if (x>P+) x-=P+;return x;}
void spfa()
{
    memset(flag,,sizeof(flag));
    memset(d,,sizeof(d));d[]=;
    int head=,tail=;q[]=;
    do
    {
        int x=q[inc(head)];flag[x]=;
        for (int i=p[x];i;i=edge[i].nxt)
        if (d[x]+edge[i].len<d[edge[i].to])
        {
            d[edge[i].to]=d[x]+edge[i].len;
            if (!flag[edge[i].to])
            {
                q[inc(tail)]=edge[i].to;
                flag[edge[i].to]=;
            }
        }
    }while (head!=tail);
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("cow.in","r",stdin);
    freopen("cow.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    for (int i=;i<=n;i++)
    {
        a[i]=read(),P=min(P,max(a[i]-m,));
        for (int j=max(,a[i]-m);j<=a[i];j++)
        flag[j]=;
    }
    for (int i=;i<P;i++)
        for (int j=;j<=;j++)
        if (flag[j]) addedge(i,(i+j)%P,j);
    spfa();
    for (int i=;i<P;i++)
    ans=max(ans,d[i]-P);
    cout<<(ans==?-:ans);
    return ;
}