玲珑学院 1132 - Down the Rabbit Hole(逆序数）

5.632
我（或者是在读这篇文字的你）不属于这个世界
这是世界的界限
6.41
世界的意义必定存在于世界之外
世界中的一切事物如其所存在般而存在，如其所发生般而发生
世界之中不存在价值

果然……好女人要有的是，烟、楼顶……还有轻飘飘的衣服呀……

某天，水上由岐看见天上掉下的毛绒玩具。
“被天空接受”那是为了寻找不知何时开始在这个城市流传的“回归天空之路”的行为。

为了被天空接受而被扔出去的木偶，在空中飞舞并最终坠落。
那是为了将其本身即为世界的少女送予天空的少女的行为。
横跨天河的vega与Altair，被称为织女星与牛郎星的两颗星星,再加上北十字星的顶之星deneb，被称为夏季大三角。 那是形容三位一体的圣之图形。
只有夏季大三角在天空闪耀之时，世界才与天空相遇。

我想试一试，第一次，也是最后一次的恶作剧

您正在打galgame，然后突然发现您今天太颓了，于是想写个数据结构题练练手：
说到数据结构，当然要说红黑树啦
有一个红黑树，n个点，然而由于写丑了，所以里面的点权是乱的
红黑树点权当然不能是乱的啊，所以这只是一道题而已，没有让您写红黑树，不用担心
m个操作：
1.询问一条链x -> y中的逆序对对数
2.修改x点的点权为y
所谓逆序对就是说在一个序列中，i < j , a[i] > a[j] 则( i , j )为一个逆序对
1操作可以理解为把x -> y这条链当作一个序列，求这个序列的逆序对对数

！！！！注意，这个树是个红黑树哦，红黑树有特殊的性质哦！！！！
！！！！红黑树的树高是多少大家仔细回想一下！！！！！！！！！！
！！！！！！！！！！红黑树树高是O( logn )的！！！！！！！！！！

INPUT
第一行两个数n,m
第二行n个数，表示每个节点的点权ai
之后n-1行，每行两个点x,y表示x点和y点之间连有一条边
之后m行，每行为opt x y
如果opt为1，表示江x节点的值修改为y
如果opt为2，表示求x -> y这条链上的逆序对对数
答案对19260817取膜
OUTPUT
对于每个询问，输出一个数表示答案
SAMPLE INPUT
5 5
1 2 3 4 5
1 2
1 3
3 4
3 5
2 4 5
2 1 5
2 3 5
1 5 1
2 4 5
SAMPLE OUTPUT
1
0
0
3
题解：因为树高是logn的，所以两个点之间的链不会很长，我们就可以暴力求逆序数。
代码：

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<vector>
#define ll long long
using namespace std;
const int N=1e5+10;
ll a[N];
int cnt,u,v,k,tot,n,q,t;
int m,x,y,opt;
ll rt;
int dep[N],fa[N],siz[N],son[N];
struct node{int to,next;}edge[N<<2];
int head[N<<2];
void add(int u,int v)
{
    edge[cnt].to=v,edge[cnt].next=head[u],head[u]=cnt++;
    edge[cnt].to=u,edge[cnt].next=head[v],head[v]=cnt++;
}
void dfs1(int u,int f,int d)
{
    fa[u]=f;
    dep[u]=d,siz[u]=1,son[u]=0;
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        int v=edge[i].to;
        if(v==f)continue;
        dfs1(v,u,d+1);
        siz[u]+=siz[v];
        if(siz[son[u]]<siz[v]) son[u]=v;
    }
}
int top[N],pre[N],tree[N];
void dfs2(int u,int tp)
{
    top[u]=tp,tree[u]=++tot,pre[tree[u]]=u;
    if(!son[u]) return ;
    dfs2(son[u],tp);
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        int v=edge[i].to;
        if(v==fa[u]||v==son[u]) continue;
        dfs2(v,v);
    }
}
void solve(int x,int y)
{
vector<ll>vx,vy;
 while(x!=y)
{
    if(dep[x]>=dep[y])
    {
        vx.push_back(a[x]);
        x=fa[x];
    }
    else
    {
        vy.push_back(a[y]);
        y=fa[y];
    }
}
vx.push_back(a[x]);
int sz=vy.size();
for(int i=sz-1;i>=0;i--)
{
    vx.push_back(vy[i]);
}
sz=vx.size();
int ans=0;
for(int i=0;i<sz;i++)
{
    for(int j=i+1;j<sz;j++)
    {
        if(vx[j]<vx[i]) ans++;
    }
}
printf("%d\n",ans);
}
void init()
{
    tot=0;
    cnt=0;
    memset(head,-1,sizeof(head));
    memset(fa,0,sizeof(fa));
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
    init();
    for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
    for(int i=1;i<n;i++) scanf("%d%d",&u,&v),add(u,v);
      dfs1(1,0,1);
      while(m--)
      {
          scanf("%d",&opt);
          if(opt==1)
         {
             scanf("%d%lld",&x,&rt);
             a[x]=rt;
         }
         else
         {
             scanf("%d%d",&x,&y);
             solve(x,y);
         }
      }
    }
}