二分和前缀和专题

一、前缀和

1.一维模板

题目链接：活动 - AcWing

#include<iostream>
using namespace std;

int a[1000005];

int main()
{
    int n,m;
    int l,r;
    cin>>n>>m;
    
    for(int i=1;i<=n;++i)
    {
        cin>>a[i];
        a[i]+=a[i-1];
    }
  
    for(int i=1;i<=m;++i)
    {
        cin>>l>>r;
        cout<<a[r]-a[l-1]<<endl;
    }
}

2.二维模板

题目链接：活动 - AcWing

#include<iostream>
using namespace std;

int map[1005][1005];
int c[1005][1005];

int main()
{
    int n,m,q;
    int x1,x2,y1,y2;
    
    cin>>n>>m>>q;
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=m;++j)
        {
            cin>>map[i][j];
        }
    }
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=m;++j)
        {
            map[i][j]+=map[i-1][j]+map[i][j-1]-map[i-1][j-1];
        }
    }
    
    for(int i=1;i<=q;++i)
    {
        cin>>x1>>y1>>x2>>y2;
        cout<<map[x1-1][y1-1]-map[x1-1][y2]-map[x2][y1-1]+map[x2][y2]<<endl;
    }
}

3.应用

题目链接：99. 激光炸弹 - AcWing题库

注意点：一个点上可能有多个炸弹

AC代码

#include<iostream>
using namespace std;

typedef long long ll;
int n,r;
int a[5005][5005]={0};
int x,y,w;
int ans=0;

int main()
{
    cin>>n>>r;
    
    r=min(r,5001);
    
    for(int i=1;i<=n;++i)
    {
        cin>>x>>y>>w;
        a[x+1][y+1]+=w;//真坑
    }
    
    for(int i=1;i<=5001;++i)
    {
        for(int j=1;j<=5001;++j)
        {
            a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
        }
    }
    
    for(int i=r;i<=5001;i++)
    {
        for(int j=r;j<=5001;j++)
        {
            ans=max(ans,a[i][j]+a[i-r][j-r]-a[i-r][j]-a[i][j-r]);
        }
    }
    
    cout<<ans<<endl;
    
    return 0;
}

练习2：Acwing 1230.K倍区间_skyang.的博客-CSDN博客

二、二分

1.整数二分

题目链接：活动 - AcWing

代码：

#include<iostream>
using namespace std;

int main()
{
    int n,m,q;
    int mid;
    int a[1000005];
    
    cin>>n>>m;
    for(int i=1;i<=n;++i)
    {
        cin>>a[i];
    }
    
    for(int i=1;i<=m;++i)
    {
        cin>>q;
        int l=1,r=n;
        while(l<r)
        {
            mid=(l+r)/2;
            if(q>a[mid]) l=mid+1;
            else r=mid;
        }
        
        if(q==a[l])
        {
            cout<<l-1<<" ";
            r=n;
            while(l<r)
            {
                if(q>=a[mid]) l=mid;
                else r=mid-1;
                mid=(l+r+1)/2;
            }
            cout<<r-1<<endl;
        }
        else cout<<"-1 -1"<<endl;
    }
}

2.实数二分

题目链接：活动 - AcWing

代码：

#include<iostream>
using namespace std;

int main()
{
    double l=-10000,r=10000;
    double num,mid;
    cin>>num;
    
    while(r-l>1e-8)
    {
        mid=(l+r)/2;
        if(num>mid*mid*mid) l=mid;
        else r=mid;
    }
    
    printf("%lf",l);
    
    return 0;
}

3.Acwing 730. 机器人跳跃问题

题目链接：730. 机器人跳跃问题 - AcWing题库

注意点：判断一个e是否可以跳到终点时，若e较大，大于每一栋楼的高度，那么之后e会呈指数级增长，会溢出变为负值，因此要增设跳出条件，若e大于等于最高楼高度，那么break；

AC代码：

#include<iostream>
using namespace std;

int n;
int a[100005]={0};
int l,r,mid,t=0;

int main()
{
    cin>>n;
    for(int i=1;i<=n;++i)
    {
        cin>>a[i];
        t=max(t,a[i]);
    }
    
    l=0,r=t;
    while(l<r)
    {
        mid=(l+r)>>1;
        int e=mid;
        bool yn=true;
        for(int i=0;i<n;++i)
        {
            if(e<a[i+1]) e-=a[i+1]-e;
            else e+=e-a[i+1];
            
            if(e<0) 
            {
                yn=false;
                break;
            }
            if(e>=t) break;//不判断e会成指数倍增长，会溢出 且该处必须大于等于最高楼，不可以是r，因为高度为r通过可能是因为在此前跳低楼获取能量
        }

        if(yn) r=mid;
        else l=mid+1;
    }
    
    cout<<l<<endl;
    
    return 0;
}

4.Acwing1227. 分巧克力

题目链接：1227. 分巧克力 - AcWing题库

#include<iostream>
using namespace std;

const int MAXN= 100005;

int n,k;
int a[MAXN],b[MAXN];

int main()
{
    cin>>n>>k;
    for(int i=1;i<=n;++i) cin>>a[i]>>b[i];
    int l=1,r=100000,mid=0;
    while(l<r)
    {
        mid=(l+r)/2;
        
        //计算分成的块数
        int num=0;
        for(int i=1;i<=n;++i)
        {
            num+=(a[i]/mid)*(b[i]/mid);
        }
        if(num>=k) l=mid+1;
        else r=mid;
    }
    
    cout<<r-1<<endl;
    
    return 0;
}