JZYZOJ1544 [haoi2016T2]放棋子 错排公式 组合数学 高精度
http://172.20.6.3/Problem_Show.asp?ID=1544&a=ProbNF
看了题解才意识到原题有错排的性质(开始根本不知道错排是什么)。
十本不同的书放在书架上。现重新摆放,使每本书都不在原来放的位置。有几种摆法?
这个问题的解答就是10个元素的错排数。
错排定义:考虑一个有n个元素的排列,若一个排列中所有的元素都不在自己原来的位置上,那么这样的排列就称为原排列的一个错排。
错排数的递推式
当n个编号元素放在n个编号位置,元素编号与位置编号各不对应的方法数(错排数)用D(n)表示,
那么D(n) = (n-1) *[D(n-2) + D(n-1)];
还是数学方面的积累太少了
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
using namespace std;
const int maxn=;
struct har{
int e[maxn*],len;
har(){memset(e,,sizeof(e));len=;}
};har f[maxn];
int n;
har plu(har x,har y){
har z;z.len=max(x.len,y.len);
for(int i=;i<=z.len;i++){
z.e[i]+=x.e[i]+y.e[i];
z.e[i+]+=z.e[i]/;
z.e[i]%=;
if(z.e[i+]>&&z.len<i+){
z.len=i+;
}
}
return z;
}
har mul(har x,int k){
har z;z.len=x.len;
for(int i=;i<=z.len;i++){
z.e[i]+=x.e[i]*k;
z.e[i+]+=z.e[i]/;
z.e[i]%=;
if(z.e[i+]>&&z.len<i+){
z.len=i+;
}
}
return z;
}
void put(har x){
for(int i=x.len;i>=;i--){
printf("%d",x.e[i]);
}cout<<endl;
}
int main(){
scanf("%d",&n);int x;
for(int i=;i<=n;i++){
for(int j=;j<=n;j++){
scanf("%d",&x);
}
}
f[].e[]=;
f[].e[]=;
for(int i=;i<=n;i++){
f[i]=mul(plu(f[i-],f[i-]),i-);
}
put(f[n]);
return ;
}
JZYZOJ1544 [haoi2016T2]放棋子 错排公式 组合数学 高精度的更多相关文章
- 【BZOJ4563】[Haoi2016]放棋子 错排+高精度
		[BZOJ4563][Haoi2016]放棋子 Description 给你一个N*N的矩阵,每行有一个障碍,数据保证任意两个障碍不在同一行,任意两个障碍不在同一列,要求你在这个矩阵上放N枚棋子(障碍 ... 
- bzoj4563: [Haoi2016]放棋子(错排+高精)
		4563: [Haoi2016]放棋子 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 387 Solved: 247[Submit][Status] ... 
- hdu1465(放错信 错排公式)
		Problem Description 大家常常感慨,要做好一件事情真的不容易,确实,失败比成功容易多了! 做好“一件”事情尚且不易,若想永远成功而总从不失败,那更是难上加难了,就像花钱总是比挣钱容易 ... 
- BZOJ4517:[SDOI2016]排列计数(组合数学,错排公式)
		Description 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是 ... 
- HDU 2048:神、上帝以及老天爷(错排公式,递推)
		神.上帝以及老天爷 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total ... 
- HDU   1465   不容易系列之一   (错排公式+容斥)
		题目链接 Problem Description 大家常常感慨,要做好一件事情真的不容易,确实,失败比成功容易多了! 做好"一件"事情尚且不易,若想永远成功而总从不失败,那更是难上 ... 
- hdu 4535(排列组合之错排公式)
		吉哥系列故事——礼尚往来 Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others)Tota ... 
- HDU 2068 RPG的错排(错排公式 + 具体解释)
		RPG的错排 Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Sub ... 
- HDU 2068 RPG错排 [错排公式]
		1.题意:1到N的序列的排列中,元素位置与元素值相对应的情况(值为i的元素在某个排列中正好排在第i个位置)大于等于序列规模一半的情况,有多少个? 2.输入输出:每组数据一个数,N,规定输入以0结尾: ... 
随机推荐
- 【BZOJ】1709: [Usaco2007 Oct]Super Paintball超级弹珠
			[算法]模拟 [题解]O(n^2)预处理横线(y),纵线(x),主对角线(y-x+n),副对角线(x+y). 然后n^2枚举每个点. 
- 【BZOJ】1014 [JSOI2008]火星人prefix
			[算法]splay [题解]对于每个结点维护其子树串的hash值,前面为高位,后面为低位. sum[x]=sum[L]*base[s[R]+1]+A[x]*base[s[R]]+sum[R],其中su ... 
- NSObject class和NSObject protocol的关系(抽象基类与协议)
			[转载请注明出处] 1.接口的实现 对于接口这一概念的支持,不同语言的实现形式不同.Java中,由于不支持多重继承,因此提供了一个Interface关键词.而在C++中,通常是通过定义抽象基类的方式来 ... 
- dot.js使用心得
			一.dot.js介绍 最近用到的数据模板引擎有很多,今天讲的doT.js也是其中一种. doT.js的特点是体积小,速度快,并且不依赖其他插件. 官网下载:http://olado.github.io ... 
- 工具_HBuilder使用快捷方式
			HBuilder常用快捷键大概共9类([4 13 3]文件.编辑.插入:[4 9 8]选择.跳转.查找:[1 1 6]运行.工具.视图) 1.文件(4) 新建 Ctrl + N 关闭 Ctrl + F ... 
- 上传漏洞新姿势(限Linux)
			服务器:Linux当前环境:nginx/1.4.7PHP版本:PHP Version 7.0.0 上传情况简介:上传 111.jpg111 确实可以成功的但是上传 1.php.jpg1111.1 ... 
- skb管理函数之skb_clone、pskb_copy、skb_copy
			skb_clone--只复制skb描述符本身,如果只修改skb描述符则使用该函数克隆: pskb_copy--复制skb描述符+线性数据区域(包括skb_shared_info),如果需要修改描述符以 ... 
- MinnowBoard
			MinnowBoard https://github.com/RafaelRMachado/MinnowBoard https://github.com/RafaelRMachado https:// ... 
- HTML5API(3)
			十一.ajax同源操作 URL说明是否允许通信 同一域名下允许 http://www.a.com/a.js , http://www.a.com/b.js 同一域名下不同文件夹允许 http://ww ... 
- 12-5 NSSet
			原文:http://rypress.com/tutorials/objective-c/data-types/nsset NSSet NSSet, NSArray, and NSDictionary ... 
