2696: 航班安排 费用流

先让我重新说一下题意qwq。
有$n$个机场，机场两两之间需要的最短时间是$f_{i,j}$，花费的费用是一个固定值$c_{i,j}$，有$m$个请求，在$s$时刻从$a$机场起飞，在恰好$t$时刻到达$b$机场，可以净获利$cc$，问最大获利。
也就是说对于路线$i->j$，你可以为了恰好到达而多飞一会儿。
数据保证给出的请求是合法的，即一定能从$i$到达$j$。

然后我们就可以愉快的建图跑费用流了。
把每个请求拆为两个点$x_1,x_2$，连边$(x_1,x_2,1,cc_i)$，然后对于可达的两个请求$x,y$，连边$(x_2,y_1,1,-c_{b_x,a_y})$，然后和源点汇点连边，即$(S,x_1,1,-c_{1,a_x}),(x_2,T,1,-c_{b_x,1})$，最后汇点和超级汇点连边限制流量 k <script type="math/tex" id="MathJax-Element-1846">k</script>。
可达性根据距离判断，由于数据保证两点之间的距离一定是最短路，所以直接判就行了。

#include<iostream>
#include<cstdio>
#include<cstring>
#define inf 1000000007
#define N 405
#define M 100005
using namespace std;
int n,m,k,t,A,T,cnt=1,ans;
int head[N],dis[N],q[N],path[N];
int list[M],next[M],from[M],flow[M],cost[M];
bool vis[N];
int a[N],b[N],st[N],en[N],cc[N];
int f[N][N],c[N][N];
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void insert(int x,int y,int z,int w)
{
    next[++cnt]=head[x];
    head[x]=cnt;
    from[cnt]=x;
    list[cnt]=y;
    flow[cnt]=z;
    cost[cnt]=w;
}
inline bool spfa()
{
    for (int i=0;i<=T;i++) dis[i]=-inf;
    int t=0,w=1,x;
    q[1]=0; dis[0]=0;
    while (t!=w)
    {
        t=(t+1)%N;
        x=q[t];
        for (int i=head[x];i;i=next[i])
            if (flow[i]&&dis[list[i]]<dis[x]+cost[i])
            {
                dis[list[i]]=dis[x]+cost[i];
                path[list[i]]=i;
                if (!vis[list[i]])
                {
                    vis[list[i]]=1;
                    w=(w+1)%N;
                    q[w]=list[i];
                }
            }
        vis[x]=0;
    }
    return dis[T]!=-inf;
}
inline void mcf()
{
    int x=inf;
    for (int i=path[T];i;i=path[from[i]]) x=min(x,flow[i]);
    for (int i=path[T];i;i=path[from[i]])
        flow[i]-=x,flow[i^1]+=x,ans+=x*cost[i];
}
int main()
{
    n=read(); m=read(); k=read(); t=read(); A=m<<1|1; T=A+1;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++) 
            f[i][j]=read();
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            c[i][j]=read();
    for (int i=1;i<=m;i++)
        a[i]=read()+1,b[i]=read()+1,st[i]=read(),en[i]=read(),cc[i]=read();
    for (int i=1;i<=m;i++)
        if (f[1][a[i]]<=st[i]) 
            insert(0,i,1,-c[1][a[i]]),insert(i,0,0,c[1][a[i]]);
    for (int i=1;i<=m;i++)
        for (int j=1;j<=m;j++)
            if (en[i]+f[b[i]][a[j]]<=st[j]) 
                insert(i+m,j,1,-c[b[i]][a[j]]),insert(j,i+m,0,c[b[i]][a[j]]);
    for (int i=1;i<=m;i++)
        insert(i,i+m,1,cc[i]),insert(i+m,i,0,-cc[i]);
    for (int i=1;i<=m;i++)
        if (en[i]+f[b[i]][1]<=t) 
            insert(i+m,A,1,-c[b[i]][1]),insert(A,i+m,0,c[b[i]][1]);
    insert(A,T,k,0); insert(T,A,0,0);
    while (spfa()) mcf();
    cout << ans << endl;
    return 0;
}