位运算 - 异或求和

 异或求和

Problem's Link:  http://acm.xmu.edu.cn/JudgeOnline/problem.php?id=1254

---

 

Mean: 

Description
给你2个区间[A,B]和[C,D],现在只要求你求出区间[A,B]和[C,D]内任意2个整数异或后的和，因为答案可能会很大，你只需将结果%mod即可。

For(i:A→B)
      For(j:C→D)
           Sum += (i^j);

Input
输入第一行为T(T<=15),表示测试的数据组数，第2行到T+1行，每行5个数字，分别为A,B,C,D,mod.
(1<=A,B,C,D<2^31,A<=B,C<=D,1<=mod& lt;=1,000,000,007),即为上述的数值。
 

Output
输出有T行，每行有一个数字，代表题述的答案（即Sum % mod）。


Sample Input
3
5 11 9 12 43
9 12 5 11 83
1 1 2 2 3

Sample Output
18
24
0

 

analyse:

看到这么大的输入，再看看时限是500ms，开始有点摸不着头脑，看来只能用log(n)的方法来做了。

典型的位运算题目。

首先我们求出[A,B]区间中所有数各个位上1的个数。比如[1~8]就是1,4,4,4；

怎么求呢？

我们可以打表找一下规律：

1: 0001

2: 0010

3: 0011

4: 0100

5: 0101

6: 0110

7: 0111

8: 1000

9: 1001

.....

我们会发现，[1,n]中第k位总共1的个数就是：n/(k<<1) *(k<<1 /2 )。或者说：符合 x mod 2^k >=2^（k-1） 的 x 第k位就是1。

额，这个要自己去找规律。

要求[A~B]区间内各位上1的个数，只需要B1[i]-A1[i]就行。

同理，我们可以求出[C~D]区间上各位1的个数。

剩下的就是怎么求和了。

求和比较简单：

区间【A,B】第k位为1的个数乘以区间【C,D】第k位为0的个数+区间【A,B】第k位为0的个数乘以区间【C,D】第k位为1的个数）*2^(k-1)。

 

Time complexity: O(logn)

 

Source code: 

/*
* this code is made by crazyacking
* Verdict: Accepted
* Submission Date: 2015-08-05-11.38
* Time: 0MS
* Memory: 137KB
*/
#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#define  LL long long
#define  ULL unsigned long long
using namespace std;

void get1( LL n , LL arr []) { /**< 求1~n所有数各个位上1的个数之和 */
      for( int i = 1; i < 32; ++ i) arr [ i ] = 0;
      for( LL i = 1 , f = 2 ,b = 1; b <=n; ++ i , f <<= 1 ,b <<= 1) {
            arr [ i ] =(n / f) *( f / 2); /**< 对于每一位,1~n可分为n/t组(t是每组的01数量)，其中每组有t/2个是1 */
            if(n % f >=b) arr [ i ] +=(n % f -b + 1); /**< 加上余数部分 */
      }
}

int main() {
      ios_base :: sync_with_stdio( false);
      cin . tie( 0);
      LL n , t;
      cin >> t;
      while( t --) {
            LL A ,B , C , D , Mod , a1 [ 32 ], b1 [ 32 ], c1 [ 32 ], d1 [ 32 ];
            cin >> A >>B >> C >> D >> Mod;
            get1( A - 1 , a1); get1(B , b1); get1( C - 1 , c1); get1( D , d1);
            LL ab1 [ 32 ], ab0 [ 32 ], cd1 [ 32 ], cd0 [ 32 ];
            for( int i = 1; i < 32; ++ i) ab1 [ i ] = b1 [ i ] - a1 [ i ], cd1 [ i ] = d1 [ i ] - c1 [ i ];
            int l1 =B - A + 1 , l2 = D - C + 1;
            for( int i = 1; i < 32; ++ i) ab0 [ i ] = l1 - ab1 [ i ], cd0 [ i ] = l2 - cd1 [ i ];
            LL sum = 0 , base = 1;
            for( int i = 1; i < 32; ++ i)
                  sum =( sum +( base <<( i - 1)) % Mod *( ab0 [ i ] % Mod * cd1 [ i ] % Mod + ab1 [ i ] % Mod * cd0 [ i ] % Mod) % Mod) % Mod;
            cout << sum << endl;
      }
      return 0;
}
/*

*/