BZOJ5055: 膜*

n<=300000个数求所有i<j<k,Ai>Aj>Ak的Ai*Aj*Ak的和。

还以为什么奇怪题呢。。就是个偏序。。

i<j,Ai>Aj就是逆序对模板，把树状数组里的+1改成+Ai即可，每个数的Ai*Aj数就树状数组求个前缀和即可。

上面那个算出来存在Bi，再用类似的方法求一次即可。

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
//#include<queue>
//#include<iostream>
using namespace std;

int n;
#define maxn 300011
const int mod=19260817;
int a[maxn],id[maxn],two[maxn],thr[maxn];
int lisan[maxn];
struct BIT
{
    int a[maxn];
    void clear() {memset(a,0,sizeof(a));}
    void add(int x,int v) {for (;x<=n;x+=x&-x) a[x]=(a[x]+v)%mod;}
    int query(int x) {int ans=0;for (;x;x-=x&-x) ans=(ans+a[x])%mod;return ans;}
}t;
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]),lisan[i]=a[i];
    sort(lisan+1,lisan+1+n);
    for (int i=1;i<=n;i++) id[i]=lower_bound(lisan+1,lisan+1+n,a[i])-lisan;
    t.clear();
    for (int i=1;i<=n;i++)
    {
        t.add(id[i],a[i]);
        two[i]=1ll*a[i]*t.query(id[i]-1)%mod;
    }
    t.clear();
    int ans=0;
    for (int i=1;i<=n;i++)
    {
        t.add(id[i],two[i]);
        ans=(ans+1ll*a[i]*t.query(id[i]-1)%mod)%mod;
    }
    printf("%d
",ans);
    return 0;
}

跑得有点慢 难道被xu掉了？