Codeforces 1037F. Maximum Reduction

总感觉我这种做法会T，一直没写，看了其他人的题解也是这样，，，就果断写了，，可能数据不太深，或者玄学复杂度

题意即求xk-1长度的所有区间的最大值的和，对每一个i(数组下边)，他对答案的贡献数量就是在以ar[i]为最大值的最大子区间中所有符合条件的区间数量

求ar[i]的作用区间，即，求最小的l，对于x>=l&&x<i，ar[x]<=ar[i]；求最大的r，对于k<=r&&k>i，ar[k]<ar[i]；则ar[i]的作用区间为[l,r]；我代码中的ar[i]的作用区间为（pre[i],nex[i]）

然后遍历i左右中更小的区间，求ar[i]对答案的贡献数量即可，代码思路应该很清晰

#include<iostream>
#include<cstdio>
#include<cmath>
#include<queue>
#include<vector>
#include<string.h>
#include<cstring>
#include<algorithm>
#include<set>
#include<map>
#include<fstream>
#include<cstdlib>
#include<ctime>
#include<list>
#include<climits>
#include<bitset>
using namespace std;
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define fopen freopen("input.in", "r", stdin);freopen("output.in", "w", stdout);
#define left asfdasdasdfasdfsdfasfsdfasfdas1
#define tan asfdasdasdfasdfasfdfasfsdfasfdas
typedef long long ll;
typedef unsigned int un;
const int desll[][]={{,},{,-},{,},{-,}};
const ll mod=1e9+;
const int maxn=2e6+;
const int maxm=1e9+;
const double eps=1e-;
int m,n,k;
int ar[maxn];
int pre[maxn],nex[maxn];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++)scanf("%d",&ar[i]);
    pre[]=;
    for(int i=;i<=n;i++){
        int ins=i-;
        while(ins> && ar[ins]<ar[i])ins=pre[ins];
        pre[i]=ins;
    }
    nex[n]=n+;
    for(int i=n-;i>;i--){
        int ins=i+;
        while(ins<=n && ar[ins]<=ar[i])ins=nex[ins];
        nex[i]=ins;
    }
    //for(int i=1;i<=n;i++)cout<<pre[i]<<" ";cout<<endl;
    //for(int i=1;i<=n;i++)cout<<nex[i]<<" ";cout<<endl;
    ll ans=;
    for(int i=;i<=n;i++){
        int l=pre[i],r=nex[i];
        ll a,b;
        if(i-l<r-i){
            for(int x=l+;x<=i;x++){//x为符合条件的区间的起始点
                a=(i-x+-+m-)/(m-);//a为包含i的区间中最少的（m-1）段数
                a=max(a,1LL);
                b=(r-x-)/(m-);//b为以r-1为区间终点，可以贡献的最多(m-1)段数
                ///cout<<"x = "<<x<<" "<<a<<" "<<b<<endl;
                if(b>=a){
                    ans = (ans + (b-a+)*ar[i]%mod)%mod;
                }
            }
        }
        else{
            for(int x=i;x<r;x++){
                a=(x-i+-+m-)/(m-);
                a=max(a,1LL);
                b=(x-l-)/(m-);
                //cout<<"x = "<<x<<" "<<a<<" "<<b<<endl;
                if(b>=a){
                    ans = (ans + (b-a+)*ar[i]%mod)%mod;
                }
            }
        }
        //cout<<i<<" "<<ans<<" "<<a<<" "<<b<<endl;
    }
    printf("%I64d\n",ans);

    return ;
}