URAL 1029 Ministry

题目链接

题目大意

有一个n层楼的办公大楼，每一层楼有m个办公室，你希望通过一份文件，文件被通过当且仅当最顶层楼的人签过字。每个人签字有一定的花费，并且他签字只会当如下情况之一满足时。

• 他正楼下的人签过字。
• 他左面或右面的人签过字。

求这份文件被通过的最小花费。

思路

考虑DP，状态显然：$f[i][j]$表示签到第$i$层第$j$个办公室的最小花费。
转移时枚举三种方式转移。
有环，怎么办呢，考虑Dijkstra转移。或者想一下我当时的脑洞。

有费用啊，是不是可以费用流啊，数据范围挺好的
是可以费用流，费用流流量都是1。。。。。
那不就是个最短路吗

解法：连边跑Dijkstra即可，复杂度可以承受

另一种思路

考虑刚才的DP，有环不是嘛吗，但是发现，一个点要么是从左边转移来，在这里终止或往右边推，要么是从右边过来，在这里终止或者往左边推。
其实是不可能有说“从左来兜个圈子再回左边去”，这很明显是浪费钱的。
所以说，从左和从右各推一次即可
复杂度比最短路优秀一个log

代码（最短路）

#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<ctime>
#include<map>
#include<queue>
#include<vector>
#include<stack>
#include<set>
#include<cctype>
#define pa pair<int,int>
#define INF 0x3f3f3f3f
#define inf 0x3f
#define fi first
#define se second
#define mp make_pair
#define ll long long
#define ull unsigned long long
#define pb push_back

using namespace std;

inline ll read()
{
    long long f=1,sum=0;
    char c=getchar();
    while (!isdigit(c)) {if (c=='-') f=-1;c=getchar();}
    while (isdigit(c)) {sum=sum*10+c-'0';c=getchar();}
    return sum*f;
}
const int MAXN=50030;
const int MAXM=MAXN*10;
struct edge{
    int next,to,val,from;
}e[MAXM];
int head[MAXN],cnt;
void addedge(int u,int v,int w)
{
    e[++cnt].next=head[u];
    e[cnt].from=u;
    e[cnt].to=v;
    e[cnt].val=w;
    head[u]=cnt;
} 
queue <int> q;
ll dis[MAXN];
int from[MAXN],a[510][510];
bool inqueue[MAXN];
int n,m,s,t;
int id[510][510];
void spfa()
{
    for (int i=0;i<=n*m+10;i++)
        dis[i]=1e18;
    dis[s]=0,q.push(s),inqueue[s]=1;
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for (int i=head[x];i;i=e[i].next)
        {
            int v=e[i].to;
            if (dis[v]>dis[x]+e[i].val)
            {
                dis[v]=dis[x]+e[i].val;
                from[v]=i;
                if (!inqueue[v])
                    q.push(v),inqueue[v]=1;
            }
        }
        inqueue[x]=0;
    } 
}
vector <int> lu;
int main()
{
    scanf("%d%d",&n,&m);
    int tot=0;
    for (int i=0;i<n;i++)
        for (int j=0;j<m;j++)
            scanf("%d",&a[i][j]),id[i][j]=tot++;
    s=n*m+1,t=s+1;
    for (int i=0;i<m;i++)
        addedge(s,id[0][i],a[0][i]);
    for (int i=0;i<n;i++)
        for (int j=0;j<m;j++)
        {
            if (j) addedge(id[i][j],id[i][j-1],a[i][j-1]);
            if (j<m-1) addedge(id[i][j],id[i][j+1],a[i][j+1]);
            if (i<n-1) addedge(id[i][j],id[i+1][j],a[i+1][j]); 
        }
    for (int i=0;i<m;i++)
        addedge(id[n-1][i],t,0);
    spfa();
    for (int i=t;e[from[i]].from!=s;i=e[from[i]].from)
        lu.push_back(e[from[i]].from%m+1);
    for (int i=lu.size()-1;i>=0;i--)
        printf("%d ",lu[i]);
    return 0;
}
/*
3 4
10 10 1 10
2 2 2 10
1 10 10 10
*/

代码（普通DP）

#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<ctime>
#include<map>
#include<queue>
#include<vector>
#include<stack>
#include<set>
#include<cctype>
#define pa pair<int,int>
#define INF 0x3f3f3f3f
#define inf 0x3f
#define fi first
#define se second
#define mp make_pair
#define ll long long
#define ull unsigned long long
#define pb push_back

using namespace std;

inline ll read()
{
    long long f=1,sum=0;
    char c=getchar();
    while (!isdigit(c)) {if (c=='-') f=-1;c=getchar();}
    while (isdigit(c)) {sum=sum*10+c-'0';c=getchar();}
    return sum*f;
}
const int MAXN=510;
int a[MAXN][MAXN],f[MAXN][MAXN];
int from[MAXN][MAXN];
vector <int> lu;
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            scanf("%d",&a[i][j]);
    memset(f,inf,sizeof(f));
    memset(f[0],0,sizeof(f[0]));
    for (int i=1;i<=n;i++)
    {
        for (int j=1;j<=m;j++)
        {
            if (f[i-1][j]+a[i][j]<f[i][j])
                f[i][j]=f[i-1][j]+a[i][j],from[i][j]=j;
            if (f[i][j-1]+a[i][j]<f[i][j])
                f[i][j]=f[i][j-1]+a[i][j],from[i][j]=j-1;
        }
        for (int j=m;j>=1;j--)
        {
            if (f[i-1][j]+a[i][j]<f[i][j])
                f[i][j]=f[i-1][j]+a[i][j],from[i][j]=j;
            if (f[i][j+1]+a[i][j]<f[i][j])
                f[i][j]=f[i][j+1]+a[i][j],from[i][j]=j+1;
        }
    }
    int ans=INF,pos;
    for (int i=1;i<=m;i++)
        if (f[n][i]<ans)
            ans=f[n][i],pos=i;
    int h=n,l=pos;
    while (h!=0)
    {
        lu.push_back(l);
        int tmp=from[h][l];
        if (tmp==l) h--;
        else l=tmp;
    }
    for (int i=lu.size()-1;i>=0;i--)
        printf("%d ",lu[i]);
    return 0;
}