【BZOJ】【1520】【POI2006】Szk-Schools

网络流/费用流

---

　　比较裸的一道题

　　依旧是二分图模型，由源点S连向每个学校 i (1,0)，「注意是连向第 i 所学校，不是连向学校的标号m[i]……唉这里WA了一次」

　　然后对于每所学校 i 连接 j+n $(a[i]\leq j \leq b[i])$ 流量为1，费用为 $abs(m[i]-j)*k[i]$ ，最后对于每个标号 j 连边 j+n->T 流量为1费用为0。

　　跑完费用流以后看流量是否为n，如果不是就说明无解。

/**************************************************************
    Problem: 1520
    User: Tunix
    Language: C++
    Result: Accepted
    Time:1452 ms
    Memory:5964 kb
****************************************************************/
 
//BZOJ 1520
#include<vector>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
inline int getint(){
    int v=0,sign=1; char ch=getchar();
    while(ch<'0'||ch>'9'){ if (ch=='-') sign=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){ v=v*10+ch-'0'; ch=getchar();}
    return v*sign;
}
const int N=410,M=200000,INF=~0u>>2;
typedef long long LL;
/******************tamplate*********************/
int n,ans,flow;
struct edge{int from,to,v,c;};
struct Net{
    edge E[M];
    int head[N],next[M],cnt;
    void ins(int x,int y,int z,int c){
        E[++cnt]=(edge){x,y,z,c};
        next[cnt]=head[x]; head[x]=cnt;
    }
    void add(int x,int y,int z,int c){
        ins(x,y,z,c); ins(y,x,0,-c);
    }
    int from[N],Q[M],d[N],S,T,ed;
    bool inq[N],sign;
    bool spfa(){
        int l=0,r=-1;
        F(i,1,T) d[i]=INF;
        d[S]=0; Q[++r]=S; inq[S]=1;
        while(l<=r){
            int x=Q[l++];
            inq[x]=0;
            for(int i=head[x];i;i=next[i])
                if(E[i].v>0 && d[x]+E[i].c<d[E[i].to]){
                    d[E[i].to]=d[x]+E[i].c;
                    from[E[i].to]=i;
                    if (!inq[E[i].to]){
                        Q[++r]=E[i].to;
                        inq[E[i].to]=1;
                    }
                }
        }
        return d[T]!=INF;
    }
    void mcf(){
        int x=INF;
        for(int i=from[T];i;i=from[E[i].from])
            x=min(x,E[i].v);
        for(int i=from[T];i;i=from[E[i].from]){
            E[i].v-=x;
            E[i^1].v+=x;
        }
        flow+=x;
        ans+=x*d[T];
    }
    void init(){
        n=getint(); cnt=1;
        S=0; T=n*2+1;
        int m,x,y,k;
        F(i,1,n){
            m=getint(); x=getint(); y=getint(); k=getint();
            add(S,i,1,0);
            F(j,x,y) add(i,j+n,1,k*abs(j-m));
            add(i+n,T,1,0);
        }
        while(spfa()) mcf();
        if (flow==n) printf("%d\n",ans);
        else printf("NIE");
    }
}G1;
 
int main(){
#ifndef ONLINE_JUDGE
    freopen("1520.in","r",stdin);
    freopen("1520.out","w",stdout);
#endif
    G1.init();
    return 0;
}

 

1520: [POI2006]Szk-Schools

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 428  Solved: 220
[Submit][Status][Discuss]

Description

Input

Output

如果有可行解, 输出最小代价,否则输出NIE.

Sample Input

5
1 1 2 3
1 1 5 1
3 2 5 5
4 1 5 10
3 3 3 1

Sample Output

9

HINT

Source

[Submit][Status][Discuss]