bzoj1145: [CTSC2008]图腾totem

最新推荐文章于 2020-01-13 20:23:43 发布

原创最新推荐文章于 2020-01-13 20:23:43 发布 · 454 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#树状数组 #BZOJ #OI #数学

bzoj 同时被 3 个专栏收录

16 篇文章

订阅专栏

数学

2 篇文章

订阅专栏

树状数组

1 篇文章

订阅专栏

本文介绍了一道使用树状数组解决的复杂算法题目，并详细展示了如何通过树状数组进行有效的数据处理和查询操作。文章提供了完整的代码实现，适合希望深入了解树状数组应用的读者。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

数学+树状数组

真是一道神题233果然符合CTSC

用来练习思维与树状数组很不错的

膜拜hzwer.com/5357.html

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
using namespace std;
int op[200005];
int a[200005];
int lc[200005],rc[200005];
int n;
int xa=0,xb=0,xc=0,xd=0;
int ans=0;
const int P=16777216;
inline int lowbit(int x)
{
       return x&(-x);
}
int find(int x)
{
    int sum=0;
    while(x>0){
         sum+=op[x];
         x-=lowbit(x);
         }
    return sum;
}
void add(int x,int val)
{
     while(x<=n){
          op[x]+=val;
          x+=lowbit(x);
          }
}
int main()
{
    int i,j,k,l;
    scanf("%d",&n);
    for(i=1;i<=n;i++)scanf("%d",&a[i]);
    for(i=1;i<=n;i++){
        lc[i]=find(a[i]);
        add(a[i],1);
        }
    for(i=1;i<=n;i++)
        rc[i]=a[i]-1-lc[i];
    memset(op,0,sizeof op);
    for(i=1;i<=n;i++){
        ans=(ans+find(a[i])*(n-i-rc[i]))%P;
        add(a[i],lc[i]);
        }
    memset(op,0,sizeof op);
    for(i=1;i<=n;i++){
       ans=(ans+(lc[i]*(i-1)-find(a[i])-lc[i]*(lc[i]-1)/2)*(n-i-rc[i])%P)%P;
       add(a[i],i);
       }
    memset(op,0,sizeof op);
    for(i=1;i<=n;i++){
        ans=(ans+(lc[i]*(a[i]-1)-find(a[i])-lc[i]*(lc[i]-1)/2)*(n-i-rc[i])%P)%P;
        add(a[i],a[i]);
        }
    for(i=1;i<=n;i++){
        xa=n-i-rc[i];
        ans=(ans-((long long)xa*(xa-1)*(xa-2)/6)%P+P)%P;
        }
    printf("%d",ans);
    return 0;
}