HDU 3663 舞蹈链之不可重复覆盖

点击打开链接

题意：给一个无向图，然后n个城市的供电范围，每个现在要求每一个城市的D天都可以有电，对于城市A发电，那么与它相邻的所有城市都会有电，但是问题是每个城市一天内只可以被供电一次，否则会坏掉，并且每个城市的供电天数有范围而且每个城市只能开启开关一次，之后不能在使用，也就是说城市A的供电时间必须是连续的，还有就是可以不用这个城市

思路：因为每个城市一天只能被供电一次，那么就是不可重复覆盖的舞蹈链模型了，那么考虑的就是把什么作为列和行了，很明显这道题就是N*D+N作为列，N*D是n个城市的D天都要供电，然后后面的N是保证每个城市直发电一次，然后考虑什么作为行，对于一个城市的发电范围，我们可以截取其中所有的小区间来作为行，当然不用的情况要再加一个0 0，比如1到5的范围便可以变成16行，0 0，1 1，1 2，1 3，1 4，1 5，2 2，2 3，2 4，2 5，3 3，3 4，3 5，4 4，4 5，5 5，对于每一个范围它们连的列就是自己对应的天数的城市，就是一些小细节的处理了，大致就是这样了，对了千万不要忘记重边，重边，重边，主要的事说三遍

#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=1010;
int L[maxn*400],R[maxn*400],U[maxn*400],D[maxn*400];//节点的上下左右四个方向的链表
int C[maxn*400],H[maxn],cnt[maxn],ans[maxn],Row[maxn*400];//C列H行cnt列链表中元素个数
int n,m,id,len,d;
vector<int>G[maxn];
void init(int lll){
    for(int i=0;i<maxn;i++) G[i].clear();
    for(int i=0;i<=lll;i++){
        cnt[i]=0;U[i]=D[i]=i;
        L[i+1]=i;R[i]=i+1;
    }
    R[lll]=0;id=lll+1;
    memset(H,-1,sizeof(H));
}
void Link(int r,int c){
    cnt[c]++;C[id]=c;Row[id]=r;
    U[id]=U[c];D[U[c]]=id;
    D[id]=c;U[c]=id;
    if(H[r]==-1) H[r]=L[id]=R[id]=id;
    else{
        L[id]=L[H[r]];R[L[H[r]]]=id;
        R[id]=H[r];L[H[r]]=id;
    }
    id++;
}
void Remove(int Size){
    L[R[Size]]=L[Size];
    R[L[Size]]=R[Size];
    for(int i=D[Size];i!=Size;i=D[i]){
        for(int j=R[i];j!=i;j=R[j]){
            U[D[j]]=U[j];D[U[j]]=D[j];
            cnt[C[j]]--;
        }
    }
}
void Resume(int Size){
    for(int i=D[Size];i!=Size;i=D[i]){
        for(int j=R[i];j!=i;j=R[j]){
            U[D[j]]=j;D[U[j]]=j;
            cnt[C[j]]++;
        }
    }
    L[R[Size]]=Size;R[L[Size]]=Size;
}
int Dance(int k){
     int c;
     if(!R[0]){
         len=k;
         return 1;
     }
     for(int Min=inf,i=R[0];i;i=R[i])
         if(Min>cnt[i]) Min=cnt[i],c=i;
     Remove(c);
     for(int i=D[c];i!=c;i=D[i]){
         ans[k]=Row[i];
         for(int j=R[i];j!=i;j=R[j]) Remove(C[j]);
         if(Dance(k+1)) return 1;
         for(int j=L[i];j!=i;j=L[j]) Resume(C[j]);
    }
    Resume(c);
    return 0;
}
int A[70],B[70],tmp[maxn][2];
void Linkline(int l,int r,int x,int kkk){
    Link(kkk,n*d+x);
    for(int i=l;i<=r;i++){
        Link(kkk,(i-1)*n+x);
        for(unsigned j=0;j<G[x].size();j++){
            int t=G[x][j];
            Link(kkk,(i-1)*n+t);
        }
    }
}
int visedge[70][70];
int main(){
    int u,v;
    while(scanf("%d%d%d",&n,&m,&d)!=-1){
        init(n*d+n);
        memset(visedge,0,sizeof(visedge));
        for(int i=1;i<=m;i++){
            scanf("%d%d",&u,&v);
            if(visedge[u][v]) continue;
            G[u].push_back(v);
            G[v].push_back(u);
            visedge[u][v]=1;visedge[v][u]=1;
        }
        for(int i=1;i<=n;i++) scanf("%d%d",&A[i],&B[i]);
        int kkk=1;
        for(int i=1;i<=n;i++){
            tmp[kkk][0]=0;tmp[kkk++][1]=0;
            Link(kkk-1,n*d+i);
            for(int j=A[i];j<=B[i];j++){
                for(int l=j;l<=B[i];l++){
                    tmp[kkk][0]=j;tmp[kkk++][1]=l;
                    Linkline(j,l,i,kkk-1);
                }
            }
        }
        int flag=Dance(0);
        if(flag==0) printf("No solution\n");
        else{
            sort(ans,ans+len);
            for(int i=0;i<len;i++) printf("%d %d\n",tmp[ans[i]][0],tmp[ans[i]][1]);
        }
        printf("\n");
    }
    return 0;
}