hdu6044 Limited Permutation【读入优化+dfs】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6044
题意：有长度为n的序列，分别表示l[i]和r[i]，然后让你构造一个a序列，a[i]满足在[l[i],r[i]]区间里面，a[i]为最小值，满足这样条件的序列有多少种
解析：首先根据题面的限定，肯定存在一个区间i是[1,n]，对于这个区间，可以说第i个值已经确定了，因为这个区间是当前的整个区间，对于剩下的来说，肯定从剩下的n-1个数中选i-1个数出来填在左边，而剩下的填在右边，接下来在乘上左右区间的各自的方法数，总的计数方案应该为ans = dfs(1,i-1)*dfs(i+1,n) *C(r-l,i-l)，所以需要递归往下划分，此时再往下划分的话，就会把这个区间分为[1,i-1]和[i+1,n]，依次递归向下，就可以吧所有的区间扫一遍

#include <bits/stdc++.h>
using namespace std;
const int maxn = 3e6+100;
const int mod = 1e9+7;
typedef long long ll;
namespace fastIO
{
    #define BUF_SIZE 100000
    //fread -> read
    bool IOerror = 0;
    inline char nc()
    {
        static char buf[BUF_SIZE], *p1 = buf + BUF_SIZE, *pend = buf + BUF_SIZE;
        if (p1 == pend)
        {
            p1 = buf;
            pend = buf + fread(buf, 1, BUF_SIZE, stdin);
            if (pend == p1)
            {
                IOerror = 1;
                return -1;
            }
        }
        return *p1++;
    }
    inline bool blank(char ch)
    {
        return ch == ' ' || ch == '\n' || ch == '\r' || ch == '\t';
    }
    inline void read(int &x)
    {
        char ch;
        while (blank(ch = nc()));
        if (IOerror)
            return;
        for (x = ch - '0'; (ch = nc()) >= '0' && ch <= '9'; x = x * 10 + ch - '0');
    }
    #undef BUF_SIZE
};
using namespace fastIO;
struct node
{
    int l,r,id;
    bool operator < (const node &b)const
    {
        if(l==b.l)
            return r>b.r;
        return l<b.l;
    }
}a[maxn];
ll h[maxn];
ll qpow(ll x,ll n)
{
    ll res = 1;
    while(n)
    {
        if(n&1)
            res = res*x%mod;
        x = x*x%mod;
        n >>= 1;
    }
    return res;
}
ll C(int n,int m)
{
    ll a = h[n];
    ll b = h[m]*h[n-m]%mod;
    return a*qpow(b,mod-2)%mod;
}
int id;
ll dfs(int l,int r)
{
    if(a[id].l!=l || a[id].r!=r)
        return 0;
    int last = a[id].id;
    id++;
    ll ansl = 1,ansr = 1;
    if(l<=last-1)
        ansl = dfs(l,last-1);
    if(r>=last+1)
        ansr = dfs(last+1,r);
    return ansl*ansr%mod*C(r-l,last-l)%mod;
}
int main(void)
{
    h[0] = 1;
    for(int i=1; i<maxn; i++)
        h[i] = (h[i-1]*i)%mod;
    int n;
    int case_t = 1;
    while(read(n),!fastIO::IOerror)
    {
        for(int i=1;i<=n;i++)
            read(a[i].l);
        for(int i=1;i<=n;i++)
        {
            read(a[i].r);
            a[i].id = i;
        }
        sort(a+1,a+n+1);
        id = 1;
        ll ans = dfs(1,n);
        printf("Case #%d: %I64d\n",case_t++,ans);
    }
    return 0;
}