HDU 6058 枚举 Kanade's sum

题目

Kanade's sum

Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 2496    Accepted Submission(s): 1037

Problem Description

Give you an array A[1..n]  of length n. 

Let f(l,r,k) be the k-th largest element of A[l..r].

Specially , f(l,r,k)=0 if r−l+1<k.

Give you k , you need to calculate ∑nl=1∑nr=lf(l,r,k)

There are T test cases.

1≤T≤10

k≤min(n,80)

A[1..n] is a permutation of [1..n]

∑n≤5∗105

 

Input

There is only one integer T on first line.

For each test case,there are only two integers n,k on first line,and the second line consists of n integers which means the array A[1..n]

 

Output

For each test case,output an integer, which means the answer.

 

Sample Input

15 21 2 3 4 5

 

Sample Output

题目大意

   给我们一个数组A，给了我们一个区间，变成for循环就是for(l从1到n) for(r从l到n) 求区间内第K大的数，把这些第K大的数累计求和

解题思路

  采用枚举的思想，从每一个数组中的元素遍历他的左右边看分别有多少个比他大的数，这样我们就可以知道该元素可以有多少个第K大的组合可能

#include<iostream>#include<cstdio>using namespace std;#define read(a) scanf("%d",&a)#define LL long longconst int maxn=500000+10;int a[maxn];int l[maxn],r[maxn];int main(){       //freopen("1003.in", "r", stdin);       //freopen("data.out", "w", stdout);    int t;    read(t);    while(t--)    {        int n,k;        read(n);        read(k);        for(int i=0;i<n;i++)        {            read(a[i]);        }        LL ans=0;        for(int i=0;i<n;i++)        {            int lcnt=1,rcnt=1,j;//lcnt代表元素a[i]左边比他大的数有多少个，rcnt同理            for( j=i+1;j<n;j++)            {                if(rcnt>k)                    break;                if(a[j]>a[i])                {                    r[rcnt++]=j-i;//r[rcnt]代表右边第rcnt个比a[i]大的数距离a[i]的距离，这个是方便计算的，可以等于j,                                  //但是计算的时候要特殊处理右边只有一个比a[i]大的时候，下方的rnum=1，比较麻烦,                                  //原来是那样做的，不建议                }            }            if(j>=n)                r[rcnt]=n-i; //如果a[i]右边比他大的数没超过k个，                             //那么我们知道a[i]右边比他大的数只有rcnt-1个，                             //我们假设距离a[i]最远的比他大的那个数为righht，                             //（程序中没有right这个变量，这里就是为了方便理解）                             //这里的r[rcnt]就是为了方便后面统计right右边有多少个比a[i]小的数            for(j=i-1;j>=0;j--)            {                if(lcnt>k)                    break;                if(a[j]>a[i])                {                    l[lcnt++]=i-j;//同理上面                }            }            if(j<=0)                l[lcnt]=i+1;//同理上面            for(j=0;j<lcnt;j++)            {                if(k-j-1>=rcnt)                    continue;                int lnum=l[j+1]-l[j];                int rnum=r[k-j]-r[k-j-1];                ans+=(LL)a[i]*lnum*rnum;            }        }        printf("%lld\n",ans);    }    return 0;}

链表的代码，用链表省去了每次去寻找和比较比他大的数的时间，随着N的增大节省时间越多

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;#define LL long long#define read(a) scanf("%d",&a)const int maxn=5*1e5+10;int pre[maxn],nex[maxn],pos[maxn],f[maxn],b[maxn],a[maxn],n,k;void dele(int p)//更新前驱和后继{    if(p==1)        pre[nex[p]]=pre[p];    else if(p==n)        nex[pre[p]]=nex[p];    else    {        pre[nex[p]]=pre[p];        nex[pre[p]]=nex[p];    }    pre[p]=nex[p]=0;}void work(){  //  int n,k;    memset(f,0,sizeof(f));    memset(b,0,sizeof(b));    read(n);read(k);    for(int i=1;i<=n;i++)    {        read(a[i]);        pos[a[i]]=i;//记录a[i]的位置        pre[i]=i-1;//记录前驱        nex[i]=i+1;//记录后继    }    LL ans=0;    for(int i=1;i<=n-k+1;i++)//本代码的核心所在，因为是1-n的数所以可以是第k大的只可能是1-(n-k+1)这几个数                             //从1开始找，因为1最小所以它的两边都是比他大的数，找完1之后，在dele()函数里面把1从链表里面拿下来                             //这样继续找2，因为1不在链表里面了，所以2的两边都是比它大的数，以此类推直到n-k+1    {        int fro=0,beh=0;//front behind        int p=pos[i];        for(int v=p;v&&fro<=k;v=pre[v]) f[fro++]=v;        for(int v=p;v<=n&&beh<=k;v=nex[v]) b[beh++]=v;        f[fro]=0,b[beh]=n+1;        for(int j=0;j<fro&&j<k;j++)//fro有可能是k+1,要有j<k限制他，MMP，找了半个多小时，烦躁        {            if(k-1-j>beh-1)                continue;           {                int lnum=f[j]-f[j+1];                int rnum=b[k-1-j+1]-b[k-1-j];                ans+=1ll*i*lnum*rnum;           }        }        dele(p);    }    printf("%lld\n",ans);}int main(){      //freopen("1003.in", "r", stdin);      //freopen("data.out", "w", stdout);    int T;    read(T);    while(T--)    {        work();    }    return 0;}