zoukankan      html  css  js  c++  java
  • 【2019多校】Make Rounddog Happy(分治)

    Rounddog always has an array a1,a2,,ana1,a2,⋯,an in his right pocket, satisfying 1ain1≤ai≤n .

    A subarray is a non-empty subsegment of the original array. Rounddog defines a good subarray as a subsegment al,al+1,,aral,al+1,⋯,ar that all elements in it are different and max(al,al+1,,ar)(rl+1)kmax(al,al+1,…,ar)−(r−l+1)≤k .

    Rounddog is not happy today. As his best friend, you want to find all good subarrays of aa to make him happy. In this case, please calculate the total number of good subarrays of aa .

    InputThe input contains several test cases, and the first line contains a single integer T (1T20)T (1≤T≤20) , the number of test cases.

    The first line of each test case contains two integers n (1n300000)n (1≤n≤300000) and k (1k300000)k (1≤k≤300000) .

    The second line contains nn integers, the ii -th of which is ai (1ain)ai (1≤ai≤n) .

    It is guaranteed that the sum of nn over all test cases never exceeds 10000001000000 .
    OutputOne integer for each test case, representing the number of subarrays Rounddog likes.

    Sample Input

    2
    5 3
    2 3 2 2 5
    10 4
    1 5 4 3 6 2 10 8 4 5

    Sample Output

    7
    31
    SOLUTION:
    很容易想到考虑每一个数作为最大值的情况,然后处理出他所能想两边扩展的区间
    但是向两边的很好搞,但是必须跨区间,然后就自闭了
    后来知道分治可以很好的搞这个问题
    题解:https://www.cnblogs.com/Dillonh/p/11390927.html
    启发式的所有东西都是玄学的复杂度。。。。。qwq

    #include<bits/stdc++.h>
    #define ll long long
    #define rep(i,a,b) for(int i=a;i<=b;i++)
    using namespace std;
    const int maxn=1000010;
    int a[maxn],st[maxn][21],lg[maxn],K,N; ll ans;
    int A[maxn],B[maxn],vis[maxn];
    int get(int L,int R)
    {
        int k=lg[R-L+1];
        return a[st[L][k]]>=a[st[R-(1<<k)+1][k]]?st[L][k]:st[R-(1<<k)+1][k];
    }
    void solve(int L,int R)
    {
        if(L>R) return ;
        int pos=get(L,R);
        if(pos-L<R-pos){
            rep(i,L,pos){
                int t=a[pos]-K,lR=i+t-1;
                int fcy=min(R,B[i]);
                lR=max(lR,pos);
                if(lR>fcy) continue;
                ans+=fcy-lR+1;
            }
        }
        else {
            rep(i,pos,R){
                int t=a[pos]-K,rL=i-t+1;
                int fcy=max(L,A[i]);
                rL=min(rL,pos);
                if(rL<fcy) continue;
                ans+=rL-fcy+1;
            }
        }
        solve(L,pos-1); solve(pos+1,R);
    }
    int main()
    {
        int T;
        lg[0]=-1; rep(i,1,maxn-1) lg[i]=lg[i>>1]+1;
        scanf("%d",&T);
        while(T--){
            scanf("%d%d",&N,&K); ans=0;
            rep(i,1,N) scanf("%d",&a[i]);
            rep(i,1,N) st[i][0]=i;
            rep(i,1,20) {
                rep(j,1,N+1-(1<<i))
                 st[j][i]=a[st[j][i-1]]>=a[st[j+(1<<(i-1))][i-1]]?st[j][i-1]:st[j+(1<<(i-1))][i-1];
            }
        
            rep(i,1,N) vis[i]=0; A[1]=1; vis[a[1]]=1;
            rep(i,2,N){
                 if(vis[a[i]])  A[i]=max(A[i-1],vis[a[i]]+1);
                 else A[i]=A[i-1];
                 vis[a[i]]=i;
            }
            rep(i,1,N) vis[i]=0; B[N]=N; vis[a[N]]=N;
            for(int i=N-1;i>=1;i--){
                if(vis[a[i]]) B[i]=min(B[i+1],vis[a[i]]-1);
                else B[i]=B[i+1];
                vis[a[i]]=i;
            }
            solve(1,N);
            printf("%lld
    ",ans);
        }
        return 0;
    } 
    

      












  • 相关阅读:
    [原创]ExtAspNet秘密花园(二) — 一切从头开始
    使用领域定义建模语言来提高生产力
    报表控件团队成立了!!
    常见的程序语言分类的区别(简要说明)
    怎样做一名“专业”的程序员?
    推荐一个计数器网站
    在c中嵌入sql操作PostgreSQL
    让Visual C++ 2005 Express Edition也可以开发windows应用程序
    使用PowerDesigner进行代码生成
    sdmagazine与Dr. Dobb's Journal合并了?
  • 原文地址:https://www.cnblogs.com/zhangbuang/p/11604831.html
Copyright © 2011-2022 走看看