[bzoj3979] [WF2012]infiltration

　　膜了半天ccz大爷的题解...随机调整好神啊。。

　　大概就是把点随机一下，按随机后的顺序控制，顺便删掉无用点。

　　做完一次后，把一个已被控制的点丢回去，然后继续随机顺序做下去。。。

　　为啥我要尝试的次数比ccz大爷多那么多啊QAQ

　　实测一次丢多个点回去反而不好...

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
const int maxn=79;
int st[maxn],b[maxn];
int mp[maxn][maxn],cd[maxn],deg[maxn];
int i,j,k,n,m,ans;
char s[maxn];
 
 
inline void run(){
    int i,top=n,now,numb=0;register int j;
    memset(deg,0,(n+1)<<2);
    for(i=1;i<=n;i++)
        for(st[i]=i,j=1;j<=cd[i];j++)deg[mp[i][j]]++;
    for(int t=20;t;t--){
        for(i=2;i<=top;i++)swap(st[i],st[rand()%(i-1)+1]);
        for(i=1;i<=top;i++){
            now=st[i];
            for(j=1;j<=cd[now];j++)if(deg[mp[now][j]]<=1)break;
            if(j>cd[now]){
                for(j=1;j<=cd[now];j++)deg[mp[now][j]]--;
                b[++numb]=now;
                st[i]=st[top],i--,top--;
            }
        }
        for(i=1;i<=n;i++)if(!deg[i])return;
        if(top<ans)ans=top;
        if(numb){
            now=rand()%numb+1;
            for(j=1;j<=cd[b[now]];j++)deg[mp[b[now]][j]]++;
            st[++top]=b[now];
            b[now]=b[numb],numb--;
        }
    }
}
int main(){
    int T=0;srand(1844677);
    while(scanf("%d",&n)!=EOF){
        T++;
        for(i=1;i<=n;i++){
            scanf("%s",s+1);s[i]='1';
            cd[i]=0;
            for(j=1;j<=n;j++)if(s[j]=='1')mp[i][++cd[i]]=j;
        }
        ans=n;
        for(i=100;i;i--)run();
        printf("Case %d: %d\n",T,ans);
    }
}