题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6058
题外话:这场多校,真心感觉自己的无力,全队都发挥的很差,结束的时候排名掉到了90多,后期没做出字典树那个题,直到现在看到标程也依然不懂那个题。。。状态非常不好,力求可以好好调整一下。
题意:题目给了一个求和式子,求所有l,r区间的第K大的和。K<=80。
解法:
全队,开始讨论如何维护这K个比他大的值,然后队友去写主席树,我实在没想到这题的做法,之后队友想到了直接用双链表来求出对于一个数X左边最近的kk个比他大的和右边最近k个比他大的,扫一下就可以知道有几个区间的k大值是X,然后我们从小到大枚举X,每次维护一个链表,链表里面只有大于等于X的数,然后从左往右扫一遍就可以了。链表的删除是O(1)的。
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
struct FastIO
{
static const int S = 1310720;
int wpos;
char wbuf[S];
FastIO() : wpos(0) {}
inline int xchar()
{
static char buf[S];
static int len = 0, pos = 0;
if (pos == len)
pos = 0, len = fread(buf, 1, S, stdin);
if (pos == len) return -1;
return buf[pos ++];
}
inline int xuint()
{
int c = xchar(), x = 0;
while (c <= 32) c = xchar();
for (; '0' <= c && c <= '9'; c = xchar()) x = x * 10 + c - '0';
return x;
}
inline int xint()
{
int s = 1, c = xchar(), x = 0;
while (c <= 32) c = xchar();
if (c == '-') s = -1, c = xchar();
for (; '0' <= c && c <= '9'; c = xchar()) x = x * 10 + c - '0';
return x * s;
}
inline void xstring(char *s)
{
int c = xchar();
while (c <= 32) c = xchar();
for (; c > 32; c = xchar()) * s++ = c;
*s = 0;
}
inline void wchar(int x)
{
if (wpos == S) fwrite(wbuf, 1, S, stdout), wpos = 0;
wbuf[wpos ++] = x;
}
inline void wint(LL x)
{
if (x < 0) wchar('-'), x = -x;
char s[24];
int n = 0;
while (x || !n) s[n ++] = '0' + x % 10, x /= 10;
while (n--) wchar(s[n]);
}
inline void wstring(const char *s)
{
while (*s) wchar(*s++);
}
~FastIO()
{
if (wpos) fwrite(wbuf, 1, wpos, stdout), wpos = 0;
}
} io;
const int maxn = 5e5+10;
int T,n,k,pre[maxn],nxt[maxn],pos[maxn];
inline void del(int i){
pre[nxt[i]] = pre[i];
nxt[pre[i]] = nxt[i];
}
inline LL solve(int x){
static int a[82], b[82];
int ca=0,cb=0;
for(int i=x;i;i=pre[i]){
a[++ca]=i-pre[i];
if(ca==k) break;
}
for(int i=x;i<=n;i=nxt[i]){
b[++cb]=nxt[i]-i;
if(cb==k) break;
}
LL ret=0;
for(int i=1; i<=ca; i++)if(k-i+1<=cb) ret+=(LL)a[i]*b[k-i+1];
return ret;
}
int main()
{
T = io.xint();
while(T--){
n = io.xint();
k = io.xint();
for(int i=1; i<=n; i++){
int x;
x = io.xint();
pos[x]=i;
}
for(int i=0; i<=n+1; i++) pre[i]=i-1,nxt[i]=i+1;
pre[0]=0;
nxt[n+1]=n+1;
LL ans=0;
for(int i=1; i<=n; i++){
int x = pos[i];
ans += (LL)solve(x)*i;
del(x);
}
printf("%lld\n", ans);
}
return 0;
}