首页 > 代码库 > lowbit

lowbit

树状数组(lowbit)

Time Limit:1000ms   Memory Limit:128MB

 

题目描述

这天,LYK在学习树状数组。

当它遇到一个叫lowbit的函数时有点懵逼。lowbit(x)的意思是将x分解成二进制,它的值就是,其中k是最小的满足(x & )>0的数。(&是二进制中的and运算)

LYK甚至知道lowbit(x)=(x&-x)。但这并没什么用处。

现在LYK有了n个数字,为了使自己更好的理解lowbit是什么意思。它想对所有n^2个二元组求lowbit。具体的,对于一个二元组(ai,aj),它的值为lowbit(ai xor aj) (xor表示异或的意思),那么总共有n^2对二元组,LYK想知道所有二元组的值加起来是多少。

这个答案可能很大,你只需输出这个值对1000000007取模后的结果就可以了。

 

输入格式(lowbit.in)
第一行一个数n,表示有n个这样的数字。

第二行n个数ai。

 

输出格式(lowbit.out)

一个数表示答案。

 

输入样例

5

1 2 3 4 5

 

输出样例

32

数据范围

对于30%的数据n<=1000。

对于另外10%的数据ai<=1。

对于再另外10%的数据ai<=3。

对于再再另外20%的数据ai<1024。

对于100%的数据1<=n<=100000,0<=ai<2^30。

  先说一下暴力的做法(只能拿70分);

分类处理

  对于前30%的数据,直接按照题目去做就行。

  接下来40%的,a[i]范围较小,开一个桶数组,记录每个数出现了多少次(因为n<=100000,而a[i]<=1024,每个数一定出现过多次,我们只算一次就行了)

  

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<math.h>
using namespace std;
#define M 1000000007
long long n,c;
int f[10205];
long long ans,o,maxn;
long long last=0,x,y,t;
long long     lowbit(long long x)
{
    return x&(-x);
} 
int main()
{
    freopen("lowbit.in","r",stdin);
    freopen("lowbit.out","w",stdout);
    scanf("%lld",&n);
    if(n<=1000) //这里掉了个等号,调了好久。。。。
    {
        for(int i=1;i<=n;i++)
        scanf("%d",&f[i]);
        for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
        o=lowbit((long long)f[i] ^ f[j]),ans=(ans+o)%M;
        ans=ans*2%M;
        cout<<ans;
        return 0;
    }
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&c);
        f[c]++;maxn=max(maxn,c);
    }    
    for(int i=0;i<=maxn;i++)
    for(int j=i+1;j<=maxn;j++)
    {        
        t=(i ^ j);
        o=lowbit(t);
        ans=(ans+o*f[i]*f[j]%M)%M;
    }
    ans=ans*2%M; 
    cout<<ans;
    return 0;
} 

满分的做法就要考虑lowbit 和二进制数的性质了。

  lowbit是指 某数在二进制中最后一个“1”所在位置的权值。

  两个数 xor(异或)的lowbit与 他们从低位到高位第一个不同位的权值有关。

 

思路

  (1)那么,我们就从低位向高位注意比较有多少个在这个位上是“1”,多少个是“0”。分别放在集合X,集合Y中。

  (2)那么对于任意a∈X,b∈Y,,    lowbit(a,b)就是这个位的权值。

  这样集合X和Y之间的lowbit就求出来了(集合X,Y内部的数的lowbit还没求),接下来就分别对集合X,Y内部的数,对下一位进行同样操作 (1)(2)

实现的算法可以是分治,也可以是trie树。

时间原因,分治和trie树的代码先待续。。

lowbit