PKU3067 Japan - 树状数组

题目描述：

城市被分为东西两个部分，两岸城市之间有高速公路连接。给出高速公路连接情况，判断有公路一共有多少个交点。

分析：

类似二部图匹配的模型。

设公路左边城市为x，右边城市为y。首先按照x城市的编号从大到小将公路a[]排序。

对于公路a[i]交点个数为，0到i-1之间的公路y<a[i].y的个数。（描述不清￥%#&……不说了，代码很明白）

这里是典型的树状数组的应用。

要注意两个地方：第一，公路可能有很多条，数组要开到10^6；第二，求和可能超int，要用int64。

 

1. /*
2. PKU3067 Japan
3. */
5. #include <stdio.h>
6. #include <stdlib.h>
7. #include <memory.h>
9. #define clr(a) memset(a,0,sizeof(a))
10. #define int64 __int64
11. #define I64d "%I64d"
13. #define M 1000000
14. #define N 3000
16. typedef struct{
17.     int x,y;
18. }Node;
20. int n,m;
21. Node a[M];
22. int c[N];
24. int cmp(const void *A,const void *B){
25.     Node a = *(Node*)A;
26.     Node b = *(Node*)B;
27.     if(b.x==a.x) return b.y-a.y;
28.     return b.x-a.x;
29. }
31. int lowbit(int x){
32.     return x&(-x);
33. }
35. void Add(int i){
36.     while(i<=n){
37.         c[i]++;
38.         i+=lowbit(i);
39.     }
40. }
42. int64 Sum(int i){
43.     int64 s=0;
44.     while(i>0){
45.         s+=c[i];
46.         i-=lowbit(i);
47.     }
48.     return s;
49. }
51. int main()
52. {
53.     int T,Tn;
54.     
55.     scanf("%d",&Tn);
56.     
57.     for(T=1;T<=Tn;T++){
58.         int i,j,k;
59.         //input
60.         scanf("%*d%d%d",&n,&m);
61.         for(i=0;i<m;i++){
62.             scanf("%d%d",&a[i].x,&a[i].y);
63.         }
64.         //sort
65.         qsort(a,m,sizeof(a[0]),cmp);
66.         
67.         //work
68.         int64 sum=0;
69.         clr(c);
70.         
71.         for(k=0;k<m;k++){
72.             i=a[k].y;
73.             Add(i);
74.             sum+=Sum(i-1);
75.         }
76.         //output
77.         printf("Test case %d: "I64d"/n",T,sum);
78.     }
79.     
80.     return 0;
81. }