数论三之排列组合Ⅱ——Virus Tree 2，RGB Coloring，123 Triangle，排列计数，排队，卡农

最新推荐文章于 2022-12-28 11:53:39 发布
最新推荐文章于 2022-12-28 11:53:39 发布
阅读量228 收藏 1
点赞数 1
分类专栏： 数学 # 排列组合 文章标签： 数学 排列组合
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Emm_Titan/article/details/119221609
版权

丝且人一口

Virus Tree 2

description

solution

距离不超过 2 2 2的点，可能是儿子/孙子/父亲/爷爷

考虑从上到下，由上面的颜色选取决定下面的颜色

显然，第一个点有 K K K种选择，第二个点有 K − 1 K-1 K−1种，每次不同都要 − 1 -1 −1

答案就是每个点的选择的乘积

所以只需要把这种过程的递减通过dfs树来实现

对于现在的子树根节点 u u u，如果 v v v是第一个儿子，那么需要考虑 v v v有没有爷爷

如果不是第一个儿子，那么就是前一个儿子的方案数再 − 1 -1 −1

code

#include <cstdio>
#include <vector>
using namespace std;
#define mod 1000000007
#define int long long
#define maxn 1000005
vector < int > G[maxn];
int n, k;
int w[maxn];

void dfs( int u, int fa ) {
	int r = 0;
	for( auto v : G[u] ) {
		if( v == fa ) continue;
		else {
			if( r ) w[v] = r - 1;
			else if( ! fa ) w[v] = k - 1;
			else w[v] = k - 2;
			r = w[v]; 
		}
	}
	for( auto v : G[u] ) 
		if( v ^ fa ) dfs( v, u );
}

signed main() {
	scanf( "%lld %lld", &n, &k );
	for( int i = 1, u, v;i < n;i ++ ) {
		scanf( "%lld %lld", &u, &v );
		G[u].push_back( v );
		G[v].push_back( u );
	}
	w[1] = k;
	dfs( 1, 0 );
	int ans = 1;
	for( int i = 1;i <= n;i ++ )
		if( w[i] <= 0 ) return ! printf( "0\n" );
		else ans = ans * w[i] % mod;
	printf( "%lld\n", ans );
	return 0;
}

RGB Coloring

description

solution

将绿色 A + B A+B A+B看成既涂了红色 A A A，又涂了蓝色 B B B，则红色和蓝色格子彼此独立，就算涂在同一个格子上也没关系

枚举红色格子的数量 i i i，计算得出蓝色格子数量 j = K − A × i B j=\frac{K-A\times i}{B} j=BK−A×i​

判断格子数量为整数且都在 n n n以内即可，然后用组合计数从 n n n个格子中选红/蓝色

∑ C n i × C n j \sum C_n^i\times C_n^j ∑Cni​×Cnj​

code

#include <cstdio>
#define int long long
#define mod 998244353
#define maxn 300005
int n, A, B, k;
int fac[maxn], inv[maxn];

int qkpow( int x, int y ) {
	int ans = 1;
	while( y ) {
		if( y & 1 ) ans = ans * x % mod;
		x = x * x % mod;
		y >>= 1;
	}
	return ans;
}

void init() {
	fac[0] = inv[0] = 1;
	for( int i = 1;i <= n;i ++ )
		fac[i] = fac[i - 1] * i % mod;
	inv[n] = qkpow( fac[n], mod - 2 );
	for( int i = n - 1;i;i -- )
		inv[i] = inv[i + 1] * ( i + 1 ) % mod;
}

int C( int n, int m ) {
	if( m < 0 || n < m ) return 0;
	else return fac[n] * inv[m] % mod * inv[n - m] % mod;
}

signed main() {
	scanf( "%lld %lld %lld %lld", &n, &A, &B, &k );
	init();
	int ans = 0; 
	for( int i = 0, j;i <= n;i ++ ) {
		if( ( k - A * i ) % B ) continue;
		else j = ( k - A * i ) / B;
		ans = ( ans + C( n, i ) * C( n, j ) ) % mod;
	}
	printf( "%lld\n", ans );
	return 0;
}

123 Triangle

description

solution

N > 1 N>1 N>1式子是差分形式，后面序列只有能是 0 / 1 / 2 0/1/2 0/1/2

