[Luogu P3330] [BZOJ 2227] [ZJOI2011]看电影

最新推荐文章于 2019-08-23 12:09:14 发布

LPA20020220

最新推荐文章于 2019-08-23 12:09:14 发布

阅读量213

点赞数

分类专栏：数学排列组合文章标签：数学排列组合

本文链接：https://blog.csdn.net/LPA20020220/article/details/84964044

版权

数学同时被 2 个专栏收录

122 篇文章 0 订阅

订阅专栏

排列组合

22 篇文章 0 订阅

订阅专栏

洛谷传送门

BZOJ传送门

题目描述

到了难得的假期，小白班上组织大家去看电影。但由于假期里看电影的人太多，很难做到让全班看上同一场电影，最后大家在一个偏僻的小胡同里找到了一家电影院。但这家电影院分配座位的方式很特殊，具体方式如下：

电影院的座位共有 $K$ 个，并被标号为 $1 . . . K$ ，每个人买完票后会被随机指定一个座位，具体来说是从 $1 . . . K$ 中等可能的随机选取一个正整数，设其为 $L$ 。
如果编号 $L$ 的座位是空位，则这个座位就分配给此人，否则将 $L$ 加一，继续前面的步骤。
如果在第二步中不存在编号 $L$ 的座位，则该人只能站着看电影，即所谓的站票。

小白班上共有 $N$ 人（包括小白自己），作为数学爱好者，小白想知道全班都能够有座位的概率是多少。

输入输出格式

输入格式：

输入文件第一行有且只有一个正整数 $T$ ，表示测试数据的组数。第 $2\sim T+1$ 行，每行两个正整数 $N, K$ ，用单个空格隔开，其含义同题目描述。

输出格式：

输出文件共包含 $T$ 行。第 $i$ 行应包含两个用空格隔开的整数 $A, B$ ，表示输入文件中的第 $i$ 组数据的答案为 $A / B$ 。（注意，这里要求将答案化为既约分数）

输入输出样例

输入样例#1：

输出样例#1：

1 1
0 1
3 4

数据范围

对于100%的数据， $T\le 50,\ n,k\le 200$

解题分析

神仙题，感觉思路和THU的你的生命已如风中残烛差不多，然而还是不会做…

考虑在最后一个位置的后面加一个位置，并让超出这个位子而站着的人回到 $1$ 号位置，这样只要 $k\ge n$ 就一定有解。

这样形成了一个环，所有人都可以选 $1\to k+1$ 的所有位置，所以方案是 $k+1)^n$ ，然而事实上是存在以 $1, 2, 3, . . ., k + 1$ 开头的同构的情况的，所以需要除以 $k + 1$ 。

那么合法方案当且仅当 $k + 1$ 号位置不被占，所以我们直接枚举空的位置变成 $k + 1$ 号位置即可。

最后答案就是 $\frac{(k+1)^{n-1}\times (k+1-n)}{k^n}$ 。

注意 $\bot k$ ，所以 $g c d$ 算 $k + 1 - n$ 和 $k^n$ 的就好了，否则可能会 $T$ …

不得不吐槽这数据真水，只有一组数据有 $0\ 1$ 的情况（如果想看蒟蒻博主怎么知道的请看下一篇博客）

代码如下：

#include <cstdio>
#include <cmath>
#include <cctype>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#define R register
#define IN inline
#define W while
using namespace std;
namespace BigInteger
{
#define maxn 1005

struct Big_integer{  
	int d[maxn], len;  

	void clean() { while(len > 1 && !d[len-1]) len--; }  

	IN Big_integer()          { memset(d, 0, sizeof(d)); len = 1; }  
	IN Big_integer(int num)   { *this = num; }   
	IN Big_integer(char* num) { *this = num; }  
	IN Big_integer operator = (const char* num){  
		memset(d, 0, sizeof(d)); len = strlen(num);  
		for(R int i = 0; i < len; i++) d[i] = num[len-1-i] - '0';  
		clean();  
		return *this;  
	}  

	IN Big_integer operator = (int num){  
		char s[10005]; sprintf(s, "%d", num);  
		*this = s;  
		return *this;  
	}  

	IN Big_integer operator + (const Big_integer& b){  
		Big_integer c = *this; R int i;  
		for (i = 0; i < b.len; i++){  
			c.d[i] += b.d[i];  
			if (c.d[i] > 9) c.d[i]%=10, c.d[i+1]++;  
		}  
		while (c.d[i] > 9) c.d[i++]%=10, c.d[i]++;  
		c.len = max(len, b.len);  
		if (c.d[i] && c.len <= i) c.len = i+1;  
		return c;  
	}  

