首先找每一列的周期规律,将矩阵按正常二进制来看是这样的:
然后我们把最右边看着第一列,依次向左递增
每一列的周期为 2 i 2^i 2i,且前 2 i − 1 2^{i-1} 2i−1个为 0 0 0,后 2 i − 1 2^{i-1} 2i−1个为 1 1 1
只要给出起始的第 u u u行和第 d d d行,那么第 i i i列在这部分区间的1的个数等于 s u m ( d ) − s u m ( u − 1 ) sum(d)-sum(u-1) sum(d)−sum(u−1)
而 s u m ( x ) sum(x) sum(x)的计算方法为:
先计算出该行前有 x x x个周期,接着对周期取模,看答案 r e s res res超出半个周期 2 i − 1 2^{i-1} 2i−1多少,就有 x ∗ 2 i − 1 + ( r e s > 2 i − 1 ? r e s − 2 i − 1 : 0 ) {x*2^{i-1} } + (res>2^{i-1}?res-2^{i-1}:0) x∗2i−1+(res>2i−1?res−2i−1:0)
但是可能有 1 e 9 1e9 1e9列,枚举肯定超时
那么我们冷静分析一下,不难发现周期是按2的次幂增长的,到了第31行,周期为 2 31 2^{31} 231,半个周期为 2 30 > 1 e 9 2^{30} > 1e^9 230>1e9,因为前半个周期为0,那么当列数超过31之后,已经不会累计答案了,那么特判一下即可
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <math.h>
#include <cstdio>
#include <string>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <unordered_map>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;
const double eps=1e-8;
const double pi=acos(-1.0);
const int inf=0x3f3f3f3f;
const ll INF=1e18;
const int Mod=1e9+7;
const int maxn=2e5+10;
int main(){
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
int t,l,r,u,d;
//cout<<(1<<30)<<endl;
cin>>t;
while(t--){
cin>>l>>r>>u>>d;
if(l>=31) cout<<"0\n";
else{
r=min(r,30);
int ans=0;
for(int i=l;i<=r;i++){
int T=(1<<i),num=(1<<(i-1));
int res1=(u-1)%T,p=(u-1)/T;
int sum1=num*p+(res1>num?res1-num:0);
int res2=d%T,q=d/T;
int sum2=num*q+(res2>num?res2-num:0);
int sum=sum2-sum1;
//cout<<sum<<endl;
if(sum&1) ans++;
}
cout<<ans<<"\n";
}
}
return 0;
}