（除非 N = 1 N=1 N=1的情况答案有可能是 3 3 3）

如果此时序列有 1 1 1，那么答案一定只能是 0 / 1 0/1 0/1，因为 0 / 2 0/2 0/2碰到 1 1 1都会变成 1 1 1

如果序列没有 1 1 1，则答案只会是 0 / 2 0/2 0/2，对于这种情况，将每个数 / 2 /2 /2得到结果后再 × 2 \times 2 ×2是等价的

所以现在的答案都只会是 0 / 1 0/1 0/1，考虑在 ( m o d 2 ) \pmod 2 (mod2)意义下做

在 ( m o d 2 ) \pmod 2 (mod2)的情况下，加减法可以看成二进制的异或

所求即为一段序列相邻两个数异或直到最后成一个数的答案

考虑每个数被异或的次数，这其实是个倒杨辉三角 ( 0 , 0 ) (0,0) (0,0)开始

长为 n n n的序列的第 i i i个元素被异或的次数为 C n − 1 i − 1 C_{n-1}^{i-1} Cn−1i−1​

计算在模 2 2 2意义下（除法不一定有逆元）的组合数

  • 计算每个阶乘中分解 2 2 2的次数，将除法变成减法，对于组合数最后剩下的 2 2 2的次数为 0 0 0，则为 1 ( m o d 2 ) 1\pmod 2 1(mod2)，反之 0 ( m o d 2 ) 0\pmod 2 0(mod2)
  • 阶乘的分解，拆分成对每个数求出分解中 2 2 2的次数，做前缀和

code

#include <cmath>
#include <cstdio>
#define maxn 1000005
int n;
char s[maxn];
int x[maxn], cnt[maxn]; 

int main() {
	scanf( "%d %s", &n, s + 1 );
	n --;
	for( int i = 1;i <= n;i ++ )
		x[i] = fabs( s[i] - s[i + 1] );
	bool flag = 1;
	for( int i = 1;i <= n;i ++ )
		if( x[i] == 1 ) { flag = 0; break; }
	if( flag ) 
		for( int i = 1;i <= n;i ++ )
			x[i] >>= 1;
	for( int i = 1;i <= n;i ++ ) {
		int t = i;
		while( ! ( t & 1 ) ) t >>= 1, cnt[i] ++;
		cnt[i] += cnt[i - 1];
	}
	int ans = 0;
	for( int i = 1;i <= n;i ++ )
		ans ^= cnt[n - 1] - cnt[i - 1] - cnt[n - i] ? 0 : ( x[i] & 1 );
	if( flag ) ans <<= 1;
	printf( "%d", ans );
	return 0;
}

[SDOI2016]排列计数

description

solution

组合数从 n n n个中选 m m m个强制 a i = i a_i=i ai​=i，剩下的 n − m n-m n−m个 a i ≠ i a_i≠i ai​​=i，就是 n − m n-m n−m的错排数

code

#include <cstdio>
#define mod 1000000007
#define int long long
#define maxn 1000005
int fac[maxn], inv[maxn], D[maxn];
int T, n, m;

int qkpow( int x, int y ) {
	int ans = 1;
	while( y ) {
		if( y & 1 ) ans = ans * x % mod;
		x = x * x % mod;
		y >>= 1;
	}
	return ans;
}

void init() {
	fac[0] = inv[0] = 1;
	for( int i = 1;i < maxn;i ++ )
		fac[i] = fac[i - 1] * i % mod;
	inv[maxn - 1] = qkpow( fac[maxn - 1], mod - 2 );
	for( int i = maxn - 2;i;i -- )
		inv[i] = inv[i + 1] * ( i + 1 ) % mod;
	D[0] = 0, D[1] = 0, D[2] = 1;
	for( int i = 3;i < maxn;i ++ )
		D[i] = ( D[i - 1] + D[i - 2] ) * ( i - 1 ) % mod;
}

int C( int n, int m ) {
	return fac[n] * inv[m] % mod * inv[n - m] % mod;
}

signed main() {
	init();
	scanf( "%lld", &T );
	while( T -- ) {
		scanf( "%lld %lld", &n, &m );
		if( n == m ) printf( "1\n" ); 
		else printf( "%lld\n", C( n, m ) * D[n - m] % mod );
	}
	return 0;
}

[HNOI2012]排队

