2nd SCPC - K 无限矩阵(二进制规律+前缀和)

在这里插入图片描述
首先找每一列的周期规律,将矩阵按正常二进制来看是这样的:
在这里插入图片描述
然后我们把最右边看着第一列,依次向左递增

每一列的周期为 2 i 2^i 2i,且前 2 i − 1 2^{i-1} 2i1个为 0 0 0,后 2 i − 1 2^{i-1} 2i1个为 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(u1)

s u m ( x ) sum(x) sum(x)的计算方法为:

先计算出该行前有 x x x个周期,接着对周期取模,看答案 r e s res res超出半个周期 2 i − 1 2^{i-1} 2i1多少,就有 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) x2i1+(res>2i1?res2i1: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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值