HDU 5036 Explosion

考虑爆破顺序为一个序列，
对于第i扇门，有x扇门的爆破能够直接或间接打开第i扇门，
那么第i扇门被爆破，当且仅当序列中i排在这x个数的最前面，
因此第i扇门的爆破期望值是1/x，
然后每一扇门的期望都是独立的，
那么把所有1/x加起来就是答案...
那么具体实现需要统计每一个xi，
那么可以先强连通缩点得到一个DAG，
然后bitset暴力搞一下，
时间复杂度大约n^3/64...

#include <cstdio>
#include <bitset>
#include <algorithm>
#include <iostream>

using namespace std;

const int MAX = 1010;

int T,k,N,num;
bitset<MAX> key[MAX];

int main(void)
{
//freopen("input.txt","r",stdin);
int cas = 1;
scanf("%d", &T);
while(T--){
scanf("%d", &N);
for(int i = 0; i < N; ++i){
key[i].reset();
key[i][i] = true;
}
for(int i = 0; i < N; ++i){
scanf("%d", &k);
for(int j = 0; j < k; ++j){
scanf("%d", &num);
key[i][num-1] = true;
}
}

for(int j = 0; j < N; ++j)
for(int i = 0; i < N; ++i)
if(key[i][j])
key[i] |= key[j];

double ans = 0.0;
for(int j = 0; j < N; ++j){
int cnt = 0;
for(int i = 0; i < N; ++i)
if(key[i][j]) cnt++;
ans += 1.0 / cnt;
//printf("%f\n",ans);
}

printf("Case #%d: %.5f\n",cas++,ans);
}
return 0;
}