description

solution

不给模数还明示答案可能很大，赤裸裸的大整数胁迫！！居心叵测！！

两名老师真的是很烦了，没有这么讨人厌的老师，可以直接女生插板法了

  • 两名老师之间是男生

    先排男生 A n n A_n^n Ann​，产生 n + 1 n+1 n+1个空

    再插板两名老师， A n + 1 2 A_{n+1}^2 An+12​，产生 n + 3 n+3 n+3个空

    最后插板 m m m名女生， A n + 3 m A_{n+3}^m An+3m​

  • 两名老师之间是女生

    先排男生 A n n A_n^n Ann​

    再打包两名老师放同一个隔板里，老师间还有顺序， A 2 2 C n + 1 1 A_2^2C_{n+1}^1 A22​Cn+11​

    先强制选一名女生放进老师中间 m m m，剩下的就隔板插

    男生 n n n个，老师和强制女生打包装成一个，总共是 n + 1 n+1 n+1个，产生 n + 2 n+2 n+2个空， A n + 2 m − 1 A_{n+2}^{m-1} An+2m−1​

综上， a n s = A n n A n + 1 2 A n + 3 m + A 2 2 C n + 1 1 A n + 2 m − 1 ans=A_n^nA_{n+1}^2A_{n+3}^m+A_2^2C_{n+1}^1A_{n+2}^{m-1} ans=Ann​An+12​An+3m​+A22​Cn+11​An+2m−1​

同样，老师不相邻 = = =不考虑老师 − - −老师相邻

  • 不考虑老师

    把老师当成男的，女的直接插， A n + 2 n + 2 A n + 3 m A_{n+2}^{n+2}A_{n+3}^m An+2n+2​An+3m​

  • 老师相邻

    捆绑法强制两名老师一起当成一个男的 A 2 2 A n + 1 n + 1 A n + 2 m A_2^2A_{n+1}^{n+1}A_{n+2}^m A22​An+1n+1​An+2m​

二者做差就是答案
在这里插入图片描述

code

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int n, m;

struct Int {
	int g[20000], len;
	
	Int() {
		memset( g, 0, sizeof( g ) );
		len = 0;
	}
	
	Int( int x ) {
		memset( g, 0, sizeof( g ) );
		len = 0;
		if( ! x ) len = 1;
		else while( x ) g[len ++] = x % 10, x /= 10;
	}
	
	void clean() {
		while( len > 1 && ! g[len - 1] ) len --;
	}
	
	Int operator - ( Int t ) {
		Int ans;
		ans.len = len;
		for( int i = 0;i < t.len;i ++ ) {
			ans.g[i] = g[i] - t.g[i];
			if( ans.g[i] < 0 ) ans.g[i] += 10, g[i + 1] --;
		}
		for( int i = t.len;i < len;i ++ )
			ans.g[i] = g[i]; 
		ans.clean();
		return ans;
	}
	
	Int operator * ( Int t ) {
		Int ans;
		ans.len = len + t.len;
		for( int i = 0;i < len;i ++ )
			for( int j = 0;j < t.len;j ++ )
				ans.g[i + j] += g[i] * t.g[j];
		for( int i = 0;i < ans.len;i ++ )
			ans.g[i + 1] += ans.g[i] / 10, ans.g[i] %= 10;
		ans.len ++;
		ans.clean();
		return ans;
	} 
	
	void print() {
		for( int i = len - 1;~ i;i -- )
			printf( "%d", g[i] );
	}
	
};

Int calc( int n, int m ) {
	Int ans = ( 1 );
	for( int i = n - m + 1;i <= n;i ++ )
		ans = ans * i;
	return ans;
}

int main() {
	scanf( "%d %d", &n, &m );
	Int ans;
	ans = calc( n + 2, n + 2 ) * calc( n + 3, m ) - calc( 2, 2 ) * calc( n + 1, n + 1 ) * calc( n + 2, m );	
	ans.print();
	return 0;
}

[HNOI2011]卡农

description

solution

同种音乐其实通过除以 m ! m! m!就可以消掉

相当于在 S S S中选 m m m个子集，满足

  • 子集非空
  • 不存在两个完全一样的子集
  • ∀ i , i ∈ [ 1 , n ] \forall_{i,i\in[1,n]} ∀i,i∈[1,n]​每个元素出现次数为偶数

