题意：
有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 ",a) #define pb push_back #define mkp make_pair template<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 222 int 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 ",y,all-isap(0,2*n+1,n)); } return 0; }