树状数组求逆序对 基础版+强化版本(去重+离散化)

 

是getsum(a[i]-1)

不是getsum(i-1);

20210920

复制代码
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5;
int a[maxn],n,c[maxn],ans[maxn];
int lowbit(int x)
{
    return x&(-x);
}
int update(int id,int val)
{
    a[id]=val;
    while(id<=n)
     {  
         c[id]+=val;
         id+=lowbit(id);
     }
}
int getsum(int x)
{ int ans=0;
    while(x>0)
     {
         ans+=c[x];
         x-=lowbit(x);
     }
     return ans;
}
int main( )
{

   cin>>n;
   for(int i=1;i<=n;i++)
    {
        cin>>a[i];
    }
    

    
    for(int i=1;i<=n;i++)
    {
        update(a[i],1);
        ans[i]=i-getsum(a[i]-1);
        cout<<ans[i]<<" ";
    }    
    
}
复制代码

 

强化版本,自己想的+手写的,因为看不懂也不想看洛谷那群家伙写的代码

 

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6;
long long b[maxn],n,c[maxn],vis[maxn];
struct node{
    long long val,id;
}a[maxn];
bool cmp(node a,node b)
{
    return a.val<b.val;
}
int lowbit(int x)
{
    return x&(-x);
}
void update(int id,int val)
{
    while(id<=n)
     
        c[id]+=val;
        id+=lowbit(id);
     }
}
long long getsum(int x)
{ long long ans=0;
    while(x>0)
     {
        ans+=c[x];
        x-=lowbit(x);
     }
     return ans;
}
void pre( )
{
  memset(vis,0,sizeof(vis));
  memset(c,0,sizeof(c));
}
int main( )
{
 
 
    
     
//freopen("lyslys.in","r",stdin);
int t;
cin>>t;
while(t--)
{
    pre( );
   cin>>n;
    for(int i=1;i<=n;i++)
    {   cin>>a[i].val;
        a[i].id=i;
        //cout<<a[i]<<endl;
    /*  update(a[i],1);
        int x=i-getsum(a[i]-1)-1;
       /* int y;
        y=i-1-x;
        if(x>y)
        {
            ans[i]=1;
            //逆序的更多,放前面贪心
        }
        else {
            ans[i]=0;
        }
        //这边要不要讨论==0的情况
        */
    }
         
    sort(a+1,a+1+n,cmp);
     
     
     int i,j;
     long long cnt=0;
      
    for(i=1;i<=n;)
    {   cnt++;
        for(j=i+1;j<=n;j++)
         {
            if(a[j].val==a[i].val)
             {
               a[j].val=cnt;
             }
             else break;
         }
         a[i].val=cnt;
         i=j;
    }
     
    for(int k=1;k<=n;k++)
    {
        b[a[k].id]=a[k].val;
    }
     
     
    cnt=0;
     
     
    for(int i=1;i<=n;i++)
    {  
        //cout<<a[i]<<endl;
        update(b[i],1);
        vis[b[i]]++;
         
        long long x,y;
         
        x=i-getsum(b[i]-1)-vis[b[i]];
          
        y=getsum(b[i]-1);
         
       // cout<<x<<" "<<y<<endl;
         
        cnt+=min(x,y);
       /* int y;
        y=i-1-x;
        if(x>y)
        {
            ans[i]=1;
            //逆序的更多,放前面贪心
        }
        else {
            ans[i]=0;
        }
        //这边要不要讨论==0的情况
        */
       //cout<<x<<" ";
    }
    printf("%lld\n",cnt);
}
 
}

  

 

posted @   liyishui  阅读(45)  评论(0编辑  收藏  举报
编辑推荐:
· Linux系列:如何用heaptrack跟踪.NET程序的非托管内存泄露
· 开发者必知的日志记录最佳实践
· SQL Server 2025 AI相关能力初探
· Linux系列:如何用 C#调用 C方法造成内存泄露
· AI与.NET技术实操系列(二):开始使用ML.NET
阅读排行:
· 无需6万激活码!GitHub神秘组织3小时极速复刻Manus,手把手教你使用OpenManus搭建本
· C#/.NET/.NET Core优秀项目和框架2025年2月简报
· Manus爆火,是硬核还是营销?
· 终于写完轮子一部分:tcp代理 了,记录一下
· 【杭电多校比赛记录】2025“钉耙编程”中国大学生算法设计春季联赛(1)
点击右上角即可分享
微信分享提示