首页 > 代码库 > upcoj2673 It Can Be Arranged(isap)
upcoj2673 It Can Be Arranged(isap)
题意:
有N节课,每节课有起止时间和学生数
然后给你M是每个教室容纳的学生数
然后给你n*n的矩阵表示上完第i节课然后上第j节课需要a[i][j]的时间调整
第i节课结束时间加上调整时间要严格小于第j节课的开始时间
问你最少需要多少个教室
思路:
很简单的一个最大流
每节课分为两个节点,0为超级源点,2*n+1为超级汇点
先算出每节课需要的教室数,加到all里
然后从0到1-n加一个流量为教室数的边
从n+1——2*n到2*n+1加一个流量为教室数的边
然后将满足条件的i和j课程,从i到j+n连一条流量为inf的边
答案就是all-最大流
/* ***********************************************Author :devil************************************************ */#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>#include <vector>#include <queue>#include <set>#include <stack>#include <map>#include <string>#include <cmath>#include <stdlib.h>#define inf 0x3f3f3f3f#define LL long long#define rep(i,a,b) for(int i=a;i<=b;i++)#define dec(i,a,b) for(int i=a;i>=b;i--)#define ou(a) printf("%d\n",a)#define pb push_back#define mkp make_pairtemplate<class T>inline void rd(T &x){char c=getchar();x=0;while(!isdigit(c))c=getchar();while(isdigit(c)){x=x*10+c-‘0‘;c=getchar();}}#define IN freopen("in.txt","r",stdin);#define OUT freopen("out.txt","w",stdout);using namespace std;const int mod=1e9+7;const int N=1e5+10;int t,n,m;#define maxm 500010#define maxn 222int a[maxn],b[maxn],s[maxn],all;int head[maxn],eid;int dis[maxn];//残量网络中节点i到汇点t的最短距离int num[maxn];//和t的最短距离等于i的节点数量int cur[maxn];//当前弧下标int pre[maxn];//可增广路上的上一条弧的编号struct node{ int v,cap,next;} edge[maxm];void add(int u,int v,int cap){ edge[eid].v=v; edge[eid].cap=cap; edge[eid].next=head[u]; head[u]=eid++; edge[eid].v=u; edge[eid].cap=0; edge[eid].next=head[v]; head[v]=eid++;}void bfs(int source,int sink)//预处理,利用反向BFS,更新dis数组{ queue<int>q; while(!q.empty()) q.pop(); memset(num,0,sizeof(num)); memset(dis,-1,sizeof(dis)); q.push(sink); dis[sink]=0; num[0]=1; while(!q.empty()) { int u=q.front(); q.pop(); for(int i=head[u]; i!=-1; i=edge[i].next) { int v=edge[i].v; if(dis[v]==-1) { dis[v]=dis[u]+1;//找允许弧 num[dis[v]]++; q.push(v); } } }}int isap(int source,int sink,int n)//n为残量网络中的节点到汇点的最大距离,通常节点的个数,即上限{ memcpy(cur,head,sizeof(cur)); int flow=0,u=pre[source]=source; bfs(source,sink);//更新dis数组 while(dis[source]<n) { if(u==sink) { int df=inf,pos; for(int i=source;i!=sink;i=edge[cur[i]].v)//追踪增广路路径,最小残量df { if(df>edge[cur[i]].cap) { df=edge[cur[i]].cap; pos=i; } } for(int i=source;i!=sink;i=edge[cur[i]].v)//更新流量 { edge[cur[i]].cap-=df; edge[cur[i]^1].cap+=df; } flow+=df; u=pos; } int st; for(st=cur[u];st!=-1;st=edge[st].next)//从当前弧开始查找允许弧 if(dis[edge[st].v]+1==dis[u]&&edge[st].cap)//找到允许弧跳出 break; if(st!=-1) { cur[u]=st; pre[edge[st].v]=u; u=edge[st].v; } else { if((--num[dis[u]])==0) break;//GAP优化,出现断层结束 int mind=n; for(int id=head[u];id!=-1;id=edge[id].next)//retreat操作:更新dis数组 { if(mind>dis[edge[id].v]&&edge[id].cap) { cur[u]=id;//修改标号的同时修改当前弧 mind=dis[edge[id].v]; } } dis[u]=mind+1; num[dis[u]]++; if(u!=source) u=pre[u];// 回溯继续寻找允许弧 } } return flow;}void init(){ memset(head,-1,sizeof(head)); eid=0; all=0;}int main(){ rd(t); rep(y,1,t) { init(); rd(n),rd(m); rep(i,1,n) rd(a[i]),rd(b[i]),rd(s[i]); rep(i,1,n) { int w=s[i]/m; if(s[i]%m) w++; all+=w; add(0,i,w); add(i+n,n*2+1,w); } rep(i,1,n) rep(j,1,n) { rd(m); if(b[i]+m<a[j]) add(i,j+n,inf); } printf("Case %d: %d\n",y,all-isap(0,2*n+1,n)); } return 0;}
upcoj2673 It Can Be Arranged(isap)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。