2017 计蒜之道 初赛 第三场 腾讯狼人杀（简单)(暴搜)

传送门

不久前，腾讯推出了一款改进后的狼人杀游戏。如下是这款游戏的背景和详细规则：

很久很久以前，在莱茵河畔一座岸壁陡峭的山顶上，有一个名叫“杜斯特瓦德”的小村庄。不久前，这个小村庄每晚都会受到狼人的侵袭。每个夜晚，狼人都会在村中进行抢劫，并且会有一个村民成为这群狼人的牺牲品。然而村民们不会坐以待毙，他们试图在白天找到这些狼人并处决他们。

作为杜斯特瓦德村的长老，你需要组织村民来对抗狼人，村民包括许多有超能力的人，比如预言家可以检验一个人是否是隐藏在神民阵营的狼人、猎人可以开枪打死一名狼人。因为狼人白天隐藏的非常深，仅仅单人的行动往往不能限制狼人的行动，多人行动将会起到更好的效果。比如预言家和猎人一起行动将被检验出是狼人的“村民”直接打死，村庄将回归和平。

我们将这样的一对人所产生的战斗力记作 $w$，值得注意的是，一个人可以被多次计算，比如预言家可以先后和猎人、女巫（或其他人）组合，即如果预言家和猎人的战斗力为 $w$，预言家和女巫的战斗力为 w'w​′​​，这三个人的战斗力之和为 w+w'w+w​′​​。

不过因为经过“新月”事件，村庄贫苦破败，你需要支付一定的金额给各个被你聘用的村民，假设你邀请了k个村民，则需要 $k(2n-k)$块钱（$n$ 为村民总数），你的资金也不太足够，所以你希望资金的利用率最大，资金的利用率定义为：$($被选择的人战斗力之和 $/$ 你支付的资金$)$。除此之外有些人基于他们的重要性，是必须被选择的。你必须邀请他们来参加战斗（保证至少一人）。

输入格式

第一行包括两个整数 $n,m$ 分别表示 $n$ 个村民和有 $m$ 对人之间有战斗力。

接下来 $m$ 行，每行有三个整数 u_i,v_i,w_iu​i​​,v​i​​,w​i​​ (1 \le u_i,v_i \le n,(1≤u​i​​,v​i​​≤n, 1 \le w_i \le 100,1≤w​i​​≤100, u_i \neq v_i)u​i​​≠v​i​​)表示 这对人可以产生的战斗力 w_iw​i​​。

接下来 $n$ 个数，对于第 $i$ 个数，$0$ 表示第 $i$ 个人可参加战斗，$1$ 表示第 $i$ 个人必须参加战斗（保证至少一个 $1$）。

对于简单版本，$1\le n\le 20,$ $0\le m\le n(n-1)/2$；

对于中等版本，$1\le n\le 300,$ $0\le m\le \min (1000,n(n-1)/2)$；

对于困难版本，$1\le n\le 400,$ $0\le m\le \min (10000,n(n-1)/2)$。

输出格式

对于简单版本和中等版本，输出一个 $4$ 位小数，表示资金利用率的最大值，结果四舍五入；

对于困难版本，输出一个 $6$ 位小数，表示资金利用率的最大值，结果四舍五入。

样例输入1

3 11 2 31 0 0

样例输出1

0.3750

样例输入2

3 11 2 30 0 1

样例输出2

0.3333

#include <iostream>#include <stdio.h>#include <cstring>#include <string>#include <algorithm>#include <queue>#include <stack>#include <cmath>#include <map>#include <bitset>#include <set>#include <vector>#include <functional>using namespace std;#define pi acos(-1)#define endl '\n'#define rand() srand(time(0));#define me(x) memset(x,0,sizeof(x));#define foreach(it,a) for(__typeof((a).begin()) it=(a).begin();it!=(a).end();it++)#define close() ios::sync_with_stdio(0);typedef long long LL;const int INF=0x3f3f3f3f;const LL LINF=0x3f3f3f3f3f3f3f3fLL;//const int dx[]={-1,0,1,0,-1,-1,1,1};//const int dy[]={0,1,0,-1,1,-1,1,-1};const int maxn=1e4+5;const int maxx=1e6+3;const double EPS=1e-7;const int MOD=10000007;typedef pair<int, int>P;#define mod(x) ((x)%MOD);template<class T>inline T min(T a,T b,T c) { return min(min(a,b),c);}template<class T>inline T max(T a,T b,T c) { return max(max(a,b),c);}template<class T>inline T min(T a,T b,T c,T d) { return min(min(a,b),min(c,d));}template<class T>inline T max(T a,T b,T c,T d) { return max(max(a,b),max(c,d));}//typedef tree<pt,null_type,less< pt >,rb_tree_tag,tree_order_statistics_node_update> rbtree;/*lch[root] = build(L1,p-1,L2+1,L2+cnt);    rch[root] = build(p+1,R1,L2+cnt+1,R2);中前*//*lch[root] = build(L1,p-1,L2,L2+cnt-1);    rch[root] = build(p+1,R1,L2+cnt,R2-1);中后*/long long gcd(long long a , long long b){if(b==0) return a;a%=b;return gcd(b,a);}#define FOR(x,n,i) for(int i=x;i<=n;i++)#define FOr(x,n,i) for(int i=x;i<n;i++)#define W whileint n,m,sum=0,num=0;int vis[maxn],mapp[30][30],cnt[30];double ans=0;void dfs(int x){    if(x==n+1)    {       // cout<<num<<endl;        double val=sum*1.0/(double)(num*(2*n-num));//取效率大值        if(ans<val) ans=val;        return ;    }    FOR(0,1,i)//0表示不取 1表示取    {        if(vis[x]&&i==0) continue;        cnt[x]=i;        if(cnt[x])        {            FOr(1,x,j)//只取在x前面的农民            {                if(cnt[j])                    sum+=mapp[j][x];              //  cout<<j<<"  "<<x<<endl;            }            num++;        }        dfs(x+1);        if(cnt[x])//回溯        {            FOr(1,x,j)                if(cnt[j])                    sum-=mapp[j][x];            num--;        }    }}int main(){    cin>>n>>m;    me(mapp);me(vis);me(cnt);    FOR(1,m,i)    {        int x,y,z;        cin>>x>>y>>z;        mapp[x][y]=z;mapp[y][x]=z;    }    FOR(1,n,i)        cin>>vis[i];    dfs(1);    printf("%.4lf\n",ans);}