思路:
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][j−1]−sum[i−1][j]+sum[i−1][j−1]
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;
}