[HNOI2017]影魔

Posted D O Time

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[HNOI2017]影魔相关的知识,希望对你有一定的参考价值。

标签:线段树

题解:

  首先对于题目的条件进行分析,p1的条件是一个区间的两个端点必须是这一个区间的最大与次大值。p2的条件是一个区间的一个端点是最大值,而另一个端点不是次大值。显然,他们都需要两个条件(对于两个端点都有要求,这样的话不太好操作)。感受一下,p1的条件更为严苛,于是我们这样考虑:
  一个区间的一个端点是最大值(两个条件的公共部分),另一个不管,此时加上p2的贡献。
  我们可以通过线段树求解出这一个东西。对于i,找到右边大于他的第一个数,如果没有,自然是最后一个数,但是这样不好,于是我们加一个第n+1个数为+∞。然后设ri为右边大于他的第一个数,那么[i,i+1]、[i,i+2]、......[i,ri]都符合条件,于是在线段树中对于[i+1,ri]都+=p2。对于全部询问排序,给以i为左端点的区间加上贡献,查询[i+1,R]的和。反过来再搞一遍,注意询问区间也要反过来。
  然后我们发现题目中的对另一个端点的限制是互补的,我们的区间[i,ri]满足p1,而[i,i+1~ri-1]满足p2。因为只有一个满足p1,而且一个i对应一个ri,于是我们在ri处+p1-p2-p2即可,完美解决这一个问题。
小结:此题关键在于两个条件有一半相同,而另一半互补,所以可以不管互补的那一半到时候再考虑。

  1 #include<cstdio>
  2 #include<cstring>
  3 #include<iostream>
  4 #include<algorithm>
  5 #define ls k*2
  6 #define rs (k*2+1)
  7 #define LL long long
  8 using namespace std;
  9 const int MAXN=210000;
 10 int n,m,p1,p2,tp;
 11 int v[MAXN],aft[MAXN],st[MAXN];
 12 LL ans[MAXN],sum[MAXN*5],lz[MAXN*5];
 13 struct ed
 14 {
 15   int L,R,id;
 16 }ask[MAXN];
 17 #define down() 18 { 19   sum[ls]+=lz[k]*(mid-ll+1); 20   sum[rs]+=lz[k]*(rr-mid); 21   lz[ls]+=lz[k]; 22   lz[rs]+=lz[k]; 23   lz[k]=0; 24 }
 25 inline int gi() {int res; scanf("%d",&res); return res;}
 26 bool comp(ed x,ed y){return x.L<y.L;}
 27 int find(int x)
 28 {
 29   int L=1,R=m,mid;
 30   while(L<=R)
 31     {
 32       mid=(L+R)/2;
 33       if(ask[mid].L<x)
 34         L=mid+1;
 35       else
 36         R=mid-1;
 37     }
 38   return L;
 39 }
 40 void add(int k,int ll,int rr,int L,int R,LL Val)
 41 {
 42   if(ll==L && rr==R)
 43     {
 44       sum[k]+=Val*(rr-ll+1);
 45       lz[k]+=Val;
 46       return;
 47     }
 48   int mid=(ll+rr)/2;
 49   down();
 50   if(R<=mid)
 51     add(ls,ll,mid,L,R,Val);
 52   else if(mid<L)
 53     add(rs,mid+1,rr,L,R,Val);
 54   else
 55     {
 56       add(ls,ll,mid,L,mid,Val);
 57       add(rs,mid+1,rr,mid+1,R,Val);
 58     }
 59   sum[k]=sum[ls]+sum[rs];
 60 }
 61 LL query(int k,int ll,int rr,int L,int R)
 62 {
 63   if(ll==L && rr==R) return sum[k];
 64   int mid=(ll+rr)/2;
 65   down();
 66   if(R<=mid)
 67     return query(ls,ll,mid,L,R);
 68   else if(mid<L)
 69     return query(rs,mid+1,rr,L,R);
 70   else
 71     return query(ls,ll,mid,L,mid)+query(rs,mid+1,rr,mid+1,R);
 72 }
 73 void work()
 74 {
 75   tp=0; st[0]=n+1;
 76   for(int i=n;i>=1;i--)
 77     {
 78       while(tp && v[st[tp]]<v[i]) tp--;
 79       aft[i]=st[tp];
 80       st[++tp]=i;
 81     }
 82   memset(sum,0,sizeof sum);
 83   memset(lz,0,sizeof lz);
 84   for(int i=n;i>=1;i--)
 85     {
 86       if(aft[i]>i)
 87         {
 88           add(1,1,n+1,i+1,aft[i],p2);
 89           add(1,1,n+1,aft[i],aft[i],p1-2*p2);
 90         }
 91       int p=find(i);
 92       while(ask[p].L==i)
 93         {
 94           ans[ask[p].id]+=query(1,1,n+1,i+1,ask[p].R);
 95           p++;
 96         }
 97     }
 98 }
 99 int main()
100 {
101   n=gi();m=gi();p1=gi();p2=gi();
102   for(int i=1;i<=n;i++) v[i]=gi();
103   for(int i=1;i<=m;i++) ask[i].L=gi() , ask[i].R=gi() , ask[i].id=i;
104   sort(ask+1,ask+1+m,comp);
105   work();
106   reverse(v+1,v+1+n);
107   for(int i=1;i<=m;i++)
108     {
109       swap(ask[i].L,ask[i].R);
110       ask[i].L=n-ask[i].L+1;
111       ask[i].R=n-ask[i].R+1;
112     }
113   sort(ask+1,ask+1+m,comp);
114   work();
115   for(int i=1;i<=m;i++)
116     printf("%lld\n",ans[i]);
117   return 0;
118 }

以上是关于[HNOI2017]影魔的主要内容,如果未能解决你的问题,请参考以下文章

[BZOJ4826][HNOI2017]影魔(主席树)

[HNOI2017]影魔

[HNOI2017]影魔

[AH/HNOI2017]影魔

bzoj4826 [Hnoi2017]影魔

HNOI2017影魔