FNOI_TEST_1 c

本文介绍了一种解决区间子序列计数问题的有效算法。针对给定数列和区间限制,快速计算出满足条件的子区间数量。通过预处理和二分查找优化,实现了O(n+mlogn)的时间复杂度。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

http://www.elijahqi.win/archives/603
3.区间(c.cpp/c/pas)
【问题描述】 给你一个 n 个数的数列 A,其中每个数都是正整数。 给出一个正整数 k。 有 m 个询问,每个询问给出一个区间[L,R],询问有多少个区间[l,r],满足: 1.L<=l<=r<=R,即[l,r]是[L,R]的子区间 2.对于数列 A,区间[l,r]中所有数的和<=k 【输入数据】 第一行三个正整数 n、m、k。 第二行 n 个正整数表示数列 A。 接下来 m 行,第 i 行给出两个正整数 L、R,表示第 i 个询问。 【输出数据】 输出 m 行,第 i 行表示第 i 个询问的答案。 【输入输出样例】 c.in c.out 536 12345 15 12 34 8 3 2 【输入输出样例说明】 对于第 1 个询问:合法区间为[1,1][2,2][3,3][4,4][5,5][1,2][2,3][1,3],共 8 个 对于第 2 个询问:合法区间为[1,1][2,2][1,2],共 3 个 对于第 3 个询问:合法区间为[3,3][4,4],共 2 个 【数据规模和约定】 测试点编号 n m 特殊性质 1 <=10 <=10 2 <=100 <=100 3 4 <=1000 <=1000 A 中所有数均为 1 5 6
<=100000
=1 A 中所有数均为 1
7 8
<=100000
A 中所有数均为 1
9 10 对于 100%的数据: 1<=n<=100000, 1<=m<=100000, 1<=k<=10^9, 1<=L<=R<=n, 1<=A 中的数<=10000

还有o(m+n)的算法,还有莫队,明早刚

数据真的超级大啊,(long long)括号加错地方只有%50的得分

由于 p[i]递增,对于每个询问,前面一段 l 对应的 r 的最大值=p[l],后面一段 l 对 应的 r 的最大值为 R。二分出两段的分界点 x: ①L<=l<=x 的部分对答案的贡献为   x i 1)i-( p[i] L ,对 p[i]-i+1 维护前缀和即可 O(1)计算。 ②x+1<=l<=R 的部分对答案的贡献为    R x iR 1i )1( ,问题转化为等差数列求和, O(1)计算。 时间复杂度 O(n+mlogn)。

#include<cstdio>
#include<algorithm>
#define N 110000
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x;
}
int b[N],a[N],n,m,k;
long long s[N];
int main(){
    freopen("c.in","r",stdin);
    n=read();m=read();k=read();
    for (int i=1;i<=n;++i) a[i]=read();
    int sum=0,p=1;//p为指针 为在a数组中的指针 
    for (int i=1;i<=n;++i){
        if (a[i]>k){
            for (int z=p;z<=i;++z) b[z]=i-1;// i前面的都只可以延伸到i-1;其他情况在下一层已经处理 
            sum=0;p=i+1;
            continue;
        }
        while (sum+a[i]>k) sum-=a[p],b[p++]=i-1;sum+=a[i];
    }
    while (p<=n) b[p++]=n;//处理剩余一段的 
    //利用s处理B的前缀和
    for (int i=1;i<=n;++i) s[i]=b[i]+s[i-1];
    for (int i=1;i<=m;++i){
        int l=read(),r=read();
        int  x=std::lower_bound(b+l,b+r+1,r)-b-1;//lowerbound 二分查找(l,r)区间内值比r大or等于的地方,返回指针 
        long long ans=s[x]-s[l-1]-(long long)(l+x)*(x-l+1)/2+x-l+1;
        ans+=(long long)r*(r-x)-(long long)(x+1+r)*(r-x)/2+r-x;
        printf("%lld\n",ans);
    } 
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值