数位dp

洛谷 P2602 [ZJOI2010]数字计数
给定两个正整数a和b，求在[a,b]中的所有整数中，每个数码(digit)各出现了多少次。
输入格式：
输入文件中仅包含一行两个整数a、b，含义如上所述。

输出格式：
输出文件中包含一行10个整数，分别表示0-9在[a,b]中出现了多少次。
输入样例#1： 1 99
输出样例#1： 9 20 20 20 20 20 20 20 20 20

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#if 0
Writers: Goes && G.S.M.
Have falled in love
#endif
#define ll long long
const int N=15;
ll f[N][2][N][2];
int num[N];
/*
f[len][1/0][sum][1/0]
len 从高位到低位，当前是第几位
当前位置和num[len]是否相等 
sum 当前数字出现的次数
之前是否有前导0
*///dfs 中 d 表示当前在算的数码   0 true  1 false
ll dfs(int len,bool fl,int sum,bool z,int d)
{
    if(!len)return sum;ll ret=0;
    if(f[len][fl][sum][z]!=-1) return f[len][fl][sum][z];
    for(int i=0;i<=9;i++)
    if(!fl&&i>num[len])break;
    else ret+=dfs(len-1,fl||(i<num[len]),
        sum+((!z||i)&&(i==d)),z&&i==0,d);
    return f[len][fl][sum][z]=ret;
}
inline ll solve(ll x,int d)
{
    int len=0;
    while(x)    num[++len]=x%10, x/=10;
    memset(f,-1,sizeof(f));
    return dfs(len,0,0,1,d);
}

int main()
{
    ll a,b;scanf("%lld%lld",&a,&b);
    for(int i=0;i<=9;i++)
        printf("%lld ",solve(b,i)-solve(a-1,i));
    return 0;
}

洛谷 * P2606 [ZJOI2010]排列计数*

称一个1,2,…,N的排列P1,P2…,Pn是Magic的，当且仅当2<=i<=N时，Pi>Pi/2. 计算1，2，…N的排列中有多少是Magic的，答案可能很大，只能输出模P以后的值
输入样例#1： 复制
20 23
输出样例#1： 复制
16

100%的数据中，1 ≤N ≤ 10^6, P≤ 10^9，p是一个质数。

【此题题解为粘，不过能看懂，借用了。】
题意：求1到n的所有排列中，满足小根堆性质的排列的个数。

这里建立DP模型：f[i]表示i个不同的数的所有排列中满足小根堆性质的排列的个数。

对于转移，首先计算出ii个节点的完全二叉树中，根节点的左子树包含的节点数l，右子树包含的节点数r。

首先，根节点的值必须为最小值。再考虑剩下的i-1i−1个节点。很容易想到，可以在这i-1i−1个节点中选出ll个节点作为左子树，剩下的rr个节点作为右子树。所以得出转移：
f [ i ] = C ( l , i - 1 ) * f [ l ] * f [ r ]

同时注意坑点(在BZOJ上)：n可以大于p，所以求组合数要用到LucasLucas定理。

#include<iostream>
#include<cstring>
#include<cstdio>
#define ll long long
using namespace std;
inline int read(){
    char ch='*';int sum=0,f=1;
    while(ch>'9'||ch<'0'){
        if(ch=='-')f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0')
        sum*=10,sum+=ch-'0',ch=getchar();
    return sum*f;
}
#if 0
Writers: G.S.M. && Goes.
Fall in love.
#endif
const int N=1e6+5;
int log[N],inv[N];
int f[N],fac[N],n,MO;
int qpow(int a,int b){
    int res=1;
    while(b){
        if(b&1)res=(ll)1*res*a%MO;
        a=(ll)1*a*a%MO;b>>=1;
    }return res;
}
int C(int x,int y){
    if(!y)return 1;int u=C(x/MO,y/MO),v=x%MO,w=y%MO,z;
    if(v<w)z=0;else z=(ll)1*((ll)1*
        fac[v]*inv[w]%MO)*inv[v-w]%MO;
    return (ll)1*u*z%MO;
}
int main()
{
    n=read();MO=read();log[0]=-1;fac[0]=1;
    for(int i=1;i<=n;i++)
        fac[i]=(ll)1*fac[i-1]*i%MO,
        log[i]=log[i>>1]+1;
    int kx=min(MO-1,n);inv[kx]=qpow(fac[kx],MO-2);
    for(int i=kx-1;i>=0;i--)
        inv[i]=(ll)1*inv[i+1]*(i+1)%MO;
    f[1]=f[2]=1;f[3]=2;int l=1,r=1;
    for(int i=4;i<=n;i++){
        if(i-(1<<log[i])+1<=(1<<log[i]-1))l++;
        else r++;
        f[i]=(ll)1*((ll)1*C(i-1,l)*f[l]%MO)*f[r]%MO;
    }printf("%d\n",f[n]);return 0;
}