洛谷P4009汽车加油行驶问题——网络流24题(最短路)

题目：https://www.luogu.org/problemnew/show/P4009

网络流24题中不是网络流的最短路题；

把每个点拆成各个油量上的点，根据要求连边即可；

注意：点数最大为100*100*11，因为虽然k<=10，但还有k=0的状态！(竟然因为边界调了一晚上)

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
queue<int>q;
int const MAXN=,inf=1e9;//MAXN不可为100005
int n,k,a,b,c,nd[][][],cnt,ct=,head[MAXN];
int dis[MAXN];
bool vis[MAXN];
struct N{
    int to,next,w;
    N(int t=,int n=,int w=):to(t),next(n),w(w) {}
}edge[MAXN*];
void spfa()
{
    memset(dis,,sizeof dis);
    int s=nd[k][][];
    dis[s]=;vis[s]=;q.push(s);
    while(q.size())
    {
        int x=q.front();q.pop();vis[x]=;
//        cout<<x<<endl;
        for(int i=head[x];i;i=edge[i].next)
        {
            int u=edge[i].to;
            if(dis[u]>dis[x]+edge[i].w)
            {
//                cout<<dis[x]<<endl;
                dis[u]=dis[x]+edge[i].w;
                if(!vis[u])vis[u]=,q.push(u);
            }
        }
//        cout<<q.size()<<endl;
    }
//    cout<<"spfa"<<endl;
//    ll mn=inf;
//    for(int l=0;l<=k;l++)
//        mn=min(mn,dis[nd[l][n][n]]);
    printf("%d",dis[cnt+]);
}
int main()
{
    scanf("%d%d%d%d%d",&n,&k,&a,&b,&c);//a加油 b返回 c建库
    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
            for(int l=;l<=k;l++)
                nd[l][i][j]=++cnt;
    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
        {
            bool d;
            scanf("%d",&d);
            if(!d)
            {
                edge[++ct]=N(nd[k][i][j],head[nd[][i][j]],a+c),head[nd[][i][j]]=ct;
                for(int l=;l<=k;l++)
                {
                    int x=nd[l][i][j];
                    if(i>)edge[++ct]=N(nd[l-][i-][j],head[x],b),head[x]=ct;
                    if(j>)edge[++ct]=N(nd[l-][i][j-],head[x],b),head[x]=ct;
                    if(i<n)edge[++ct]=N(nd[l-][i+][j],head[x],),head[x]=ct;
                    if(j<n)edge[++ct]=N(nd[l-][i][j+],head[x],),head[x]=ct;
//                    if(l<k)edge[++ct]=N(nd[k][i][j],head[x],a+c),head[x]=ct;//不优！
                }
            }
            if(d)
            {
                int x=nd[k][i][j];
                for(int l=;l<k;l++)
                    edge[++ct]=N(x,head[nd[l][i][j]],a),head[nd[l][i][j]]=ct;
                if(i>)edge[++ct]=N(nd[k-][i-][j],head[x],b),head[x]=ct;
                if(j>)edge[++ct]=N(nd[k-][i][j-],head[x],b),head[x]=ct;
                if(i<n)edge[++ct]=N(nd[k-][i+][j],head[x],),head[x]=ct;
                if(j<n)edge[++ct]=N(nd[k-][i][j+],head[x],),head[x]=ct;
            }
        }
    for(int l=;l<=k;l++)
        edge[++ct]=N(cnt+,head[nd[l][n][n]],),head[nd[l][n][n]]=ct;
//    cout<<cnt<<endl;
    spfa();
    return ;
}