[CSP-S模拟测试]:幻魔皇(数学)

题目描述

  幻魔皇拉比艾尔很喜欢斐波那契树,他想找到神奇的节点对。
  所谓斐波那契树,根是一个白色节点,每个白色节点都有一个黑色节点儿子,而每个黑色节点则有一个白色和一个黑色节点儿子。神奇的节点对则是指白色节点对。
  请问对于深度为$n$的斐波那契树,其中距离为$i$的神奇节点对有多少个?拉比艾尔需要你对于$1\leqslant i\leqslant 2n$的所有$i$都求出答案。


输入格式

一行一个正整数$n$。


输出格式

一行$2n$个整数表示答案,对$123456789$取模。


样例

样例输入:

5

样例输出:

0 2 3 3 1 1 1 0 0 0 0


数据范围与提示

对于$20\%$的数据$n\leqslant 10$;
对于$40\%$的数据$n\leqslant 20$;
对于$60\%$的数据$n\leqslant 30$;
对于$80\%$的数据$n\leqslant 400$;
对于$100\%$的数据$n\leqslant 5,000$。


题解

先来看一下这棵树:

首先,来看几点性质:

  $\alpha.$每一行的点数呈斐波那契数列。

  $\beta.$从第二行开始,每行的黑点个数呈斐波那契数列。

  $\gamma.$从第三行开始,每行的白点个数呈斐波那契数列。

  $\delta.$对于每一个黑(白)点,如果它有一棵深度为$k$的子树,那么所有这些子树都是一样的。

得到了这些性质,再来看这道题:

设$w[i]$表示第$i$行的白点的个数,$b[i]$表示第$i$行黑点的个数,$sw[i]$表示小于等于$i$行的所有白点的个数,$sb[i]$表示小于等于$i$行的所有黑点的个数(注意此处的行数其实也是深度,根节点的深度为$1$)。

思考如何计数,根据性质$\delta$,可以考虑将每一棵不同的子树内部贡献,再乘上这样子树的数量即可得到答案。

再来考虑如何统计子树内部的贡献,可以大致分为两种情况:

  $\alpha.$神奇节点对的$LCA$为白色点。

  $\beta.$神奇节点对的$LCA$为黑色点。

观察上面的图可知,对于情况$\alpha$,$LCA$一定是其中一个白点,就是说这两个白点一定是祖先和后代的关系。

那么,对于一个深度为$k$的子树,其贡献为$w[i+1]$(因为祖先白点的行数是$1$),而这样的子树有$sw[n-i]$棵。

所以,对于情况$\alpha$,其对答案的贡献即为:

$$\sum \limits_{i=1}^n sw[n-i]\times w[i+1]$$

对于情况$\beta$,就显得麻烦很多了,两个白点必须在这个黑点($LCA$)的两个子节点的子树中。

考虑枚举两侧的距这个黑点的距离$i,j$(白色子节点的子树中的白点的距离为$i$,黑色子节点的子树中的白点的距离为$j$),并将这棵树进行如下图所示划分:

惊喜的发现,如果将右侧树的根节点由黑变白对答案没有影响,而右侧树就变成了一棵深度为$j+1$的普通树,那么其距离为$j$的白点就有$w[j+1]$个。

左侧深度为$i+1$,则对于左侧树其深度为$i$,所以距离为$i$的白点个数为$w[i]$。

那么对于这棵树,其对答案的贡献就是$w[i]\times w[j+1]$,而这样的子树的数量取决与$i,j$中更大的那一个,所以这样的子树有$sb[n-\max(i,j)]$棵。

所以,对于情况$\beta$,其对答案的贡献即为:

$$\sum \limits_{i=1}^n \sum \limits_{j=1}^n sb[n-\max(i,j)]\times w[i]\times w[j+1]$$

时间复杂度:$\Theta(n^2)$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
using namespace std;
const int mod=123456789;
int n;
long long w[5001],b[5001],sw[5001],sb[5001];
long long ans[10001];
int main()
{
	scanf("%d",&n);
	w[1]=w[3]=1;
	sw[1]=sw[2]=1;sw[3]=2;
	for(int i=4;i<=n;i++)
	{
		w[i]=(w[i-1]+w[i-2])%mod;
		sw[i]=(sw[i-1]+w[i])%mod;
	}
	b[2]=1;sb[2]=1;
	for(int i=3;i<=n;i++)
	{
		b[i]=(b[i-1]+b[i-2])%mod;
		sb[i]=(sb[i-1]+b[i])%mod;
	}
	for(int i=1;i<=n;i++)
		ans[i]=sw[n-i]*w[i+1]%mod;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			ans[i+j]=(ans[i+j]+sb[n-max(i,j)]*w[i]%mod*w[j+1])%mod;
	for(int i=1;i<=2*n;i++)printf("%lld ",ans[i]);
	return 0;
}


rp++

原文地址:https://www.cnblogs.com/wzc521/p/11658504.html

时间: 2024-10-05 21:05:25

[CSP-S模拟测试]:幻魔皇(数学)的相关文章

[CSP-S模拟测试]:Race(数学+Trie树)

题目描述 一年一度的运动会开始了.有$N$个选手参赛,第$i$个选手有一个能力值(保证$A[i]$两两不同),比赛一共进行了天.在第$j$天($0\leqslant j\leqslant 2^{m-1}$)的比赛中,第$i$个选手的得分为$A[i]\ xor\ j$,然后从大到小排名,排名为$x$($x$从$0$开始)的同学会获得的积分,你需要求出每个同学最后总的积分和$q[i]$模${10}^9+7$的结果$p[i]$.为了避免输出文件过大,你只要输出$p[i]$的异或和即可. 输入格式 第一

