2021-08-09校网比赛B题

在这里插入图片描述
在这里插入图片描述
思路:

30pts:
预处理出二维前缀和,循环枚举左上端点和右下端点,
判断区间和是否在闲置内即可。
区间和求值,设sum[i][j]表示以(i,j)为右下端点的二维前缀和
d i s = s u m [ i ] [ j ] − s u m [ i ] [ j − 1 ] − s u m [ i − 1 ] [ j ] + s u m [ i − 1 ] [ j − 1 ] dis = sum[i][j]-sum[i][j-1] - sum[i-1][j]+sum[i-1][j-1] dis=sum[i][j]sum[i][j1]sum[i1][j]+sum[i1][j1]

100pts:
对于一个区间和,设左上端点为(x1,y1),右下端点为(x2,y2)
它的权值计算为:
sum=
                 \;\;\;\;\;\;\;\; a[x1]*b[y1]+a[x1]*b[y1+1]+a[x1]*b[y1+…]+a[x1]*b[y2]      + \;\;+ +
                 \;\;\;\;\;\;\;\; a[x1+1]*b[y1]+a[x1+1]*b[y1+1]+a[x1+1]*b[y1+…]+a[x1+1]*b[y2]      + \;\;+ +
                 \;\;\;\;\;\;\;\; a[x1+2]*b[y1]+a[x1+2]*b[y1+1]+a[x1+2]*b[y1+…]+a[x1+1]*b[y2]      + \;\;+ +
                 \;\;\;\;\;\;\;\; a[x1+…]*b[y1]+a[x1+…]*b[y1+1]+a[x1+…]*b[y1+…]+a[x1+…]*b[y2]    + \;+ +
                 \;\;\;\;\;\;\;\; a[x2]*b[y1]+a[x2]*b[y1+1]+a[x2]*b[y1+…]+a[x2]*b[y2]

提取一下公因式:
( a [ x 1 ] + a [ x 1 + 1 ] + a [ x 1 + 2 ] + . . . + a [ x 2 ] ) ∗ ( b [ y 1 ] + b [ y 1 + 1 ] + b [ y 1 + 2 ] + . . . + b [ y 2 ] ) (a[x1]+a[x1+1]+a[x1+2]+...+a[x2])*(b[y1]+b[y1+1]+b[y1+2]+...+b[y2]) (a[x1]+a[x1+1]+a[x1+2]+...+a[x2])(b[y1]+b[y1+1]+b[y1+2]+...+b[y2])

发现x跟y没有必然的关系,压根不冲突。
考虑分开去求。
预处理出x和y的前缀和suma和sumb
枚举a所有可能的区间(x1,x2)
对于限制L,R来说,
设tmp=suma[x2]-suma[x1-1]
b的限制则为(L/tmp)=l,(R/tmp)=r
考虑建立一个数组c,表示对于y来说的某一段累加和。
因为c可以具有单调性(将c进行升序排序),我们用二分来求y的两个点。
那我们用二分求出第一个大于等于l的c(设其为x),
并求出第一个大于r的c(设其为y).
那可供我们选择的区间数就为y-x.
最后用ans将二分出所有区间数加总,就是答案。

时间复杂度:O( n 2 l o g 2 ( m 2 ) n^2log2(m^2) n2log2(m2))

#include <cstdio>
#include <iostream>
#include <algorithm>

using namespace std;

typedef long long ll;

const int N = 1e3 + 10;
ll n, m, L, R, a[N], b[N], suma[N], sumb[N], c[N * N], ans, cnt; 

ll binarysearch(ll l, ll r)
{
	ll x = lower_bound(c + 1, c + 1 + cnt, l) - c;
	ll y = upper_bound(c + 1, c + 1 + cnt, r) - c;
	return y - x;
}
int main()
{
	scanf("%lld%lld%lld%lld", &n, &m, &L, &R);
	for(ll i = 1; i <= n; i++) scanf("%lld", &a[i]), suma[i] = suma[i - 1] + a[i];
	for(ll i = 1; i <= m; i++) scanf("%lld", &b[i]), sumb[i] = sumb[i - 1] + b[i];
	for(ll i = 1; i <= m; i++)
	 for(ll j = i; j <= m; j++)
	 	c[++cnt] = sumb[j] - sumb[i - 1];
	sort(c + 1, c + 1 + cnt);
	for(ll i = 1; i <= n; i++)
	 for(ll j = i; j <= n; j++)
		{
			ll tmp = suma[j] - suma[i - 1];
			ll l = L / tmp, r = R / tmp;
			if(l % tmp) ++l; 
			ans += binarysearch(l, r);
 		}
	printf("%lld",ans);
	return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值