码迷,mamicode.com
首页 > 其他好文 > 详细

hdu4000 && hrbust1625

时间:2015-10-05 21:59:32      阅读:209      评论:0      收藏:0      [点我收藏+]

标签:

ikki的数字
Time Limit: 1000 MS Memory Limit: 32768 K
Total Submit: 22(12 users) Total Accepted: 9(8 users) Rating: 技术分享技术分享技术分享技术分享 Special Judge: No
Description

ikki 近期对数字颇感兴趣。如今ikki在纸上写了连续的N个数字,每一个数字都是[1,N]之间随意的一个数并且不反复,即这串数字

是数字1~N的一个排列。数字的序号从1到N,如今ikki想考你一下:

在这N个数字中能找出多少个3个数的组合满足:num[x]<num[z]<num[y]且x<y<z,当中x,y,z为这三个数字的下标。

Input
多组測试数据,第一行一个整数T 表示測试数据的个数。

对于每组数据,第一行输入一个整数N表示数列中数字的个数(1<=N<=5000)

第二行输入N个数字表示一个1~N的排列。


Output

对于每组数据,输出”Case #k: p” ,k表示第k组例子,p表示满足要求的3个数字的组合数目,每组输出占一行。

因为结果可能比較大,结果需对100000007取模。

Sample Input
2
6
1 3 2 6 5 4

3 5 2 4 1


Sample Output
Case #1: 10
Case #2: 1

Author

周洲@hrbust


隐藏着树状数组~~~根本没看出来,事实上主要是没思路,思路出来了才干用树状数组求解

推断满足i<j<k且num[i]<num[k]<num[j]的总组数
利用树状数组能够求出一个数前面比它小的数的个数,进而能够知道前面比它大的数的个数,总的比它大的个数减去前面比它大的个数等于后面比它大的个数,Cn2 = x*(x-1)/2;
然后肯定要减去后面的全部组成的i<j<k且a[i]<a[j]<a[k]的个数;

  1. 能够先求出(xyz,xzy)的总数量 
  2. 仅仅需出去x后面多少个比它大的个数n,C(n,2)就是了 
  3. 然后求出xyz的个数, 
  4. 对于a,求出比a小的个数low[a],比a大的个数high[a],low[a]*high[a]就是答案 
  5. 能够借助树状数组求上述个数 
注意了,总体求解!

并非说单独考虑某个数

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
typedef long long LL;
const int mod = 100000007;
int T,n;
int c[5010];
int lowbit(int x)
{
    return x & -x;
}

LL getsum(int x)
{
    LL sum = 0;
    while(x > 0)
    {
        sum += c[x];
        x -= lowbit(x);
    }
    return sum;
}


void update(int x , int val)
{
    while(x <= n)
    {
        c[x] += val;
        x += lowbit(x);
    }
}

int main()
{
#ifdef xxz
    freopen("in.txt","r",stdin);
#endif // xxz
    cin>>T;
    int Case = 1;
    while(T--)
    {

        memset(c,0,sizeof(c));
        cin>>n;
        LL ans = 0;
        for(int i = 1; i <= n; i++)
        {
            int temp;
            cin>>temp;
            update(temp,1);
            LL presmaller = getsum(temp - 1);
            LL prebigger = i-1 - presmaller;
            LL totbigger = n - temp;
            LL afterbigger = totbigger - prebigger;

            ans -= presmaller * afterbigger;

            if(afterbigger >= 2)
            {
                ans += afterbigger*(afterbigger - 1)/2;
            }

        }
        cout<<"Case #"<<Case++<<": "<<ans%mod<<endl;


    }
    return 0;
}

版权声明:本文博主原创文章,博客,未经同意不得转载。

hdu4000 &amp;&amp; hrbust1625

标签:

原文地址:http://www.cnblogs.com/hrhguanli/p/4856264.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!