【组合计数问题小总结】

鉴于自己每次见到有关组合计数（求方案数之类）的题都要跪，所以专门开个网页记录总结一下。。。不定期更新

【BZOJ1079】记忆化搜索，由于每一种颜色可以涂得次数比较小，所以可以将涂的次数相同的颜色看成一个等价类，记忆化搜索转移即可。

【BZOJ2111】这道题的条件其实已经给的非常裸了。。。等价于给定一棵二叉树，左儿子的编号是当前编号的两倍，右儿子的编号是左儿子编号加一，要求给每个点标一个权值，使得任何一棵子树中根节点编号最小，从叶子节点开始合并即可，注意有部分情况需要用LUCAS（变态数据）

【BZOJ3505】从总数中减去三点共线的方案，枚举三点共线情况中最外面两点横纵坐标差的绝对值，依次剪去即可。

/**************************************************************
    Problem: 2111
    User: RicardoWang
    Language: C++
    Result: Accepted
    Time:968 ms
    Memory:40340 kb
****************************************************************/
 
#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
#define MAXN 1000005
long long N,P,sz[MAXN*2],inv[MAXN],fac[MAXN],f[MAXN];
void ready()
{
    fac[0]=fac[1]=1;
    for(int i=2;i<MAXN;i++)fac[i]=fac[i-1]*i%P;
    inv[0]=inv[1]=1;
    for(int i=2;i<P && i<=N;i++)inv[i]=(P-P/i)*inv[P%i]%P;
    for(int i=1;i<P && i<=N;i++)inv[i]=inv[i]*inv[i-1]%P;
    return ;
}
long long C(long long x,long long y)
{
    if(x<y)return 0;
    if(x<P && y<P)
    {
        return fac[x]*inv[y]%P*inv[x-y]%P;
    }
    else return C(x/P,y/P)*C(x%P,y%P)%P;
}
int main()
{
    //freopen("in.txt","r",stdin);
    cin>>N>>P;
    ready();
    for(int i=N;i;i--)
    {
        sz[i]=sz[i<<1]+sz[(i<<1)|1]+1;
        f[i]=C(sz[i]-1,sz[i<<1])*((i<<1) > N ? 1 : f[i<<1])%P*( (i<<1|1) > N ?1 : f[i<<1|1]) % P;
    }
    cout<<f[1]<<endl;
    return 0;
}

/**************************************************************
    Problem: 3505
    User: RicardoWang
    Language: C++
    Result: Accepted
    Time:344 ms
    Memory:1272 kb
****************************************************************/
 
#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
using namespace std;
#define MAXN 1005
int N,M;
long long ans,t;
int GCD(int x,int y)
{
    int t;
    while(y){t=x%y; x=y; y=t;}
    return x;
}
int main()
{
    scanf("%d%d",&N,&M); N++; M++;
    ans=N*M;
    ans=ans*(ans-1)*(ans-2)/6;
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<M;j++)if(i|j)
        {
            t=GCD(i,j);
            if(t>1)ans-=(long long)(t-1)*(N-i)*(M-j)*( i&&j ? 2 : 1);
        }
    }
    cout<<ans<<endl;
    return 0;
}