【ZJOI2013】bzoj3112 防守战线

首先可以列出形如

$$\sum_{j=l_i}^{r_i} X_j\ge d_i\ \ \ \ (1\le i \le m)$$
$$\text{Minimize} \sum_{i=1}^n X_ic_i$$
的线性规划式子。但是注意到 $X_i$的分布很杂乱，不好处理。考虑转化为对偶问题。
$$\sum_{l_j\le i \le r_j} Y_j\le c_i\ \ \ \ (1\le i\le n)$$
$$\text{Maximize} \sum_{i=1}^mY_id_i$$
这样就可以差分建图跑最大费用流了。但是这样是会T的，spfa加了SLF以后在bzoj上卡过了。

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=1020,maxe=200010,maxm=10010,p=1011,s=1005,t=1006,oo=0x3f3f3f3f;
int fir[maxn],ne[maxe],to[maxe],w[maxe],c[maxe],
que[maxn],in[maxn],dis[maxn],mnw[maxn],pre[maxn],
a[maxn],l[maxm],r[maxm],d[maxm],
n,m,num;
int rd()
{
    int x=0;
    char c=getchar();
    while (c<'0'||c>'9') c=getchar();
    while (c>='0'&&c<='9')
    {
        x=x*10+c-'0';
        c=getchar();
    }
    return x;
}
void check()
{
    for (int i=0;i<=t;i++)
        for (int j=fir[i];j;j=ne[j])
            if (w[j])
                printf("%d->%d:%d,%d\n",i,to[j],w[j],c[j]);
}
void add(int u,int v,int x,int y)
{
    num++;
    ne[num*2]=fir[u];
    fir[u]=num*2;
    to[num*2]=v;
    w[num*2]=x;
    c[num*2]=y;
    ne[num*2+1]=fir[v];
    fir[v]=num*2+1;
    to[num*2+1]=u;
    w[num*2+1]=0;
    c[num*2+1]=-y;
}
void inc(int &x)
{
    x=x+1==p?0:x+1;
}
int find(int &ret)
{
    int hd=0,tl=1,u,v;
    que[0]=s;
    for (int i=0;i<=n;i++) dis[i]=-oo;
    dis[t]=-oo;
    dis[s]=0;
    mnw[s]=oo;
    in[s]=1;
    while (hd!=tl)
    {
        u=que[hd];
        inc(hd);
        for (int i=fir[u];i;i=ne[i])
            if (w[i]&&dis[v=to[i]]<dis[u]+c[i])
            {
                dis[v]=dis[u]+c[i];
                pre[v]=i;
                mnw[v]=min(mnw[u],w[i]);
                if (!in[v])
                {
                    in[v]=1;
                    que[tl]=v;
                    if (dis[que[tl]]>dis[que[hd]]) swap(que[tl],que[hd]); 
                    inc(tl);
                }
            }
        in[u]=0;
    }
    if (dis[t]==-oo) return 0;
    for (int i=pre[t];i;i=pre[to[i^1]])
    {
        w[i]-=mnw[t];
        w[i^1]+=mnw[t];
    }
    ret=dis[t]*mnw[t];
    return 1;
}
int main()
{
    int ans=0,x;
    n=rd();
    m=rd();
    for (int i=1;i<=n;i++)
    {
        a[i]=rd();
        if (a[i]>a[i-1]) add(i-1,t,a[i]-a[i-1],0);
        else add(s,i-1,a[i-1]-a[i],0);
        add(i,i-1,oo,0);
    }
    add(s,n,a[n],0);
    for (int i=1;i<=m;i++)
    {
        l[i]=rd();
        r[i]=rd();
        d[i]=rd();
        add(r[i],l[i]-1,oo,d[i]);
    }
    //check();
    while (find(x)) ans+=x;
    printf("%d\n",ans);
}