[CSP-S模拟测试]:Equation(数学+树状数组)

题目描述 有一棵$n$个点的以$1$为根的树,以及$n$个整数变量$x_i$.树上$i$的父亲是$f_i$,每条边$(i,f_i)$有一个权值$w_i$,表示一个方程$x_i+x_{f_i}=w_i$,这$n−1$个方程构成了一个方程组.现在给出$q$个操作,有两种类型:$\bullet 1\ u\ v\ s$,表示询问加上$x_u+x_v=s$这个方程后,整个方程组的解的情况.具体来说,如果方程有唯一解,输出此时$x_1$的值:如果有无限多个解,输出$inf$:如果无解,输出$none$.注意

[CSP-S模拟测试]:Divisors(数学)

题目描述 给定$m$个不同的正整数$a_1,a_2,...,a_m$,请对$0$到$m$每一个$k$计算,在区间$[1,n]$里有多少正整数是$a$中恰好$k$个数的约数. 输入格式 第一行包含两个正整数$n,m$,分别表示区间范围以及$a$数组的大小.第二行包含$m$个不同的正整数$a_1,a_2,...,a_m$,表示$a$数组. 输出格式 输出$m+1$行,每行一个整数,其中第$i$行输出$k=i$的答案. 样例 样例输入1: 10 34 6 7 样例输出1: 4411 样例输入2: 5

[CSP-S模拟测试]:Silhouette(数学)

题目描述 有一个$n\times n$的网格,在每个格子上堆叠了一些边长为$1$的立方体. 现在给出这个三维几何体的正视图和左视图,求有多少种与之符合的堆叠立方体的方案.两种方案被认为是不同的,当且仅当某个格子上立方体的数量不同. 输出答案对$10^9+7$取模的结果. 输入格式 从文件$silhouette.in$中读入数据. 第一行一个整数$n$. 第二行$n$个整数,第$i$个表示正视图中从左到右第$i$个位置的高度$A_i$. 第三行$n$个整数,第$i$个表示左视图中从左到右第$i$个

[CSP-S模拟测试]:sum(数学+莫队)

题目传送门(内部题63) 输入格式 第一行有一个整数$id$,表示测试点编号.第一行有一个整数$q$,表示询问组数.然后有$q$行,每行有两个整数$n_i,m_i$. 输出格式 一共有$q$行,每行一个整数表示每组询问的答案$S_{n_i,m_i}$对$10^9+7$取模的结果. 样例 样例输入: 151 12 13 24 35 5 样例输出: 2371532 数据范围与提示 对于所有数据,$1\leqslant q,n_i,m_i\leqslant 10^5$. 题解 考场上把$80$分部分分

[CSP-S模拟测试]:棋盘(数学+高精度)

题目描述 在一个大小为$N\times N$的棋盘上,放置了$N$个黑色的棋子.并且,对于棋盘的每一行和每一列,有且只有一个棋子.现在,你的任务是再往棋盘上放置$N$个白色的棋子.显然,白色棋子不能与黑色棋子重合.在此基础上,放置的方式还需要满足:对于棋盘的每一行和每一列,有且只有一个白色棋子.当然,放置的方式有很多种,你只需要输出不同的放置方案数即可. 输入格式 输入文件为$board.in$.第一行包含一个正整数$N$.接下来$N$行,每行$N$个整数用于描述棋盘.$0$表示这个位置是空的,

[CSP-S模拟测试]:chess(数学)

题目描述 $dirty$在一个棋盘上放起了棋子. 棋盘规格为$n\times m$,他希望任意一个$n\times n$的区域内都有$C$个棋子.$dirty$很快就放置好了一个满足条件的棋盘方案,但是他认为这样过于简单了,他希望知道有多少个满足条件的方案. 输入格式 输入三个整数$n,m,C$,含义如题所述. 输出格式 输出一行一个整数,表示答案对$10^9+7$取模的结果. 样例 样例输入: 2 3 1 样例输出: 6 数据范围与提示 对于$20\%$的数据,$n,K\leqslant 4$

[CSP-S模拟测试]:Walker(数学)

题目传送门(内部题86) 输入格式 第一行$n$接下来$n$行,每行四个浮点数,分别表示变换前的坐标和变换后的坐标 输出格式 第一行浮点数$\theta$以弧度制表示第二行浮点数$scale$第三行两个浮点数$d_x,d_y$我们将用$SPJ$以$10^{-3}$的绝对误差来判断变换结果是否正确,建议输出$10$位小数以上. 样例 样例输入1: 50 0 -1 10 1 -2 11 0 -1 21 1 0 02 1 1 0 样例输出1: 1.57079632681-1 1 样例输入2: 50 0

[CSP-S模拟测试]:数列(数学)

题目传送门(内部题95) 输入格式 第一行三个整数$n,a,b$,第二行$n$个整数$x_1\sim x_n$表示数列. 输出格式 一行一个整数表示答案.无解输出$-1$. 样例 样例输入:2 2 31 2 样例输出: 3 数据范围与提示 对于$10\%$的数据,$n,a,b,|x_i|\leqslant 1,000$. 对于$30\%$的数据,$n,a,b\leqslant 1,000$. 对于另外$10\%$的数据,$a=1$. 对于另外$10\%$的数据,$a=2,b=3$. 对于$100