BZOJ 3745

分治.....

之前就了解过这种分治统计答案的算法,对于当前的区间[l,r],我们考虑过中间的那条线的区间,这种题往往都存在单调性,我们发现min和max都是随位置单调的,我们枚举左端点x,然后维护两个指针p1,p2,表示[mid+1,p1/p2]这个区间的最值大于/小于[x,mid]的最值的最远的p1/p2,那么答案就可以分3段统计,一段是[mid+1,min(p1,p2)],右端点在这一部分的区间的min和max都一样,比较容易求出。另一段是[max(p1,p2)+1,r],在这一部分的区间min和max为[mid+1,y]的最值,考虑答案的式子可以拆成:∑min*max*r-∑min*max*(l-1),那么我们就可以预处理出min*max*r和min*max的前缀和,然后就可以算出这一部分对答案的贡献。最后一段是[min(p1,p2)+1,max(p1,p2)],这个对答案的贡献和上一个和像,维护2个前缀和即可。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 using namespace std;
 5 #define maxn 500005
 6 typedef long long LL;
 7 const int mo=1000000000;
 8 int n,a[maxn];
 9 LL ans,sum1[maxn][2],sum2[maxn][2],sum3[maxn][2];
10
11 inline int read(void)
12 {
13     int x=0;
14     char ch=getchar();
15     while (ch>‘9‘||ch<‘0‘) ch=getchar();
16     while (ch>=‘0‘&&ch<=‘9‘)
17     {
18         x=x*10+ch-‘0‘;
19         ch=getchar();
20     }
21     return x;
22 }
23
24 void solve(int l,int r)
25 {
26     if (l==r)
27     {
28         ans=(ans+(LL)a[l]*a[l])%mo;
29         return;
30     }
31     int mid=(l+r)>>1;
32     int mn=1e9,mx=-1;
33     sum1[mid][0]=sum1[mid][1]=0;
34     sum2[mid][0]=sum2[mid][1]=0;
35     sum3[mid][0]=sum3[mid][1]=0;
36     for (int i=mid+1;i<=r;i++)
37     {
38         mn=min(mn,a[i]);
39         mx=max(mx,a[i]);
40         sum1[i][0]=(LL)mn*mx%mo*i%mo;sum1[i][1]=(LL)mn*mx%mo;
41         sum2[i][0]=(LL)mn*i%mo;sum2[i][1]=mn;
42         sum3[i][0]=(LL)mx*i%mo;sum3[i][1]=mx;
43     }
44     for (int i=mid+2;i<=r;i++)
45     {
46         (sum1[i][0]+=sum1[i-1][0])%=mo;(sum1[i][1]+=sum1[i-1][1])%=mo;
47         (sum2[i][0]+=sum2[i-1][0])%=mo;(sum2[i][1]+=sum2[i-1][1])%=mo;
48         (sum3[i][0]+=sum3[i-1][0])%=mo;(sum3[i][1]+=sum3[i-1][1])%=mo;
49     }
50     mn=1e9,mx=-1;
51     int t1=a[mid+1],t2=a[mid+1];
52     int p1=mid,p2=mid;
53     for (int i=mid;i>=l;i--)
54     {
55         mn=min(mn,a[i]);mx=max(mx,a[i]);
56         while (t1>=mn&&p1<r) p1++,t1=min(t1,a[p1+1]);
57         while (t2<=mx&&p2<r) p2++,t2=max(t2,a[p2+1]);
58         int tmp=min(p1,p2);
59         LL w=(LL)(mid+2-i+tmp-i+1)*(tmp-mid)/2;w%=mo;
60         ans=(ans+w*mn%mo*mx%mo)%mo;
61         tmp=max(p1,p2)+1;
62         w=(sum1[r][0]-sum1[tmp-1][0]+mo)%mo;
63         w=(w-(sum1[r][1]-sum1[tmp-1][1]+mo)%mo*(i-1)%mo+mo)%mo;
64         ans=(ans+w)%mo;
65         if (p1<=p2)
66         {
67             w=(sum2[p2][0]-sum2[p1][0]+mo)%mo*mx%mo;
68             w=(w-(sum2[p2][1]-sum2[p1][1]+mo)%mo*(i-1)%mo*mx%mo+mo)%mo;
69             ans=(ans+w)%mo;
70         }
71         else
72         {
73             w=(sum3[p1][0]-sum3[p2][0]+mo)%mo*mn%mo;
74             w=(w-(sum3[p1][1]-sum3[p2][1]+mo)%mo*(i-1)%mo*mn%mo+mo)%mo;
75             ans=(ans+w)%mo;
76         }
77     }
78     solve(l,mid);
79     solve(mid+1,r);
80 }
81
82 int main()
83 {
84     n=read();
85     for (int i=1;i<=n;i++) a[i]=read();
86     solve(1,n);
87     printf("%lld\n",ans);
88     return 0;
89 }
时间: 2024-09-28 16:13:53

BZOJ 3745的相关文章

