HDU 6411 带劲的and和【枚举贡献】

HDU 6411 带劲的and和；
分析：写了好久了，突然想补一下博客。。。
公式看完过后，就是求一个联通图内任意两点之间的max(vi,vj)×(vi&vj)。直观感受是n*n的做法，从枚举贡献的角度出发。对于任意的vi，它所带来的贡献是：1，比它小的vj的max(vi,vj)=vi，(vi&vj)=排序后的二进制累加和；2，比它小的vj的max(vi,vj)=vj，对应的vi二进制位累加，对之后的vj进行贡献；
你把一个联通图的节点权值用二进制表示画出来，考虑二进制贡献。

#pragma GCC optimize ("O3")
#pragma GCC optimize ("O2")
#include <bits/stdc++.h>
#include <ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define met(s) memset(s, 0, sizeof(s))
#define rep(i, a, b) for(int i = a; i <= b; ++i)
template <class T> inline void scan_d(T &ret) {
char c; ret = 0;
while ((c = getchar()) < '0' || c > '9');
while (c >= '0' && c <= '9') { 
ret = ret * 10 + (c - '0'), c = getchar();}}
typedef long long LL;
typedef unsigned long long ull;
typedef pair<int, int> pii;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const int MAXN = 1e5 + 10;
vector<int> G[MAXN], V[MAXN];
int dp[35], vi[MAXN], vis[MAXN];
int p;

inline void dfs(int x) {
    V[p].push_back(vi[x]);
    for(int i = 0; i < G[x].size(); ++i) {
        int u = G[x][i];
        if(vis[u]) continue;
        vis[u] = 1;
        dfs(u);
    }
}

inline void init() {
    for(int i = 0; i < MAXN; ++i) G[i].clear(), V[i].clear();
    memset(vis, 0, sizeof(vis));
}

int main() {
    int T;
    scanf("%d", &T);
    while(T--) {
        int n, m; init();
        scanf("%d %d", &n, &m);
        for(int i = 1; i <= n; ++i) scanf("%d", &vi[i]);
        for(int i = 1; i <= m; ++i) {
            int u, v;
            scanf("%d %d", &u, &v);
            G[u].push_back(v);
            G[v].push_back(u);
        }
        p = 0;
        for(int i = 1; i <= n; ++i) {
            if(vis[i]) continue;
            vis[i] = 1;
            dfs(i);
            sort(V[p].begin(), V[p].end());
            p++;
        }
        LL ans = 0;
        for(int i = 0; i < p; ++i) {
            memset(dp, 0, sizeof(dp));
            for(int j = 0; j < V[i].size(); ++j) {
                int u = V[i][j];
                for(int k = 0; k <= 30; ++k) {
                    if((1 << k) & u) {
                        ans += (1 << k) * 1LL * u % mod * 1LL * dp[k];
                        ans %= mod;
                        dp[k]++;
                    }
                }
            }
        }
        printf("%lld\n", ans % mod);
    }
    return 0;
}