438. The Glorious Karlutka River =)

题目链接：http://acm.sgu.ru/problem.php?contest=0&problem=438

题目大意：河中有一些漂浮物，每个漂浮物有一个容量。求m个人过河的最短时间。

题目思路：经典题目,加上时间限定条件,并且容量在点上,需要拆点,假设答案时间是T,那么每个点都要拆成2*T个点具体建法:S->SS为人数,SS->V无穷,UT,1->UT,2为点的容量,UT-1,2->VT,1无穷,UT-1,2->UT,1无穷,UT-1,2->T无穷效率关键是看怎么建图个方案和答案T的枚举,我用了四个方法,时间效率比较如下:sap(邻接表)按时间从小到大遍历,每次都重新建图重新增广2173MSsap(邻接表)二分时间,每次都重新建图重新增广409MSEK(邻接表)按时间从小到大,但是每次增加2*n个点后在原来的基础上增广47MSsap(邻接表)按时间从小到大,但是每次增加2*n个点后在原来的基础上增广(但是dis和gap初始化)772MS(这里sap比EK还慢,也许是因为每次推的时候dis和gap清空了(不初始化的话答案不出来),所以推的很慢,不知道有没有更好的算法)（以上思路摘自网上）

开始的时候我想到的就是拆点加二分，不敢确定，还是看了题解，看到题解有这总方法，试了一下，一直tle,于是还是只能用逐步加点的方式了，用dinic的模板过了。

#include<stdio.h>#include<stdlib.h>#include<string.h>#include<string>#include<queue>#include<algorithm>#include<vector>#include<stack>#include<list>#include<iostream>#include<map>#include<math.h>using namespace std;#define inf 0x3f3f3f3f#define Max 55555int max(int a,int b){    return a>b?a:b;}int min(int a,int b){    return a<b?a:b;}int x[100],y[100],c[100],m,n,d,w;int dis[Max],p[Max];int src,dest,eid,ans;int q[Max];struct node {     int to,c,next; }e[2222222];void addedge(int u,int v,int c){    e[eid]=(node){v,c,p[u]};    p[u]=eid++;    e[eid]=(node){u,0,p[v]};    p[v]=eid++;}int bfs(){int i,u,v,head=0,tail=0;for(i=0;i<=dest+1;i++) dis[i]=0;dis[src]=1;q[tail++]=src;while(head<tail){u=q[head++];for(i=p[u];~i;i=e[i].next)if(e[i].c&&dis[v=e[i].to]==0){dis[v]=dis[u]+1;if(v==dest)return 1;q[tail++]=v;}}return 0;}int dfs(int u,int limit){if(u==dest)return limit;int v,tmp,cost=0;for(int i=p[u];~i;i=e[i].next)if(e[i].c&&dis[u]+1==dis[v=e[i].to]){tmp=dfs(v,min(limit-cost,e[i].c));if(tmp){e[i].c-=tmp;e[i^1].c+=tmp;cost+=tmp;}elsedis[v]=-1;}return cost;}int Dinic(){while(bfs())        ans+=dfs(src,inf);return ans;}bool find(int u){    int i;    dis[u]=1;    if(w-y[u]<=d)    {        dis[n+1]=1;        return true;    }    for(i=1;i<=n;i++)    {        if(!dis[i]&&((x[u]-x[i])*(x[u]-x[i])+(y[u]-y[i])*(y[u]-y[i]))<=d*d)        {           if(find(i))                return true;        }    }    return false;}int main(){    int i,j,k;    while(scanf("%d%d%d%d",&n,&m,&d,&w)!=EOF)    {        int cnt=0;        for(i=1;i<=n;i++)        {            scanf("%d%d%d",&x[i],&y[i],&c[i]);            if(c[i]>0)            {                x[++cnt]=x[i];                y[cnt]=y[i];                c[cnt]=c[i];            }        }        n=cnt;        if(d>=w)        {            printf("1\n");            continue;        }        for(i=0;i<=n+1;i++) dis[i]=0;        for(i=1;i<=n;i++)        {            if(!dis[i]&&y[i]<=d)            {                if(find(i)) break;            }        }        if(dis[n+1]==0)        {            puts("IMPOSSIBLE");            continue;        }        memset(p,-1,sizeof(p));        eid=0;        ans=0;        src=0;dest=2*(n+m+2)*n+1;        for(k=1;k<=n+m+2;k++)        {            for(i=1;i<=n;i++)            {                if(y[i]<=d)                    addedge(0,(2*k-2)*n+i,inf);                if(w-y[i]<=d)                    addedge((2*k-1)*n+i,dest,inf);                addedge((2*k-2)*n+i,(2*k-1)*n+i,c[i]);            }            if(k>1)            {                for(i=1;i<=n;i++)                    for(j=1;j<=n;j++)                    {                        if((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j])<=d*d)                            addedge((2*k-3)*n+i,(2*k-2)*n+j,inf);                    }            }            if(Dinic()>=m) break;        }            printf("%d\n",k+1);    }}