POJ 3067 Japan(树状数组)

题目链接:

http://poj.org/problem?id=3067

题意：

给你两个集合S,R，分别有标号为1-n和标号1-m的点，并且告诉你k条边是由S与R中点连接得到的，问你这些边共多少个交点；

分析：

看了别人说是逆序数，写的时候也没注意到；首先将边按S集合点标号升序排序，对于排序后边集，当遍历第一条边(a,b)时，我们将a与b连接，那么显然如果对于下一条边(c,d)，若d<b那么显然这时候交点是一个(即连接到$[1,b-1]$这个区间里的任何一个点都会产生一个交点)，我们再连接下一条边(e,f)，那么连接到$[1,f-1]$上所产生的交点会在原来基础上增加1，相当于连接$[1,f-1]$上的每个点产生的交点数都加上了1，依次类推，就是一个区间修改，区间求和问题了；建议画图模拟一下；注意(a,b)和(a,c)以及(a,b)和(c,b)这两种情况是不产生交点的；
这里和求逆序数的考虑的不一样,仅是个人片面见解

#include<algorithm>
#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#include<vector>
#include<queue>
#include<stack>
#include<cmath>
#include<set>
#include<map>

using namespace std;

#define  inf 0x7f7f7f7f
#define  maxn 1104100
#define  N 2
#define mod 100003

typedef long long ll;
typedef struct {
    int u, v, next, w;
} Edge;
Edge e[N];
int cnt, head[N];

inline void add(int u, int v) {
    e[cnt].u = u;
    e[cnt].v = v;
    //e[cnt].w=w;
    // e[cnt].f=f;
    e[cnt].next = head[u];
    head[u] = cnt++;
    //e[cnt].u=v;
//    e[cnt].v=u;
//    e[cnt].w=0;
//    e[cnt].f=-f;
//    e[cnt].next=head[v];
//    head[v]=cnt++;
}

inline void write(int x) {
    if (x < 0)
        putchar('-'), x = -x;
    if (x > 9)
        write(x / 10);
    putchar(x % 10 + '0');
}

inline int read() {
    int x = 0, f = 1;
    char c = getchar();
    while (c < '0' || c > '9') {
        if (c == '-')
            f = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x * f;
}
struct point{
    int x,y;
    bool operator <(const point &a)const{
        if(a.x!=x)return a.x>x;
        else return a.y>y;
    }
};

point p[maxn];
ll maxm,t,n,m,k,a[1330];
void insert(int x,int val){
    while(x<=1328){
        a[x]+=val;
        x+=x&-x;
    }
}
ll sum(int x){
    ll res=0;
    while(x){
        res+=a[x];
        x-=x&-x;
    }
    return res;
}

int main() {
    cin>>t;
    for(int i=1;i<=t;i++){
        cin>>n>>m>>k;
        memset(a,0,sizeof(a));
        for(int j=0;j<k;j++){
            scanf("%d%d",&p[j].x,&p[j].y);
        }
        sort(p,p+k);
        ll ans=0;
        for(int j=0;j<k;j++){
            ans+=sum(p[j].y);
            //cout<<sum(1)<<" "<<sum(2)<<" "<<sum(3)<<" "<<sum(4)<<endl;
            insert(1,1);//这两句就是树状
            insert(p[j].y,-1);//数组的区间修改

        }
        cout<<"Test case "<<i<<": ";
        cout<<ans<<endl;
    }
    return 0;
}

我们坚持一件事情，并不是因为这样做了会有效果，而是坚信，这样做是对的。
——哈维尔