zoukankan      html  css  js  c++  java
  • 【雅礼联考DAY02】Magic

    #include<cstdio>
    #include<iostream>
    #include<cstring>
    #include<algorithm>
    using namespace std;
    typedef long long LL;
    
    const int N = 1e5;
    int n , lim , x[N + 5] , y[N + 5] , z[N + 5] , sum;
    int bl[(N << 1) + 5] , c[(N << 1) + 5] , p[(N << 1) + 5];
    struct point{
    	int l , r , bl;
    }a[N + 5];
    LL ans;
    
    inline int read()
    {
    	char ch = getchar();
    	int res = 0;
    	for(; ch < '0' || ch > '9'; ch = getchar());
    	for(; ch >= '0' && ch <= '9'; ch = getchar()) res = (res << 3) + (res << 1) + ch - '0';
    	return res;
    }
    inline int lowbit(int x) {return x & (-x);}
    inline bool cmp(point x , point y){return x.r < y.r;}
    
    inline void add(int x , int v)
    {
    	for(; x <= lim; x += lowbit(x)) c[x] += v;
    }
    
    inline int query(int x)
    {
    	int res = 0;
    	for(; x; x -= lowbit(x)) res += c[x];
    	return res;
    }
    
    int main()
    {
    //	freopen("Magic.in" , "r" , stdin);
    	scanf("%d" , &n);
    	for(register int i = 1; i <= n; i++)
    	{
    		a[i].l = read() , a[i].r = read();
    		if (a[i].l > a[i].r) swap(a[i].l , a[i].r);
    		a[i].bl = bl[a[i].l] = bl[a[i].r] = i;
    		p[a[i].l] = 1 , p[a[i].r] = 2;
    	}
    	lim = n << 1;
    	
    	sum = 0;
    	for(register int i = 1; i <= lim; i++)
    	{
    		if (p[i] == 1) x[bl[i]] += sum;
    		else sum++;
    	}
    	
    	sum = 0;
    	for(register int i = lim; i; i--)
    	{
    		if (p[i] == 2) x[bl[i]] += sum;
    		else sum++;
    	}
    	
    	sort(a + 1 , a + n + 1 , cmp);
    	for(register int i = n; i; i--) 
    	{
    		x[a[i].bl] += query(a[i].l);
    		add(a[i].l , 1);
    	}
    	memset(c , 0 , sizeof c);
    	for(register int i = 1; i <= n; i++)
    	{
    		y[a[i].bl] = query(a[i].r) - query(a[i].l);
    		add(a[i].l , 1);
    	}
    	
    	for(register int i = 1; i <= n; i++) z[i] = n - x[i] - y[i] - 1;
    	for(register int i = 1; i <= n; i++)
    		ans += 1LL * x[i] * y[i] + 1LL * (x[i] + y[i]) * z[i] / 2;
    	
    	ans = 1LL * n * (n - 1) / 2 * (n - 2) / 3 - ans;
    	printf("%lld" , ans);
    }
    
  • 相关阅读:
    拥塞避免
    计算机网络常考
    [CODEVS1014]装箱问题
    [CODEVS2055]集合划分
    [CODEVS3641]上帝选人
    [GRYZ2014]递增子序列最大和
    [GRYZ2014]最大连续子序列的和
    金矿模型看动归
    [CODEVS1220]数字三角形
    [CODEVS1294]全排列
  • 原文地址:https://www.cnblogs.com/leiyuanze/p/12337062.html
Copyright © 2011-2022 走看看