A.
解
对于一个节点 \(x\) ,如果 \(x\) 的 \(2^k\) 级祖先存在,那么把 \(fa[x][k]\) 的答案-1。
Code
#include<bits/stdc++.h>
using namespace std;
const int maxn=100003,maxlog=19;
vector<int> g[maxn];
int n,k,sz[maxn],ans[maxn],fa[maxn][maxlog],dep[maxn];
void init(int u,int last){
fa[u][0]=last;
for(int i=1;(1<<i)<=dep[u];i++)fa[u][i]=fa[fa[u][i-1]][i-1];
sz[u]=1;
for(int i=0;i<int(g[u].size());i++){
int v=g[u][i];
if(v==last)continue;
dep[v]=dep[u]+1;
init(v,u);
sz[u]+=sz[v];
}
}
int Fa(int u,int k){
for(int i=maxlog-1;i>=0;i--){
if((k>>i)&1)u=fa[u][i];
}
return u;
}
int main(){
scanf("%d%d",&n,&k);
for(int i=1;i<n;i++){
int u,v;
scanf("%d%d",&u,&v);
g[u].push_back(v),g[v].push_back(u);
}
init(1,0);
for(int i=1;i<=n;i++)ans[i]=sz[i]-1;
for(int i=1;i<=n;i++){
int x=Fa(i,k+1);
if(x)ans[x]-=sz[i];
}
for(int i=1;i<=n;i++)printf("%d\n",ans[i]);
return 0;
}
B.
解 \(\text{40pts}\)
乱做一通背包。复杂度 \(O(n^2k)\) 。
Code \(\text{40pts}\)
#include<bits/stdc++.h>
using namespace std;
const int mod=1000000009;
int Plus(int x,int y){return (x+=y)>=mod?x-mod:x;}
void PlusEqual(int &x,int y){if((x+=y)>=mod)x-=mod;}
int mul(long long x,int y){return x*y%mod;}
const int maxn=13,maxm=100003;
int n,m,dp[maxn][maxm];
int main(){
while(scanf("%d%d",&m,&n),n||m){
for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)dp[i][j]=0;
dp[0][0]=1;
for(int k=1;k<=m;k++){
for(int i=n;i>=1;i--){
for(int j=k;j<=m;j++){
PlusEqual(dp[i][j],mul(dp[i-1][j-k],4));
if(i>=2&&j>=k*2)PlusEqual(dp[i][j],mul(dp[i-2][j-k*2],6));
if(i>=3&&j>=k*3)PlusEqual(dp[i][j],mul(dp[i-3][j-k*3],4));
if(i>=4&&j>=k*4)PlusEqual(dp[i][j],dp[i-4][j-k*4]);
}
}
}
int ans=0;
for(int i=0;i<=n;i++)PlusEqual(ans,dp[i][m]);
printf("%d\n",ans);
}
return 0;
}
C.
解
\(\huge{\% 10^8+7}\)
考场上打表 \(\text{10pts}\) 。
如果你写对了打表程序,并输出了所有方案,你会发现对于固定的 \(n\) ,只有两种方案是本质不同的。
第一种是把整个正方形拆成几个小正方形。
第二种如图。
然后问题就转化为对于整数 \(n\) ,把其拆成若干个 \(\geq 2\) 的正整数的本质不同的方案数。
或者可以这么考虑:建图,若某一行的两个X位于列 \(i\) 、列 \(j\) ,连边 \((i,j)\) 。最后你会发现对于所有合法方案,这个图都由几个不相邻的简单环组成,且对于所有本质相同的方案,图的形态是一样的。这样更容易得出上述结论。
所以直接 \(O(n^2)\;\text{dp}\) 。
Code
#include<bits/stdc++.h>
using namespace std;
const int mod=100000007;
int Plus(int x,int y){return (x+=y)>=mod?x-mod:x;}
void PlusEqual(int &x,int y){if((x+=y)>=mod)x-=mod;}
int mul(long long x,int y){return x*y%mod;}
const int maxn=2001;
int dp[maxn][maxn];
int main(){
dp[0][0]=1;
for(int i=1;i<maxn;i++){
for(int j=1;j<=i;j++){
PlusEqual(dp[i][j],Plus(dp[i-j][j],dp[i-1][j-1]));
}
}
int T;
scanf("%d",&T);
while(T--){
int n;
scanf("%d",&n);
int ans=0;
for(int i=1;i<=n;i++)PlusEqual(ans,dp[n-i][i]);
printf("%d\n",ans);
}
return 0;
}