	IN Big_integer operator - (const Big_integer& b){  
		Big_integer c = *this; R int i;  
		for (i = 0; i < b.len; i++){  
			c.d[i] -= b.d[i];  
			if (c.d[i] < 0) c.d[i]+=10, c.d[i+1]--;  
		}  
		while (c.d[i] < 0) c.d[i++]+=10, c.d[i]--;  
		c.clean();  
		return c;  
	}  

	IN Big_integer operator * (const Big_integer& b)const{  
		R int i, j; Big_integer c; c.len = len + b.len;   
		for(j = 0; j < b.len; j++) for(i = 0; i < len; i++)   
			c.d[i+j] += d[i] * b.d[j];  
		for(i = 0; i < c.len-1; i++)  
			c.d[i+1] += c.d[i]/10, c.d[i] %= 10;  
		c.clean();  
		return c;  
	}  

	IN Big_integer operator / (const Big_integer& b){  
		R int i, j;  
		Big_integer c = *this, a = 0;  
		for (i = len - 1; i >= 0; i--)  
		{  
			a = a*10 + d[i];  
			for (j = 0; j < 10; j++) if (a < b*(j+1)) break;  
			c.d[i] = j;  
			a = a - b*j;  
		}  
		c.clean();  
		return c;  
	}  

	IN Big_integer operator % (const Big_integer& b){  
		R int i, j;  
		Big_integer a = 0;  
		for (i = len - 1; i >= 0; i--)  
		{  
			a = a*10 + d[i];  
			for (j = 0; j < 10; j++) if (a < b*(j+1)) break;  
			a = a - b*j;  
		}  
		return a;  
	}  

	IN Big_integer operator += (const Big_integer& b){  
		*this = *this + b;  
		return *this;  
	}  

	IN bool operator <(const Big_integer& b) const{  
		if(len != b.len) return len < b.len;  
		for(R int i = len-1; i >= 0; i--)  
			if(d[i] != b.d[i]) return d[i] < b.d[i];  
		return false;  
	}  
	IN bool operator >(const Big_integer& b) const{return b < *this;}  
	IN bool operator<=(const Big_integer& b) const{return !(b < *this);}  
	IN bool operator>=(const Big_integer& b) const{return !(*this < b);}  
	IN bool operator!=(const Big_integer& b) const{return b < *this || *this < b;}  
	IN bool operator==(const Big_integer& b) const{return !(b < *this) && !(b > *this);}  

	string str() const{  
		char s[maxn]={};  
		for(int i = 0; i < len; i++) s[len-1-i] = d[i]+'0';  
		return s;  
	}  
};  

istream& operator >> (istream& in, Big_integer& x)  
{  
	string s;  
	in >> s;  
	x = s.c_str();  
	return in;  
}  

ostream& operator << (ostream& out, const Big_integer& x)  
{  
	out << x.str();  
	return out;  
}  
}
using namespace BigInteger;
Big_integer fac1, fac2, Gcd;
Big_integer gcd(Big_integer x, Big_integer y) {return y == 0 ? x : gcd(y, x % y);}
Big_integer fpow(Big_integer base, R int tim)
{
	Big_integer ret = 1;
	W (tim)
	{
		if (tim & 1) ret = ret * base;
		base = base * base, tim >>= 1;
	}
	return ret;
}
int main(void)
{
	int T, n, k;
	ios_base::sync_with_stdio(false);
	cin >> T;
	W (T--)
	{
		cin >> n >> k;
		if (n > k) {cout << "0 1" << endl; continue;}
		fac1 = fpow(k + 1, n - 1);
		fac2 = fpow(k, n);
		Gcd = gcd(k - n + 1, fac2);
		fac1 = fac1 * ((Big_integer)(k - n + 1) / Gcd), fac2 = fac2 / Gcd;
		cout << fac1 << " " << fac2 << endl;
	}
}

LPA20020220

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
[Luogu P3330] [BZOJ 2227] [ZJOI2011]看电影

洛谷传送门BZOJ传送门题目描述到了难得的假期，小白班上组织大家去看电影。但由于假期里看电影的人太多，很难做到让全班看上同一场电影，最后大家在一个偏僻的小胡同里找到了一家电影院。但这家电影院分配座位的方式很特殊，具体方式如下：电影院的座位共有KKK个，并被标号为1...K1...K1...K，每个人买完票后会被随机指定一个座位，具体来说是从1...K1...K1...K中等可能的随机选取...
复制链接

扫一扫