bzoj 3745 [Coci2015]Norma——序列分治

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3745 如果分治,就能在本层仅算过 mid 的区间了. 可以从中间到左边地遍历左边,给右边两个指针,表示第一个更新左边造成的最小值/最大值的位置. 两个位置共同的左边可以公式算长度,用左边的最值算:两个位置共同的右边可以预处理,处理出 算上长度(相对mid的)的最值乘积求和 与 不算长度的最值乘积求和(都是前缀),把前者加到答案里,后者乘上左边到mid的长度加到答案里即可:两个位置夹着的位置

bzoj 3745: [Coci2015]Norma

Description 给定序列\(a_i\) 求 \[\sum_{i=1}^n \sum_{j=i}^n (j-i+1)\max\{a_i,a_{i+1}\cdots a_j\}\min\{a_i,a_{i+1}\cdots a_j\}\] Input 第1行,一个整数N: 第2~n+1行,每行一个整数表示序列a. Output 输出答案对10^9取模后的结果. Sample Input 4 2 4 1 4 Sample Output 109 [数据范围] \(N \le 500000\) \

bzoj 3745: [Coci2015]Norma【分治】

参考:https://blog.csdn.net/lych_cys/article/details/51203960 真的不擅长这种-- 分治,对于一个(l,r),先递归求出(l,mid),(mid+1,r),然后这个区间对答案贡献的就是经过mid的区间 我们先预处理出mid为l的右端点的mx*mn*len的前缀和与mx*mn的前缀和,然后枚举左端点,右端点维护两个下标j,k,分别表示mn和mx在左端点时的合法右端点 然后分三种情况处理,假设j<k 1.右端点在(mid+1,j)时,直接计算 2

BZOJ 1013: [JSOI2008]球形空间产生器sphere

二次联通门 : BZOJ 1013: [JSOI2008]球形空间产生器sphere /* BZOJ 1013: [JSOI2008]球形空间产生器sphere 高斯消元 QAQ SB的我也能终于能秒题了啊 设球心的坐标为(x,y,z...) 那么就可以列n+1个方程,化化式子高斯消元即可 */ #include <cstdio> #include <iostream> #include <cstring> #define rg register #define Max

bzoj 3309 DZY Loves Math - 莫比乌斯反演 - 线性筛

对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0. 给定正整数a,b,求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b). Input 第一行一个数T,表示询问数. 接下来T行,每行两个数a,b,表示一个询问. Output 对于每一个询问,输出一行一个非负整数作为回答. Sample Input 4 7558588 9653114 6514903 445

【BZOJ】[HNOI2009]有趣的数列

[算法]Catalan数 [题解] 学了卡特兰数就会啦>_<! 因为奇偶各自递增,所以确定了奇偶各自的数字后排列唯一. 那么就是给2n个数分奇偶了,是不是有点像入栈出栈序呢. 将做偶数标为-1,做奇数标为+1,显然当偶数多于奇数时不合法,因为它压不住后面的奇数. 然后其实这种题目,打表就可知啦--QAQ 然后问题就是求1/(n+1)*C(2n,n)%p了,p不一定是素数. 参考bzoj礼物的解法. 看到网上清一色的素数筛+分解质因数解法,不解了好久,感觉写了假的礼物-- 后来觉得礼物的做法才比

洛谷 P2709 BZOJ 3781 小B的询问

题目描述 小B有一个序列,包含N个1~K之间的整数.他一共有M个询问,每个询问给定一个区间[L..R],求Sigma(c(i)^2)的值,其中i的值从1到K,其中c(i)表示数字i在[L..R]中的重复次数.小B请你帮助他回答询问. 输入输出格式 输入格式: 第一行,三个整数N.M.K. 第二行,N个整数,表示小B的序列. 接下来的M行,每行两个整数L.R. 输出格式: M行,每行一个整数,其中第i行的整数表示第i个询问的答案. 输入输出样例 输入样例#1: 6 4 3 1 3 2 1 1 3

BZOJ 1012: [JSOI2008]最大数maxnumber(线段树)

012: [JSOI2008]最大数maxnumber Time Limit: 3 Sec  Memory Limit: 162 MB Description 现在请求你维护一个数列,要求提供以下两种操作:1. 查询操作.语法:Q L 功能:查询当前数列中末尾L个数中的最大的数,并输出这个数的值.限制:L不超过当前数列的长度.2. 插入操作.语法:A n 功能:将n加上t,其中t是最近一次查询操作的答案(如果还未执行过查询操作,则t=0),并将所得结果对一个固定的常数D取模,将所得答案插入到数列

【BZOJ】【1016】【JSOI2008】最小生成树计数

Kruskal/并查集+枚举 唉我还是too naive,orz Hzwer 一开始我是想:最小生成树删掉一条边,再加上一条边仍是最小生成树,那么这两条边权值必须相等,但我也可以去掉两条权值为1和3的,再加上权值为2和2的,不也满足题意吗?事实上,如果这样的话……最小生成树应该是1和2,而不是1和3或2和2!!! 所以呢?所以对于一个图来说,最小生成树有几条边权为多少的边,都是固定的!所以我们可以做一遍Kruskal找出这些边权,以及每种边权出现的次数.然后,对于每种边权,比方说出现了$v_i$