TOJ 4074 Running Laps -- 树状数组

最新推荐文章于 2023-07-23 20:13:12 发布

randygx

最新推荐文章于 2023-07-23 20:13:12 发布

阅读量656

点赞数 2

分类专栏： ACM 树状数组文章标签：树状数组

本文链接：https://blog.csdn.net/randygx/article/details/28467129

版权

ACM 同时被 2 个专栏收录

56 篇文章 0 订阅

订阅专栏

树状数组

2 篇文章 0 订阅

订阅专栏

题目链接：http://acm.tzc.edu.cn/acmhome/problemdetail.do?&method=showdetail&id=4074

题目大意：有n头牛在环形跑道赛跑，最快的牛跑L圈时停止，每头牛速度恒定，问总共有多少次"crossing event“(一头牛从另一头后面超过它，简称CE)。

分析：因为速度恒定，所以真正的超越是不可能的，只可能是快的牛超过慢的牛刚好一整圈。于是问题变为快牛总共超过慢牛多少圈。

一个很自然的想法是把牛按速度从慢到快排序(或者说，把牛按跑的圈数从小到大排序)，然后记录前面的牛跑的总圈数，也就是说，记牛i跑a[i]圈，那么牛i总共超过比它慢的牛a[i] * i - sum(i - 1)圈，其中牛下标从0开始，sum(i - 1)是前i只牛跑的总圈数。

可是问题来了：只有超过的整数圈数才真正计数，比如前3头牛分别跑1, 2, 3.5圈，实际牛2超过牛0两次，超过牛1一次，总共3次，但如果按3.5 * 2 - (1 + 2)算，就会超过4次。解决的办法是：先只考虑整数部分，并假设a[i]的小数部分小于前面所有牛的小数部分，这时牛i超过前面的牛总共((int)a[i] - 1) - sum(i - 1)，其中sum(i - 1)是前i头牛圈数的整数部分求和，a[i] - 1是因为a[i]的小数部分小，故实际少超一圈(记住，这是假设)。那么这样子我们就少算一部分，也就是实际比a[i]的小数部分小的那些牛，应该再加上1。至此，问题就转变成了前i个小数部分中有几个比a[i]的小数部分小，也就是动态名次问题，用平衡树可以解决。

比较囧的是，一开始写的SBT居然超时了……一怒之下用树状数组写，把小数部分 * 1000000，然后当成整数处理，WA，改成2000000，AC。

还有一点比较坑，就是精度问题(处理实数总是会有这种问题)，按上面所说，应该把a[i]的小数部分变为tmp = (a[i] - (int)a[i]) * 2000000，但由于精度问题，应该查找比tmp + 1小的数的总数。代码中查找tmp + 2，并插入tmp + 1，是因为树状数组下标从1开始，故把所有数加1。

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
         #include 
         #include 
         
           #include 
          
            #include 
           
             #include 
            
              #include 
             
               #include 
              
                #include 
               
                 #include 
                
                  #include 
                 
                   #include 
                  
                    #define mp make_pair #define X first #define Y second #define MEMSET(a, b) memset(a, b, sizeof(a)) using namespace std; typedef unsigned int ui; typedef long long ll; typedef unsigned long long ull; typedef pair 
                   
                     pii; typedef vector 
                    
                      vi; typedef vi::iterator vi_it; typedef map 
                     
                       mii; typedef priority_queue 
                      
                        pqi; typedef priority_queue 
                       
                         , greater 
                        
                          > rpqi; typedef priority_queue 
                         
                           pqp; typedef priority_queue 
                          
                            , greater 
                           
                             > rpqp; const int MAX_N = 100000 + 2; const int MAX_V = 2000000; const double EPS = 1.0e-8; int speed[MAX_N]; ll sum[MAX_V + 2]; inline int lowbit(int x) { return x & (-x); } void update(int p) { while (p <= MAX_V) { ++sum[p]; p += lowbit(p); } } ll query(int p) { ll ret = 0; while (p) { ret += sum[p]; p -= lowbit(p); } return ret; } int main(int argc, char *argv[]) { // freopen("D:\\in.txt", "r", stdin); int n, l, c; ll ans = 0; ll sum = 0; cin >> n >> l >> c; for (int i = 0; i < n; ++i) { scanf("%d", speed + i); } sort(speed, speed + n); for (int i = 0; i < n; ++i) { double cir = (double)l * speed[i] / speed[n - 1]; int tmp = (int)((cir - (int)cir) * 2000000); ans += ((int)cir - 1) * (ll)i - sum + query(tmp + 2); update(tmp + 1); sum += (int)cir; } printf("%I64d\n", ans); return 0; }