设 f i : f_i: fi​: 考虑 i i i个子集，满足以上所有性质的方案数

  • 如果知道 i − 1 i-1 i−1个子集的，那么 i i i个子集的集合就被确定了下来 A 2 n − 1 i − 1 A_{2^n-1}^{i-1} A2n−1i−1​
  • 如果子集是空，那么去掉这个子集剩下的 i − 1 i-1 i−1个子集是合法方案 f i − 1 f_{i-1} fi−1​
  • 如果 i i i子集与 j j j子集重复， j j j有 i − 1 i-1 i−1种选择，剔除这两个子集剩下 i − 2 i-2 i−2子集是合法的 f i − 2 f_{i-2} fi−2​， i i i子集选择有 2 n − 1 − ( i − 2 ) 2^n-1-(i-2) 2n−1−(i−2)（排除掉剩下 i − 2 i-2 i−2个子集，不能与之重复）

在这里插入图片描述

最后别忘了 m ! m! m!

code

#include <cstdio>
#define mod 100000007
#define int long long
#define maxn 1000005
int n, m;
int A[maxn], f[maxn];

int qkpow( int x, int y ) {
	int ans = 1;
	while( y ) {
		if( y & 1 ) ans = ans * x % mod;
		x = x * x % mod;
		y >>= 1;
	}
	return ans;
}

