5164. 【NOIP2017模拟6.25】小A做CF

最新推荐文章于 2019-07-08 14:19:00 发布

lazy-sheep

最新推荐文章于 2019-07-08 14:19:00 发布

阅读量237

点赞数 1

分类专栏：题目 # 高精度文章标签：高精度思维

本文链接：https://blog.csdn.net/zsjzliziyang/article/details/82968056

版权

题目同时被 2 个专栏收录

246 篇文章 0 订阅

订阅专栏

高精度

12 篇文章 0 订阅

订阅专栏

Description

目标rating超过CLJ的小A最近在疯狂做CF，话说这次CF在赛前十分钟出了个相当奇葩的预选题，要求十分钟内必须做出这道题才能参加这次CF。
这道题是这样的，给你一个N*N的矩阵，每行有一个障碍，数据保证任意两个障碍不在同一行，任意两个障碍不在同一列，要求你在这个矩阵上放N枚棋子（障碍的位置不能放棋子），要求你放N个棋子也满足每行只有一枚棋子，每列只有一枚棋子的限制，求有多少种方案。
最近被各种神题折磨得死去活来，走火入魔的小A一下被这个超级水题给卡住了，但是这次是超过CLJ的最佳机会，于是，心急火燎的他找到了你，请你带他顺利参加这次CF正赛吧。

Input

第一行一个N，接下来一个N*N的矩阵。

Output

合法方案数。

Sample Input

1 0

0 1

Sample Output

Data Constraint

20%的数据保证： N<=10
60%的数据保证： N<=20
100%的数据保证： N<=200

Solution

20%爆搜+剪枝

60%爆搜+更好的剪枝来骗分（本验题人未打过这档分，所以拿不到的也没办法）

100%有两种方法，一种是原题解方法，一种是本验题人方法。

法一：原题解本题目由于障碍其实是顺序无关的（从每行每列不重复可看出），于是可以将棋盘进行转化： 1 0 0 1 0 0 1 0 0 0 0 1 → 0 1 0 → 0 1 + 0 0 0 0 0 0 然后可以通过求空白棋盘和有障碍棋盘的方案数，再相乘可得。如果将列位置看成数列的位置，行位置看做数列中的数，显然可以看出棋盘摆棋情况可映射为一个数列，于是有：空白棋盘方案：（全排列）n! 带障碍棋盘方案：用错排公式（就是由 n 个人坐座位，然后进行换位，每个人不能坐在原位置的方案数），递推式： f[i]=(i-1)*(f[i-1]+f[i-2])

法二：验题人解法基于转化后的模型，可以发现，对于空白棋盘方案数是相当容易求出来的，而带障碍棋盘方案可以通过容斥原理实现（即将总的减去放一个棋子在障碍上再加回两个棋子在障碍上……）预处理出组合数（因为放障碍时要选位置），可以得出如下式子：（设有 n 行 n 列，m 个障碍） Ans=∑(-1) i*C i m*（n-i）！（i 枚举障碍数，由 0 到 m）

Code1

#include<cstdio>
using namespace std;
int n,k,a[1005];
void cheng(int x)
{
	int t=0;
	for (int i=1;i<=a[0];i++)
	{
		a[i]=a[i]*x+t;
		t=a[i]/10;
		a[i]%=10;
	}
	while (t>0)
	{
		a[++a[0]]=t;
		t=a[a[0]]/10;
		a[a[0]]%=10;
	}
}
void jia()
{
	a[1]++;
	int t=1;
	while (a[t]==10)
	{
		a[t+1]++;
		a[t]=0;
		t++;
	}
	if (t>a[0]) a[0]=t;
}
void jian()
{
	int t=1;
	while (a[t]==0 && t<a[0]) t++;
	a[t]--;
	for (int i=t-1;i;i--) a[i]=9;
	if (t==a[0] && a[t]==0) a[0]--;
}
int main()
{
	scanf("%d",&n);
	a[0]=1;
	a[1]=0;
	k=-1;
	for (int i=2;i<=n;i++)
	{
		k*=-1;
		cheng(i);
		if (k>0) jia(); else jian();
	}
	for (int i=a[0];i;i--) printf("%d",a[i]);
	return 0;
}

Code2

#include<cstdio> 
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
#define mod 1000000000000000
using namespace std;
ll n,f[210][1000];
void update(int x){
	for(int i=1;i<=f[x][0];i++){
		f[x][i+1]+=f[x][i]/mod;
		f[x][i]%=mod;
	}
	while(f[x][f[x][0]+1]>0){
		f[x][++f[x][0]]+=f[x][f[x][0]-1]/mod;
		f[x][f[x][0]-1]%=mod;
	}
}
int getlen(int x){
	ll k=f[n][x],len=0;
	while(k>0){k/=10;len++;}
	return len;
}
int main(){
	scanf("%lld",&n);
	f[2][0]=1;f[2][1]=1;
	f[3][0]=1;f[3][1]=2;
	for(int i=4;i<=n;i++){
		f[i][0]=f[i-1][0];
		for(int j=1;j<=f[i][0];j++) f[i][j]+=f[i-1][j]+f[i-2][j];
		update(i);
		for(int j=1;j<=f[i][0];j++) f[i][j]*=(i-1);
		update(i);
	}
	printf("%lld",f[n][f[n][0]]);
	for(int i=f[n][0]-1;i>0;i--){
		for(int j=1;j<=15-getlen(i);j++) printf("0");
		printf("%lld",f[n][i]);
	}
	return 0;
}

作者：zsjzliziyang
QQ:1634151125
转载及修改请注明
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/82968056