HDU6044 Limited Permutation[快速输入][分治]

题意：给定$i$从1到n个区间左右端点，描述的第i个区间的最小值必须是第i个数字，问从1~n有多少种符合题意的排列。
题目提示输入使用快速的fread()
对于一个确定的区间$l_i, r_i$，其中$p_i$是其最小值为定数，其余的可能情况就可以按照组合数去组合放到两边，所以从1~n大区间开始分治，对于每个区间，使用stl::map记录区间的编号 i <script type="math/tex" id="MathJax-Element-243">i</script>。
对于用到的组合数结果，使用卢卡斯预处理出结果。

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1e6+10;
const int mod = 1e9+7;
int n;
int l[N], r[N];
map<pair<int ,int>, int> mp;
LL ans;
LL fac[N], inv[N];
inline char nc()
{
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}

template <class T>
inline bool scan_d(T & x)
{
    char c=nc();
    x=0;
    if(c==EOF)return false;
    for(; c>'9'||c<'0'; c=nc());
    for(; c>='0'&&c<='9'; x=x*10+c-'0',c=nc());
    return true;
}

void init()
{
    fac[0]=1;
    for(int i = 1; i < N; i++)
        fac[i] = (fac[i-1]*i) % mod;

    inv[0] = inv[1] = 1;
    for(int i=2; i<N; i++)
        inv[i]=(LL)(mod-mod/i)*inv[mod%i]%mod;
    for(int i=1; i<N; i++)
        inv[i]=(inv[i-1]*inv[i])%mod;
}

LL C(int n,int m)
{
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}

void dfs(int ll,int rr)
{
    if(ans==0 || rr<ll) return;
    int i = mp[make_pair(ll, rr)];
    if(i == 0)
    {
        ans=0;
        return;
    }
    if(ll == rr)return;
    int len = rr - ll;
    ans = (ans*C(len,i-ll)) % mod;
    dfs(ll, i-1);
    dfs(i+1, rr);
}
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("input.txt", "r", stdin);
    #endif // ONLINE_JUDGE
    init();
    for(int Case = 1; scan_d(n); ++Case)
    {
        mp.clear();
        for(int i = 1; i <= n; ++i)
            scan_d(l[i]);
        for(int i = 1; i <= n; ++i)
        {
            scan_d(r[i]);
            mp[make_pair(l[i], r[i])] = i;
        }
        ans = 1;
        dfs(1, n);
        printf("Case #%d: %I64d\n", Case, ans);
    }
    return 0;
}