题意：
给你n、m、d，代表有n个城市。m条城市之间的关系，每一个城市要在日后d天内都有电。
对于每一个城市，都有一个发电站，每一个发电站能够在[a,b]的每一个连续子区间内发电。
x城市发电了。他相邻的城市也有电。而且每一个发电站仅仅能启动一次。或者不启动。
如今问，怎样安排发电站的启动。保证每一个城市d天都有电。
输出发电方案。不发电的话输出0 0
思路：
一个简单的精确覆盖问题。就是建图比較麻烦一点。
这里考虑到每天都要得到电。所以把每一个城市每天都设为列（n*d）
然后每一个城市对于[a,b]的全部子区间，作为行。
然后对于每一个城市仅仅能启动1次，所以须要加入n行，代表每一个城市发过点没。
所以列是n*d+n。
还须要注意的是，在推断dance成立的时候，仅仅要推断R[0]==0 || R[0]>n*d
由于每一个发电站不一定要发电。
代码：
#include"stdio.h"
#include"algorithm"
#include"string.h"
#include"iostream"
#include"queue"
#include"map"
#include"vector"
#include"string"
using namespace std;
#define N 2010*(66*6)
#define M 2010
int ooo;
struct DLX
{
int n,m,C;
int U[N],D[N],L[N],R[N],Row[N],Col[N];
int H[M],S[M],cnt,ans[M];
void init(int _n,int _m)
{
n=_n;
m=_m;
for(int i=0; i<=m; i++)
{
U[i]=D[i]=i;
L[i]=(i==0?m:i-1);
R[i]=(i==m?0:i+1);
S[i]=0;
}
C=m;
for(int i=1; i<=n; i++) H[i]=-1;
}
void link(int x,int y)
{
C++;
Row[C]=x;
Col[C]=y;
S[y]++;
U[C]=U[y];
D[C]=y;
D[U[y]]=C;
U[y]=C;
if(H[x]==-1) H[x]=L[C]=R[C]=C;
else
{
L[C]=L[H[x]];
R[C]=H[x];
R[L[H[x]]]=C;
L[H[x]]=C;
}
}
void del(int x)
{
R[L[x]]=R[x];
L[R[x]]=L[x];
for(int i=D[x]; i!=x; i=D[i])
{
for(int j=R[i]; j!=i; j=R[j])
{
U[D[j]]=U[j];
D[U[j]]=D[j];
S[Col[j]]--;
}
}
}
void rec(int x)
{
for(int i=U[x]; i!=x; i=U[i])
{
for(int j=L[i]; j!=i; j=L[j])
{
U[D[j]]=j;
D[U[j]]=j;
S[Col[j]]++;
}
}
R[L[x]]=x;
L[R[x]]=x;
}
int dance(int x)
{
if(R[0]==0 || R[0]>ooo) //判定条件
{
cnt=x;
return 1;
}
int now=R[0];
for(int i=R[0]; i!=0 && i<=ooo; i=R[i]) //这里注意也要改
{
if(S[i]<S[now]) now=i;
}
del(now);
for(int i=D[now]; i!=now; i=D[i])
{
ans[x]=Row[i];
for(int j=R[i]; j!=i; j=R[j]) del(Col[j]);
if(dance(x+1)) return 1;
for(int j=L[i]; j!=i; j=L[j]) rec(Col[j]);
}
rec(now);
return 0;
}
} dlx;
struct node
{
int a,b;
} kx[66];
struct answer
{
int id,x,y;
}daan[2002];
int main()
{
int n,m,d;
while(scanf("%d%d%d",&n,&m,&d)!=-1)
{
ooo=n*d;
int v[66][66];
memset(v,0,sizeof(v));
while(m--)
{
int a,b;
scanf("%d%d",&a,&b);
v[a][b]=v[b][a]=1;
}
int cnt=0;
for(int i=1; i<=n; i++)
{
scanf("%d%d",&kx[i].a,&kx[i].b);
for(int j=kx[i].a; j<=kx[i].b; j++) for(int k=j; k<=kx[i].b; k++) cnt++;
}
dlx.init(cnt,n*d+n);
cnt=0;
for(int i=1; i<=n; i++)
{
for(int j=kx[i].a; j<=kx[i].b; j++)
{
for(int k=j; k<=kx[i].b; k++)
{
++cnt;
daan[cnt].id=i;
daan[cnt].x=j;
daan[cnt].y=k;
for(int l=j; l<=k; l++)
{
dlx.link(cnt,(i-1)*d+l);
for(int u=1;u<=n;u++)
{
if(v[i][u])
{
//printf("%d %d %d
",i,u,l);
dlx.link(cnt,(u-1)*d+l);
}
}
}
dlx.link(cnt,n*d+i);
}
}
}
int ans=dlx.dance(0);
if(ans==0) puts("No solution");
else
{
int us[66];
memset(us,0,sizeof(us));
for(int i=0;i<dlx.cnt;i++) us[daan[dlx.ans[i]].id]=dlx.ans[i];
for(int i=1;i<=n;i++)
{
if(us[i]==0) printf("0 0
");
else printf("%d %d
",daan[us[i]].x,daan[us[i]].y);
}
}
puts("");
}
return 0;
}