zoukankan      html  css  js  c++  java
  • bzoj 3787: Gty的文艺妹子序列

    (gty妹子好多啊。。)

    这个题整整调了快一天,真是累死了。。(这么神的题当然是在看某神犇博客的时候看到有这种做法,然后正好就做了233)

    其实就是块套树状数组,需要维护成吨的信息,,,,

    块与块之间的逆序对数,然后把每个块 i 到最后一个块 j 的逆序对和再搞到树状数组里,然后再来一个树状数组维护处到底 i 个块的时候,有多少个比 j 小的数,其实就是普通的数值插入就好了。(查询的话,块内的就用第一个树状数组求和,第二个就处理块外两部分的逆序对,然后在用第二个树状数组查分一下就可以得到块内有多少个比这个数小,或比这个数大的,答案就有了)

    对于修改,一个块自己是可以暴力重构的,对于上面的2个树状数组,因为一个块修改最短引起的是全部块(根号级别)的改动(这个块和其他块的逆序对,这个块之后比..小的数)

    所以每次修改的复杂度是log*根号的,每次查询的复杂度也是log*根号的,所以总的复杂度是N*logN*根号N的。

    (我的天呢,这种题,,,,,,,,,,简直是无脑。。。。就是码码码码码码码   while(1) cout<<"码";)

      1 #include<bits/stdc++.h>
      2 #define N 150005
      3 #define LL long long
      4 #define inf 0x3f3f3f3f
      5 #define lowbit(x) x&(-x)   
      6 using namespace std;
      7 inline int ra()
      8 {
      9     int x=0,f=1; char ch=getchar();
     10     while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
     11     while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
     12     return x*f;
     13 }
     14 int sta[105],top;
     15 void writeln(int x)
     16 {    
     17     if (!x) puts("0");
     18     else{
     19     while (x){ sta[++top]=x%10; x/=10;}
     20     while (top){ putchar(sta[top]+'0'); top--;}
     21     putchar('
    ');}
     22 }
     23 const int M=50005;
     24 int c[M],C[305][305],block_small[305][M];
     25 int n,m,Q,B,belong[M],sum[305],L[M],R[M],a[M],ans;
     26 void add(int x, int val){for (;x<=n;x+=lowbit(x)) c[x]+=val;}
     27 int ask(int x)
     28 {
     29     int sum=0;
     30     for (;x;x-=lowbit(x)) sum+=c[x];
     31     return sum;
     32 }
     33 
     34 void block_small_add(int pos, int x, int val){for (;x<=n;x+=lowbit(x)) block_small[pos][x]+=val;}
     35 int block_small_ask(int pos, int x)
     36 {
     37     int sum=0;
     38     for (;x;x-=lowbit(x)) sum+=block_small[pos][x];
     39     return sum;
     40 }
     41 
     42 void ADD(int pos, int x, int val){for (;x<=m;x+=lowbit(x)) C[pos][x]+=val;}
     43 int ASK(int pos, int x)
     44 {
     45     int sum=0;
     46     for (;x;x-=lowbit(x)) sum+=C[pos][x];
     47     return sum;
     48 }
     49 
     50 void pre()
     51 {
     52     for (int i=1; i<=m; i++) L[i]=(i-1)*B+1,R[i]=i*B; R[m]=n;
     53     for (int i=1; i<=m; i++)
     54     {
     55         int tot=0;
     56         for (int j=L[i]; j<=R[i]; j++)
     57             tot+=ask(n)-ask(a[j]),add(a[j],1);
     58         sum[i]=tot;
     59         for (int j=i+1; j<=m; j++)
     60         {
     61             tot=0;
     62             for (int k=L[j]; k<=R[j]; k++)
     63                 tot+=ask(n)-ask(a[k]);
     64             ADD(i,j,tot);
     65         }
     66         for (int j=L[i]; j<=R[i]; j++)
     67             add(a[j],-1);
     68     }
     69     for (int i=1; i<=m; i++)
     70     {
     71         for (int j=L[i]; j<=R[i]; j++)
     72             block_small_add(i,a[j],1);
     73         for (int j=1; j<=n; j++)
     74             block_small[i+1][j]=block_small[i][j];
     75     }
     76 }
     77 
     78 void ask(int x, int y)
     79 {
     80     ans=0;
     81     if (belong[x]==belong[y])
     82     {
     83         for (int i=x; i<=y; i++)
     84         {
     85             ans+=ask(n)-ask(a[i]);
     86             add(a[i],1);
     87         }
     88         for (int i=x; i<=y; i++) add(a[i],-1);
     89     }
     90     else
     91     {
     92         for (int i=belong[x]+1; i<belong[y]; i++)
     93             ans+=ASK(i,belong[y]-1)+sum[i];
     94         for (int i=x; i<=R[belong[x]]; i++)
     95         {
     96             ans+=ask(n)-ask(a[i]);
     97             ans+=block_small_ask(belong[y]-1,a[i]-1)-block_small_ask(belong[x],a[i]-1);
     98             add(a[i],1);
     99         }
    100         for (int i=L[belong[y]]; i<=y; i++)
    101         {
    102             ans+=ask(n)-ask(a[i]);
    103             ans+=block_small_ask(belong[y]-1,n)-block_small_ask(belong[y]-1,a[i])-block_small_ask(belong[x],n)+block_small_ask(belong[x],a[i]);
    104             add(a[i],1);
    105         }
    106         for (int i=x; i<=R[belong[x]]; i++) add(a[i],-1);
    107         for (int i=L[belong[y]]; i<=y; i++) add(a[i],-1);
    108     }
    109     writeln(ans);
    110 }
    111 
    112 void change(int x, int new_val)
    113 {
    114     if (new_val==a[x]) return;
    115     int last=0;
    116     for (int i=L[belong[x]]; i<x; i++)
    117     {
    118         if (a[i]>a[x]) last--;
    119         if (a[i]>new_val) last++;
    120     }
    121     for (int i=x+1; i<=R[belong[x]]; i++)
    122     {
    123         if (a[x]>a[i]) last--;
    124         if (new_val>a[i]) last++;
    125     }
    126     sum[belong[x]]+=last;
    127     for (int i=1; i<belong[x]; i++)
    128     {
    129         last=0;
    130         last-=block_small_ask(i,n)-block_small_ask(i-1,n)-block_small_ask(i,a[x])+block_small_ask(i-1,a[x]);
    131         last+=block_small_ask(i,n)-block_small_ask(i-1,n)-block_small_ask(i,new_val)+block_small_ask(i-1,new_val);
    132         ADD(i,belong[x],last);
    133     }
    134     for (int i=belong[x]+1; i<=m; i++)
    135     {
    136         last=0;
    137         last-=block_small_ask(i,a[x]-1)-block_small_ask(i-1,a[x]-1);
    138         last+=block_small_ask(i,new_val-1)-block_small_ask(i-1,new_val-1);
    139         ADD(belong[x],i,last);
    140     }
    141     for (int i=belong[x]; i<=m; i++)
    142         block_small_add(i,a[x],-1),block_small_add(i,new_val,1);
    143     a[x]=new_val;
    144 }
    145 
    146 int main()
    147 {
    148     n=ra(); B=sqrt(n); m=n/B+(n%B!=0);
    149     for (int i=1; i<=n; i++) a[i]=ra(),belong[i]=(i-1)/B+1;
    150     pre();    Q=ra();
    151     while (Q--)
    152     {
    153         int opt=ra(),x=ra()^ans,y=ra()^ans;
    154         if (opt==0) { ask(x,y);}
    155             else change(x,y);
    156     }
    157     return 0;
    158 }
  • 相关阅读:
    网络最大流算法—最高标号预流推进HLPP
    网络最大流算法—EK算法
    PROPAGATION_REQUIRED
    js左侧三级菜单导航代码
    Ubuntu上用premake编译GDAL
    2013数据结构课程设计之便利店选址(暴力枚举和随机函数两种做法)
    JAVA环境配置
    [K/3Cloud] 如何从被调用的动态表单界面返回数据
    document.getElementsByClassName在ie8及其以下浏览器的兼容性问题
    Java学习笔记51:数组转ArrayList和ArrayList转数组技巧
  • 原文地址:https://www.cnblogs.com/ccd2333/p/6435480.html
Copyright © 2011-2022 走看看