The Glorious Karlutka River =)

sgu438:http://acm.sgu.ru/problem.php?contest=0&problem=438

题意：有一条东西向流淌的河,宽为 W,河中有 N 块石头,每块石头的坐标(Xi, Yi)和最大承受人数 Ci 已知。现在有 M 个游客在河的南岸,他们想穿越这条河流,但是每个人每次最远只能跳 D 米,每跳一次耗时 1 秒。问他们能否全部穿越这条河流,如果能,最少需要多长时间。 <= N <= 50, 0 < M <= 50, 0 <= D <= 1000, 0 < W(0<= 1000, 0 < Xi < 1000, 0 < Yi < W, 0 <= Ci <= 1000)。刚看完这题，想当然的认为它是一道最小费用流问题。但是当WA之后我才明白，这题并不是去求一个给定网络的最大流，而是计算这个网络随着时间推移每次能够留出多少流量。我们通过枚举时间的方式来决定在什么时刻能够把所有的人全部送到对岸。注意人是可以从河这岸的任意x坐标出发的。

题解：这是一道费用流。具体的看代码吧。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 #include<cmath>
 #define inf 100000000
 using namespace std;
 const int E=;
 const int N=;
 struct Node{
    int v, cap, cost, next;    //  re记录逆边的下标。
 }edge[E];
 int n, m;
 int  ans;
 int k, head[N];
 int que[N], pre[N], dis[N];
 bool vis[N];
 void init(){//初始化
    k=ans=;
    memset(head,-,sizeof(head));
 }
 void addEdge(int u, int v, int ca, int co){
     edge[k].v = v;
     edge[k].cap = ca;
     edge[k].cost = co;
     edge[k].next = head[u];
     head[u] = k ++;
     edge[k].v = u;
     edge[k].cap = ;
     edge[k].cost = -co;
     edge[k].next = head[v];
     head[v] = k ++;
 }
 bool spfa(){                  //  源点为0，汇点为n。
     int i;
     for(i = ; i <=*m+;i++){
         dis[i] = inf;
         vis[i] = false;
     }
     queue<int>Q;
     Q.push();
     dis[]=;
     vis[] = true;
     while(!Q.empty()){       //  这里最好用队列，有广搜的意思，堆栈像深搜。
         int u = Q.front();
         Q.pop();
         vis[u]=;
         for(i = head[u]; i != -; i = edge[i].next){
             int v = edge[i].v;
             if(edge[i].cap && dis[v] > dis[u] + edge[i].cost){
                 dis[v] = dis[u] + edge[i].cost;
                 pre[v] = i;
                 if(!vis[v]){
                     vis[v] = true;
                    Q.push(v);
                 }
             }
         }
         vis[u] = false;
     }
     if(dis[*m+] == inf) return false;
     return true;
 }
 int  end(){
     int u, p, sum = inf;
     for(u = *m+; u != ; u = edge[p^].v){//0是超级源点
         p = pre[u];
         sum = min(sum, edge[p].cap);
     }
     for(u = *m+; u != ; u = edge[p^].v){
         p = pre[u];
         edge[p].cap -= sum;
         edge[p^].cap += sum;
     }
     return sum;
 }
 double d,w;
 struct  Point{
    double x;
    double y;
    int val;
 }num[];
 int main(){
    while(~scanf("%d%d%lf%lf",&m,&n,&d,&w)){
         init();//初始化
        for(int i=;i<=m;i++){
         scanf("%lf%lf%d",&num[i].x,&num[i].y,&num[i].val);
        }
        for(int i=;i<=m;i++){
           if(abs(num[i].y)<=d){
               addEdge(,i,n,);
           }
        }
        for(int i=;i<=m;i++){
           for(int j=;j<=m;j++){
              if(i==j)continue;
              double diss=sqrt((num[i].x-num[j].x)*(num[i].x-num[j].x)+(num[i].y-num[j].y)*(num[i].y-num[j].y));
              if(diss<=d){
                 addEdge(i+m,j,inf,);
              }
           }
           if(abs(num[i].y-w)<=d){
             addEdge(i+m,*m+,inf,);
           }
           addEdge(i,i+m,num[i].val,);
        }
        if(w<=d)addEdge(,*m+,inf,);
        int  s=n,now=,sum=;
        ans=inf;
        while(spfa()){
           int y=end();
            s-=(dis[*m+]-now)*sum+y;
            if(s<)s=;
            sum+=y;now=dis[*m+];
            int temp=now+(int)ceil(s*1.0/sum);
            if(temp<ans)ans=temp;
        }
        if(ans==inf)printf("IMPOSSIBLE\n");
        else
        printf("%d\n",ans);
    }
 }