7.29日结

P1637 三元上升子序列
题目

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
int n,m,t1[30010],t2[30010],a[30010],b[30010],l[30010],r[30010];
inline int lowbit(int x){
    return x&(-x);
}
inline int pos(int x){
    return lower_bound(b+1,b+m+1,x)-b;
}
void add(int *t,int x,int k){
    for(int i=x;i<=m;i+=lowbit(i))
        t[i]+=k;
}
long long ask(int *t,int x){
    long long ans=0;
    for(int i=x;i>0;i-=lowbit(i))
        ans+=t[i];
    return ans;
}
int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        b[i]=a[i];
    }
    sort(b+1,b+n+1);
    m=unique(b+1,b+n+1)-(b+1);
    for(int i=1;i<=n;i++)
    {
        add(t1,pos(a[i]),1);
        l[i]=ask(t1,pos(a[i])-1);
    }
    for(int i=n;i>=1;i--)
    {
        add(t2,pos(a[i]),1);
        r[i]=n-i+1-ask(t2,pos(a[i]));
    }
    long long ans=0;
    for(int i=2;i<n;i++)
        ans+=l[i]*r[i];
    cout<<ans;
	return 0;
}

在不是主函数的地方写函数的时候，如果该函数要应用于多个数组，可以加个*表示取自，上述代码就是例子。
inline表示内联函数，对于上述两个较短的可以使用。
lower_bound(b+1,b+m+1,x)-b;返回值是x在b 中的位置，话句话说就是b数组的下标，与lower_bound括号里包括的起始位置无关，并不是从始位置开始数到第几个。
m=unique(b+1,b+n+1)-(b+1);去重后的个数，与lower_bound相比一个是减b一个是减(b+1)
l[i]=ask(t1,pos(a[i])-1);的解释：pos(a[i])是大于等于a[i]的第一个位置，,pos(a[i])-1即比a[i]小的最大一个数。
add(t1,pos(a[i]),1);与l[i]=ask(t1,pos(a[i])-1);同时进行，保证了找到的数一定在他前面，因为更新一个计算一个，如果更新完在计算则计算了所有比他小的数没有保证这些数在他前面，右边的情况同理。
r[i]=n-i+1-ask(t2,pos(a[i]));中n-i+1表示从i 到n 数的个数，包括了第 i 个数，add从后往前更新，减去了后面的数中/小于等于这个数/的个数，得到的就是严格大于这个数/的个数。
对于第 i 个数，大于它的数的个数成小于它的数的个数得到第 i 个数处于中位数构成的三元上升子序列的个数。

P6278 [USACO20OPEN]Haircut G
题目

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
long long n,m,t[100010],a[100010],s[100010],ans;
inline long long lowbit(long long x){
    return x&(-x);
}
void add(int x,int k){
    for(int i=x;i<=n;i+=lowbit(i))
        t[i]+=k;
}
long long ask(int x){
    long long ans=0;
    for(int i=x;i>0;i-=lowbit(i))
        ans+=t[i];
    return ans;
}
int main()
{
    cin>>n;
	for(int i=1;i<=n;i++)
    {
        scanf("%lld",&a[i]);
        a[i]++;
    }
    for(int i=1;i<=n;i++)
    {
        long long x=n-a[i]+2;
        s[a[i]]+=ask(x-1);
        add(x,1);
    }
    cout<<'0'<<endl;
    for(int i=2;i<=n;i++)
    {
        ans+=s[i-1];
        printf("%lld\n",ans);
    }
	return 0;
}

题目有一个尤为重要的条件0≤Ai≤N，有Ai可能为0而树状数组的不能有0，所以要先把数组都加一。
0≤Ai≤N这个条件也决定了输出的方法，代码中输出的思路是把i 看作时间，每一秒fj的头发会长一微米，第 i 个头发在Ai 秒之和就不再生长。
这题我个人感觉有点类似桶排序的思想，但是用了x=n-a[i]+2;反向使得小的在后大的在前，因为树状数组的求和计算的是前 i 个的和，所以我们反向更新找比他大的数。
因为树状数组是个比较抽象的东西，难以理解（对我来说），我认为它的核心就是减少普通数组的单点/区间的修改/查询的时间，所以我的理解方法是就把他先看做普通的数组a[i]（即不看做子节点的和仅表示第 i 个的值），之后在求和再换成用t表示的树状数组，就比如上述的s[a[i]]+=ask(x-1);对于第 i 个头发s[a[i]]表示包括Ai的逆序数（Ai为小，另一个数为大），ask(x-1)就是比他大的数的个数。