BZOJ1070[SCOI2007]修车

AC通道：http://www.lydsy.com/JudgeOnline/problem.php?id=1070

codevs上也有哦：http://codevs.cn/problem/2436/

这题感觉还是比较好想的。

因为算的是平均等待时间，那么就要考虑分配问题和顺序问题。

首先我们从s向每个汽车连容量为1，费用为0的边。

然后我们将师傅拆点，拆成n分，可以看做是倒数多少个的状态。

那么我们将每个汽车(i)向每个状态(j,k)都连一条容量为1,费用为k*w[j][i]的边。

表示师傅j在倒数第k个做i产生的贡献值。

然后再将每个状态往t连容量为1，费用为0的边。

跑最小费用最大流就可以了。

这题，因为连边多一点，所以zkw稍慢一点点：

#include<cstdio>
#include<cstring>
#include<algorithm>

using namespace std;

const int maxn=65*65;
const int INF=0x3f3f3f3f;

struct Node{
    int data,next,low,cost;
}node[maxn*9*2];

#define www node[point].low
#define now node[point].data
#define ccc node[point].cost
#define then node[point].next

int n,m,cnt;
int s,t,ans;
int head[maxn],cur[maxn];
int dis[maxn];
int a[65][65];
bool vis[maxn];

void add(int u,int v,int w,int c){
    node[cnt].data=v;node[cnt].next=head[u];node[cnt].low=w;node[cnt].cost=c;head[u]=cnt++;
    node[cnt].data=u;node[cnt].next=head[v];node[cnt].low=0;node[cnt].cost=-c;head[v]=cnt++;
}

int dfs(int x,int low){
    if(x==t) return low;
    int Low;vis[x]=true;
    for(int &point=cur[x];point!=-1;point=then)
        if(www && !vis[now] && dis[x]==dis[now]+ccc){
            Low=dfs(now,min(low,www));
            if(Low){
                www-=Low,node[point^1].low+=Low;
                return Low;
            }
        }
    return 0;
}

bool modify(){
    int ad=INF;
    for(int i=s;i<=t;i++)
    if(vis[i])
        for(int point=head[i];point!=-1;point=then)
        if(www && !vis[now])
            ad=min(ad,dis[now]+ccc-dis[i]);
    if(ad==INF) return false;
    for(int i=s;i<=t;i++)
    if(vis[i])
        vis[i]=false,dis[i]+=ad;
    return true;
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("1070.in","r",stdin);
    freopen("1070.out","w",stdout);
#endif

    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&a[i][j]);
    t=n*(m+1)+1;
    for(int i=s;i<=t;i++) head[i]=-1;
    for(int i=1;i<=n;i++) add(s,i,1,0);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            for(int k=1;k<=n;k++)
                add(i,j*n+k,1,k*a[i][j]);
    for(int i=n+1;i<t;i++) add(i,t,1,0);
    
    int flag;
    do{
        for(int i=s;i<=t;i++) cur[i]=head[i];
        while(flag=dfs(s,INF))
            ans+=flag*dis[s];
    }while(modify());
    printf("%.2lf",(double)ans/n);
    return 0;
}

spfa似乎快了那么40ms？never mind，反正只是当练习模板了。

#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>

using namespace std;

const int maxn=65*65;
const int INF=0x3f3f3f3f;

struct Node{
    int data,next,low,cost;
}node[maxn*9*2];

#define www node[point].low
#define now node[point].data
#define ccc node[point].cost
#define then node[point].next

int n,m,cnt;
int s,t,ans;
int head[maxn],cur[maxn];
int dis[maxn],pre[maxn];
int a[65][65];
bool inq[maxn];
deque<int >q;

void add(int u,int v,int w,int c){
    node[cnt].data=v;node[cnt].next=head[u];node[cnt].low=w;node[cnt].cost=c;head[u]=cnt++;
    node[cnt].data=u;node[cnt].next=head[v];node[cnt].low=0;node[cnt].cost=-c;head[v]=cnt++;
}

bool SPFA(){
    memset(dis,0x3f,sizeof(dis));
    int x,Low=INF;
    q.push_back(s);dis[s]=0;pre[s]=pre[t]=-1;
    while(!q.empty()){
        x=q.front();q.pop_front();inq[x]=false;
        for(int point=head[x];point!=-1;point=then)
            if(www && dis[now]>dis[x]+ccc){
                dis[now]=dis[x]+ccc;pre[now]=point^1;
                if(!inq[now]){
                    if(q.empty() || dis[q.front()]>dis[now])
                        q.push_front(now);
                    else
                        q.push_back(now);
                    inq[now]=true;
                }
            }
    }
    if(pre[t]<0) return false;
    for(int point=pre[t];point!=-1;point=pre[now])
        Low=min(Low,node[point^1].low);
    for(int point=pre[t];point!=-1;point=pre[now])
        node[point^1].low-=Low,www+=Low;
    ans+=dis[t]*Low;
    return true;
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("1070.in","r",stdin);
    freopen("1070.out","w",stdout);
#endif

    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&a[i][j]);
    t=n*(m+1)+1;
    for(int i=s;i<=t;i++) head[i]=-1;
    for(int i=1;i<=n;i++) add(s,i,1,0);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            for(int k=1;k<=n;k++)
                add(i,j*n+k,1,k*a[i][j]);
    for(int i=n+1;i<t;i++) add(i,t,1,0);
    while(SPFA());
    printf("%.2lf",(double)ans/n);
    return 0;
}