signed main() {
	scanf( "%lld %lld", &n, &m );
	f[0] = 1, f[1] = 0;
	int t = ( qkpow( 2, n ) - 1 + mod ) % mod;
	A[0] = 1;
	for( int i = 1;i <= m;i ++ )
		A[i] = ( A[i - 1] * ( t - i + 1 ) % mod + mod ) % mod;
	int MS = 1;
	for( int i = 2;i <= m;i ++ ) {
		f[i] = ( A[i - 1] - f[i - 1] - f[i - 2] * ( i - 1 ) % mod * ( t - i + 2 ) % mod ) % mod;
		MS = MS * i % mod;
	}
	printf( "%lld\n", ( f[m] + mod ) % mod * qkpow( MS, mod - 2 ) % mod );
	return 0;
}
ikrvxt
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
Virus Tree 2（树）
eazo的博客
02-06 482
题目描述 You are given a tree with N vertices and N−1 edges. The vertices are numbered 1 to N, and the i-th edge connects Vertex ai and bi. You have coloring materials of K colors. For each vertex in the ...
C++基础数论————排列组合
C2020lax的博客
02-28 3798
排列组合的定义： 排列就是从n个数中找出m个进行排列，共有多少种方案？ 这一定要考虑顺序。 公式就是： 而组合就从从n个数中找出m个进行组合，共有多少种方案？ 这一定不要考虑顺序。 公式就是： 现在，我们来看一看它的性质： 而我们该如何计算组合数呢？ 计算组合数： 这一共有2种方法。 1、根据定义直接暴力求出组合数，但是其时间复杂度较长而且有溢出的风险，所以不推荐...
openwrt 水星mw4530r-v1 搞搞搞
weixin_33713503的博客
02-12 323
感觉周围最实用的搞硬件非路由器莫属，可惜配置都不咋的高，选择水星这款就是看中它的性价比和openwrt的支持，真乃刷机神器啊，还可以挂载usb就更绝了，价格没得说。 另外128mb的rom对大部分功能足够了，不行还可以在挂载的u盘上可以扩展，绝对是刷机的首选啊。 openwrt入门 没玩不知道，比树莓派更简单，直接上传bin文件就可以了，感谢openwrt。上传有很多种办法，最直接的就是在更...
数论和组合数学经典的组合习题解
01-15
组合数学中较为经典的题目,特别是对考试而言.一定要学会这几道题的计算!
AtCoder Beginner Contest 133 E - Virus Tree 2 （dfs）
pxlsdz的博客
07-14 274
E - Virus Tree 2 Time Limit: 2 sec / Memory Limit: 1024 MB Score :500500points Problem Statement You are given a tree withNNvertices andN−1N−1edges. The vertices are numbered11toNN, an...
错位排列——数论
JKdd123456的博客
07-24 4703
首先先介绍一下什么是错位排列： 错位重排是指一种比较难理解的复杂数学模型，是伯努利和欧拉在错装信封时发现的，因此又称伯努利-欧拉装错信封问题 也是离散课本上出现的戴错帽子问题 简介 表述为：编号是1、2、…、n的n封信，装入编号为1、2、…、n的n个信封，要求每封信和信封的编号不同，问有多少种装法？ 对这类问题有个   固定的递推公式  （类似于 dp、递推之类的），记n封信的错位重排数为...
第三十二章 数论——组合数详解（1）
Turing_Sheep的博客
12-28 2464
详细讲解不同情况下组合数的快速求法，超级详细地证明原理，并帮你熟练掌握精简代码模板。
用Python3实现排列组合
COCO56的博客
04-24 3535
文章目录认识排列与组合用代码实现考虑顺序的排列用代码实现不考虑顺序的组合 认识排列与组合 排列组合是组合学最基本的概念。所谓排列（Permutation），就是指从给定个数的元素中取出指定个数的元素进行排序。组合（Combination）则是指从给定个数的元素中仅仅取出指定个数的元素，不考虑排序。 比如我们要实现1，2，3的排列组合，我们可以很容易写出来，如下表： 任意选取的方式 待取的元素 任意选取1个元素 任意选取2个元素 任意选取3个元素 考虑顺序（排列） 123 1,2,3 12,2
排列组合问题
蒟蒻的算法浅析
11-25 2166
排列组合问题 ​ 排列组合是组合数学里的两大经典问题，下面我们先来看一下它的定义： ​ 排列：从n个不同元素中任取m（m≤n）个元素，按照一定的顺序排列起来，叫做从n个不同元素中取出m个元素的一个排列。 ​ 组合：一般地,从n个不同的元素中,任取m(m≤n)个元素为一组,（即不考虑顺序），叫作从n个不同元素中取出m个元素的一个组合。 一、全排列问题 ​ 首先我们可以先简化一下这个问题，就令n=mn...
初等数论习题解答第三版
10-14
初等数论习题解答第三版
数论（组合数求法+Lucas定理+线性求逆元
08-01
数论
学习笔记——数论——约数
01-06
2.一个整数N的约数个数的上界为. 3.1-N 每个数的约数个数综合大约为N * log N 个。 4.1-1e9中的数，约数个数最多的数仅有1536个约数 倍数法求1-N中所有数的约数集合：NlogN void gao() { vectorfactor[M]; for...
初等数论系列之三
01-18
初等数论是计算机编程和应用数学的技术科技课程，对于人工智能也很有帮助
NTT.rar_NTT_NTT c++_NTT%20c%2B%2B_nttc++_快速数论变换
09-20
用c++实现的快速数论变换(NTT)，利用原根性质，实现和FFT的功能，但是效率比FFT高。
ABC133 E Virus Tree 2 树形遍历+染色
ruanluyu的博客
12-08 207
题目：E - Virus Tree 2 给定一个含KKK个节点的无向树，给每个顶点赋上[1,K][1,K][1,K]一整数（染上一种颜色）。并且，若两个顶点距离≤2\le2≤2，那么两顶点值不同（颜色不同）。 1≤N,K≤1051\le N,K\le 10^51≤N,K≤105 考察 首先无向树的根可以随意设定，我们就从第1个开始查找（邻接表）。根的颜色肯定是K种可用。我们固定好根的颜色后，用D...
【DFS算方案数】AtCoder Beginner Contest 133 E - Virus Tree 2
qq_41037114的博客
07-08 350
https://atcoder.jp/contests/abc133/tasks/abc133_e E - Virus Tree 2 Time Limit: 2 sec / Memory Limit: 1024 MB Score :500500points Problem Statement You are given a tree withNNvertices andN...
【无码专区2】序列划分（数学）
热门推荐
10-19 9万+
天啦，只看20位！！
【无码专区11】异或2（结论 / 推式子 + 哈希hash + 大整数高精度 加减乘除重载考察）
11-13 8万+
直接写吐了，真的nm难调啊
三数论大小，用指针做函数参数实现（指针）
最新发布
06-09
这里使用了三个指向 int 类型的指针作为函数参数，通过比较指针所指向的值的大小来实现三个数的排序。在主函数中，通过使用 `&` 取地址符获取变量的地址，再将地址传递给函数，实现了指针作为函数参数的功能。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值