UVALive 8275 Annual Congress of MUD（网络流+建图）

题目链接：https://cn.vjudge.net/problem/UVALive-8275

题目大意：举办一个持续D天的聚会，每个人可以在[x,y]这个时间区间来，每个人都必须参加聚会，且只能参加一次一天，boss需要让聚会人最多的那一天人最少，设这个数为T，给出n个人的时间区间，如果某个人来了以后T会增加，就输出这个人的编号

一开始往网络流的方向想了，但是也没想出个所以然，自己还是菜啊...

对于每一个区间，我们将它与区间内的每一个点连边，容量为inf，比如区间[1,3]就与1,2,3相连。我们就把它当成个类似二分图的东西，左边是区间，右边是天。天与汇点相连，容量最开始为0。然后我们循环n，对于当前这个人，我们将源点与其区间相连，如果这个区间已经与源点相连了（前面有人和他是一样的区间），就直接在这条边加容量就行了，跑一遍最大流（最大流要累加起来），如果最大流不等于当前的人数，就表示必须要有人在同一天了，然后记录下这个人，将天数与汇点相连的边的容量都加一，因为要的是最大的最小，所以当有多人在同一天的时候，别的天数也能容纳这么多人了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
using namespace std;
const int maxn=1e3+7;
const int maxm=1e6+7;
const int inf=0x3f3f3f3f;
struct Node
{
	int to;
	int capa;
	int next;
}edge[maxm];
int cnt;
int source,sink;
int head[maxn];
int x[10010];
int y[10010];
int node_idx[110][110];
bool vis[110][110];
int num[110][110];//index of edge
int dep[maxn];
int num_to_sink[110][110];
void init()
{
	cnt=0;
	memset(head,-1,sizeof(head));
	memset(vis,false,sizeof(vis));
	memset(num,-1,sizeof(num));
	memset(num_to_sink,-1,sizeof(num_to_sink));
	return;
}
void add(int u,int v,int capa)
{
	edge[cnt].to=v;
	edge[cnt].capa=capa;
	edge[cnt].next=head[u];
	head[u]=cnt++;
	edge[cnt].to=u;
	edge[cnt].capa=0;
	edge[cnt].next=head[v];
	head[v]=cnt++;
	return;
}
bool bfs()
{
	queue<int> que;
	memset(dep,-1,sizeof(dep));
	dep[source]=0;
	que.push(source);
	while(!que.empty())
	{
		int node=que.front();
		que.pop();
		for(int i=head[node];~i;i=edge[i].next)
		{
			int v=edge[i].to;
			if(edge[i].capa>0&&dep[v]==-1)
			{
				dep[v]=dep[node]+1;
				if(v==sink) break;
				que.push(v);
			}
		}
	}
	return dep[sink]!=-1;
}
int dfs(int node,int minn)
{
	if(node==sink||minn==0)
	{
		return minn;
	}
	int r=0;
	for(int i=head[node];~i;i=edge[i].next)
	{
		int v=edge[i].to;
		if(edge[i].capa>0&&dep[v]==dep[node]+1)
		{
			int tmp=dfs(v,min(edge[i].capa,minn));
			if(tmp>0)
			{
				edge[i].capa-=tmp;
				edge[i^1].capa+=tmp;
				r+=tmp;
				minn-=tmp;
				if(minn==0) break;
			}
		}
	}
	if(!r) dep[node]=-1;
	return r;
}
int dinic()
{
	int maxflow=0;
	while(bfs())
	{
		maxflow+=dfs(source,inf);
	}
	return maxflow;
}
int main()
{
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	int n,d;
	while(~scanf("%d",&n))
	{
		if(!n) break;
		init();
		scanf("%d",&d);
		source=0;
		for(int i=0;i<n;i++)
		{
			scanf("%d%d",&x[i],&y[i]);
		}
		int node_cnt=1;
		for(int i=1;i<=d;i++)
		{
			for(int j=i;j<=d;j++)
			{
				node_idx[i][j]=node_cnt++;
			}
		}
		sink=node_cnt+1;
		for(int i=1;i<=d;i++)
		{
			for(int j=i+1;j<=d;j++)
			{
				for(int k=i;k<=j;k++)
				{
					if(i==j) continue;
					add(node_idx[i][j],node_idx[k][k],inf);
				}
			}
		}
		int st=cnt;
		int ldj=cnt;
		for(int i=1;i<=d;i++)
		{
			num_to_sink[i][i]=cnt;
			add(node_idx[i][i],sink,0);
		}
		vector<int> ans;
		int maxflow=0;
		for(int i=0;i<n;i++)
		{
			if(!vis[x[i]][y[i]])
			{
				vis[x[i]][y[i]]=true;
				num[x[i]][y[i]]=cnt;
				add(source,node_idx[x[i]][y[i]],1);
			}
			else
			{
				edge[num[x[i]][y[i]]].capa++;
			}
			maxflow+=dinic();
			if(maxflow!=i+1)
			{
				ans.push_back(i+1);
				for(int j=1;j<=d;j++)
				{
					edge[num_to_sink[j][j]].capa++;
				}
			}
		}
		int len=ans.size();
		printf("%d",ans[0]);
		for(int i=1;i<len;i++)
		{
			printf(" %d",ans[i]);
		}
		puts("");
	}
	return 0;
}