上机一 G D&C--玲珑数

D&C--玲珑数

时间限制：1000ms   内存限制：65536kb

通过率：9/38 (23.68%)    正确率：9/155 (5.81%)

题目描述

Bamboo听说知道了正确的玲珑数，就能开启玲珑阵，穿越时空，然后解锁一个满分玛丽苏cp……

于是决定停止幻想，现在就去算玲珑数。

玲珑数是什么呢？

在数列中任意两个数a[i],a[j],如果i<j且a[i]>2*a[j]，那么a[i]和a[j]就构成了一对儿玲珑对。玲珑数就是玲珑对的总数。

是不是听起来很简单？快帮Bamboo算一下玲珑数是多少吧！

输入

多组输入数据

输入的第一个数是正整数n，表示数组长度，0<n<=10000。

接下来n个正整数是数组an，均在int范围内。

然后是正整数t表示t次查询 ，0< t<1000

接下来t行，每行两个数p和q，代表查询的数组下标区间 0<=p<n, 0<=q<n。

输出

每次查询输出一行，为所求区间内的玲珑数。

输入样例

3
3 2 1
1
0 2

输出样例

1

解析：

此题与求解逆序数类似，同样在归并排序的代码基础上略做改动即可。

将数组以q = (p+r)/2划分，总玲珑数 = 左数组玲珑数+右数组玲珑数+跨数组玲珑数，左、右数组玲珑数可以通过递归得出，关键在于如何求出跨数组的玲珑数。

以下求解跨数组玲珑数：由于归并排序，L[]和R[]是已经排好序的两个数组，对L数组每个元素L[i]循环：若L[i] >R[j]*2，j++继续判断；直到L[i] <=R[j]*2，ans+=j。

最后输出结果ans即可。

求解逆序数：445 模式寻对

代码：

#include<cstdio>
#define maxn 10007

long long ans;
long long a[maxn],A[maxn];

void Merge(int p,int q,int r)
{
    int n1 = q-p+1;
    int n2 = r-q;
    long long L[n1],R[n2];
    for(int i = 0; i < n1; i++)
    {
        L[i] = A[p+i];
    }
    for(int i = 0; i < n2; i++)
    {
        R[i] = A[q+1+i];
    }

    for(int i = 0,j = 0;i < n1;i++)
    {
        while(j < n2&&L[i] > R[j]*2)
        {
            j++;
        }
        ans += j;
    }

    int i = 0,j = 0,k = p,l = 0;
    while(i < n1&&j < n2)
    {
        if(L[i] > R[j])
        {
            A[k++] = R[j++];
        }
        else
        {
            A[k++] = L[i++];
        }
    }
    while(i < n1)
    {
        A[k++] = L[i++];
    }
    while(j < n2)
    {
        A[k++] = R[j++];
    }
}

void Merge_Sort(int p,int r)
{
    if(p < r)
    {
        int q = (p+r)/2;
        Merge_Sort(p,q);
        Merge_Sort(q+1,r);
        Merge(p,q,r);
    }
}

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        for(int i = 0; i < n; i++)
        {
            scanf("%lld",&a[i]);
        }
        int t;
        scanf("%d",&t);
        while(t--)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            int k = 0;
            if(x > y)
            {
                int hold = x;
                x = y;
                y = hold;
            }
            for(int i = x; i <= y; i++)
            {
                A[k++] = a[i];
            }
            ans = 0;
            Merge_Sort(0,y-x);
            printf("%lld\n",ans);
        }
    }
}