鉴于自己每次见到有关组合计数(求方案数之类)的题都要跪,所以专门开个网页记录总结一下。。。不定期更新
【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;
}