P1006传纸条

菜狗不会dp 网络流复健

#include <bits/stdc++.h>
/*
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
*/
using namespace std;

const double eps = 1e-10;
const double pi = 3.1415926535897932384626433832795;
const double eln = 2.718281828459045235360287471352;

#define f(i, a, b) for (int i = a; i <= b; i++)
#define scan(x) scanf("%d", &x)
#define mp make_pair
#define pb push_back
#define lowbit(x) (x&(-x))

#define fi first
#define se second
#define SZ(x) int((x).size())
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define summ(a) (accumulate(all(a), 0ll))

typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef vector<int> vi;

using ll=long long;

int n,m,a[100][100],tot=1,hd[(int)7e3+9],cur[(int)7e3+9],dep[(int)7e3+9],res;
bitset<(int)7e3+9> vis;
struct Edge{
    int to,nxt,v,c;
}e[(int)1e5+9];
void add(int u,int v,int w,int c){
    e[++tot]={v,hd[u],w,c},hd[u]=tot;
    e[++tot]={u,hd[v],0,-c},hd[v]=tot;
}

bool bfs(int st,int en){
    memset(dep,0,sizeof(dep));
    memcpy(cur,hd,sizeof(hd));
    queue<int>q;
    q.push(st);
    dep[st]=1,vis[st]=1;
    while(!q.empty()){
        auto u=q.front();
        q.pop(),vis[u]=0;
        for(int eg=hd[u];eg;eg=e[eg].nxt){
            if(dep[e[eg].to]<dep[u]+e[eg].c&&e[eg].v){
                dep[e[eg].to]=dep[u]+e[eg].c;
                if(!vis[e[eg].to])q.push(e[eg].to),vis[e[eg].to]=1;
            }
        }
    }
    return dep[en]!=0;
}

int dfs(int u,int en,int flow){
    if(u==en)return flow;
    vis[u]=1;
    int r=flow;
    for(int eg=cur[u];eg&&r;eg=e[eg].nxt){
        cur[u]=eg;
        if(e[eg].v&&dep[e[eg].to]==dep[u]+e[eg].c&&!vis[e[eg].to]){
            int c=dfs(e[eg].to,en,min(r,e[eg].v));
            if(c)res+=c*e[eg].c;
            r-=c;
            e[eg].v-=c,e[eg^1].v+=c;
        }
    }
    vis[u]=0;
    return flow-r;
}

int dinic(int st,int en){
    int ret=0;
    while(bfs(st,en)){
        ret+=dfs(st,en,0x3f3f3f3f);
    }
    return ret;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cin>>n>>m;
    f(i,1,n)f(j,1,m)cin>>a[i][j];
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j){
            if(i!=n)add(i*m+j-m+n*m,i*m+j,1,a[i+1][j]);
            if(j!=m)add(i*m-m+j+n*m,i*m-m+j+1,1,a[i][j+1]);
            if((i==1&&j==1)||(i==n&&j==m)){
                add(i*m-m+j,i*m-m+j+n*m,2,0);
            }else add(i*m-m+j,i*m-m+j+n*m,1,0);
        }
    }
    dinic(1,2*n*m);
    cout<<res;
    return 0;
}