match

题目描述

现在有一个关于数字匹配的概念：对于两个正整数x,y，把x,y分别转换成二进制数，并将它们非前导零的部分取出当做一个字符串，把这两个字符串连起来形成一个新的字符串，如果新字符串中0的个数和1的个数之差的绝对值恰好为1，则认为这两个数字是匹配的。

现在给你下标为1到n的n个数a[1],a[2],……,a[n]，你可以随机取i,j(1≤i＜j≤n)，请问a[i]与a[j]匹配的概率是多少,把答案对于1e9+7进行取模

注：分数对1e9+7取模的含义是对于一个b≠0的不可约分数a/b,存在q使
b×qmod(109+7)=a(10^{9}+7),mod指取模运算，q即对1e9+7取模的结果

输入描述:

第一行一个数Q(1≤Q≤5)表示组数
每组第一行一个数n(2≤n≤200000)
第二行n个数，表示a[i](1≤i≤n)，1≤a[i]≤1e9

输出描述:

答案有Q行，每行表示对应的结果

示例1

输入

1
5
1 2 3 4 5

输出

800000006

说明

分数形式为2/5

#include<bits/stdc++.h>
using namespace std;
#define ll long long 
const ll int MOD=1e9+7;
int num[100];
ll qpow(ll a,ll b)//快速幂 
{
	ll result=1;
	while(b)
	{
		if(b&1)
			result=result*a%MOD; 
		a=(a*a)%MOD;
		b>>=1;
	}
	return result;
}
ll sum(ll x)//计数1，0 
{
	ll a1=0,a0=0;
	while(x)
	{
		if(x%2)
			a1++;
		else
			a0++;
		x/=2;
	}
	return a1-a0+30;
}
int main()
{
	ll q;
	cin>>q;
	while(q--)
	{
		int n,x;
		scanf("%d",&n);
		memset(num,0,sizeof(num));//初始化 
		for(int i=1;i<=n;i++)//统排 
		{
			scanf("%d",&x);
			num[sum(x)]++;
		}
		ll cy=0,c=0;//可行方案数，总数 
		for(int i=1;i<100;i++)
		{
			for(int j=i+1;j<100;j++)//注意避免数组越界！ 
			{
				if(num[i]&&num[j])
					if(fabs(i+j-60)==1)
						cy+=num[i]*num[j];
				cy%=MOD;
			}
		}
		for(int i=n-1;i>=1;i--)
			c+=i;
		c%=MOD;
		ll q=(cy*qpow(c,(MOD-2)))%MOD;//分数取模公式 
		printf("%lld\n",q);
	}
	return 0;
}