hdu 3465 Life is a Line 树状数组求逆序数-CSDN博客

本文链接：https://blog.csdn.net/lois_123/article/details/47006469

题意：给出x轴上的一个区间（l，r），给出n条直线，求在该区间内有多少交点；

思路：两条直线 a，b，与 x=l这条直线的交点为（l，ya1）（l，yb1）；与x=r的交点为（r，ya2），（r，yb2）；

若a，b在（l，r）中相交，则ya1<yb1 &&ya2>yb2 || ya1>yb1&& ya2<yb2(交叉才能相交)，就是求逆序数啦

树状数组求逆序数：将数一个一个插入数组，在插入的时候就计算在该数之前（包括该数）有多少数（假设tem个）<=该数；到该数的数目总个数-tem，就是这个数的逆序数，一个一个求，一个一个加就是答案；

code：

#include <iostream>
#include <algorithm>
#include<stdio.h>
#include<string.h>
using namespace std;

const int MAX=5e4+10;

int out[MAX],ps;
struct point
{
    int i;
    double l,r;
} p[MAX];

bool cmp1(point a,point b)
{
    if(a.l!=b.l)
        return a.l <b.l;
    return a.r<b.r;
}
bool cmp2(point a,point b)
{
    if(a.r!=b.r)
        return a.r<b.r;
    return a.l<b.l;
}

int lowbit(int x)
{
    return x&(-x);
}
void modify(int pos,int num)
{
    while(pos<=ps)
    {
        out[pos]+=num;
        pos+=lowbit(pos);
    }
}
int sum(int pos)
{
    int su=0;
    while(pos>0)
    {
        su+=out[pos];
        pos-=lowbit(pos);
    }
    return su;
}

int main()
{
   // freopen("i.txt","r",stdin);
    int n,sz;
    double ll,rr,x1,y1,x2,y2;

    while(scanf("%d",&n)!=EOF)
    {
        memset(out,0,sizeof(out));
        scanf("%lf%lf",&ll,&rr);
        sz=0,ps=0;
        for(int i=0; i<n; i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            if(x2==x1)
            {
                if(ll<x1&&x1<rr)sz++;
            }
            else
            {
                p[ps].l=y1+(y2-y1)*(ll-x1)/(x2-x1);
                p[ps].r=y1+(y2-y1)*(rr-x1)/(x2-x1);
                ps++;
            }
        }

        sort(p,p+ps,cmp1);
        for(int i=0; i<ps; i++)
            p[i].i=i+1;
        sort(p,p+ps,cmp2);

        int ans=0;
        for(int i=0; i<ps; i++)
        {
            modify(p[i].i,1);
            ans+=(i+1-sum(p[i].i));
        }
        ans += sz*ps;
        printf("%d\n",ans);
    }
    return 0;
}