题目链接:Portal Vjudge

Solution

一开始拿到这一题Sb了,把空放到dp中一起考虑了,这样计数就变得很麻烦。

其实我们可以把空位拿出来,假设它是存在的,最后再放回去。

那么就可以钦定某个人放到第一个,因为这是环排列,所以最后乘以n就可以了。

设\(Dp[i][j]\)表示来了\(i\)个人,分成\(j\)个联通块的方案数。

那么新来的人有三种选择,加入一个联通块,独立成为一个联通块,合并两个联通块:

\[Dp[i + 1][j + 1] += dp[i][j] * j\\
Dp[i + 1][j - 1] += dp[i][j] * j\\
Dp[i + 1][j] += dp[i][j] * 2j
\]

最后利用隔板法:$$ ans = n\sum_{i = 1}^{G} dp[m][i] * {n - k - 1 \choose i - 1}$$

一般环排列有空的计数,可以把空位全部抽出来计数,并且在计数过程中假定空存在,最后再用隔板法插入。

如果所求排列是环排列,那就钦定一种选法然后乘上环的大小即可。

Code

#include<bits/stdc++.h>
using namespace std;
#define rep(i, a, b) for(int i = (a), i##_end_ = (b); i <= i##_end_; ++i)
#define drep(i, a, b) for(int i = (a), i##_end_ = (b); i >= i##_end_; --i)
#define clar(a, b) memset((a), (b), sizeof(a))
#define debug(...) fprintf(stderr, __VA_ARGS__)
typedef long long LL;
typedef long double LD; const int Maxn = 4009, Mod = 1e9 + 7;
class Seatfriends {
private:
int dp[Maxn][Maxn];
int invFac[Maxn], fac[Maxn];
public:
int fpm(int base, int tims) {
int r = 1;
while (tims) {
if (tims & 1) r = r * 1ll * base % Mod;
base = base * 1ll * base % Mod;
tims >>= 1;
}
return r;
}
int C(int n, int m) {
if (m > n) return 0;
return fac[n] *1ll * invFac[m] % Mod * invFac[n - m] % Mod;
} void init(int n) {
fac[0] = 1;
rep (i, 1, n) fac[i] = 1ll * fac[i - 1] * i % Mod;
invFac[n] = fpm(fac[n], Mod - 2);
drep (i, n - 1, 0) invFac[i] = invFac[i + 1] * (i + 1ll) % Mod;
} int countseatnumb(int n, int k, int g) {
init(n + k); clar(dp, 0);
dp[1][1] = 1;
rep (i, 1, k)
rep (j, 1, g)
if (dp[i][j]) {
(dp[i + 1][j + 1] += 1ll * dp[i][j] * j % Mod) %= Mod;
(dp[i + 1][j] += 2ll * j % Mod * dp[i][j] % Mod) %= Mod;
if (j > 1) (dp[i + 1][j - 1] += 1ll * dp[i][j] * j % Mod) %= Mod;
}
LL ans = 0;
if (n == k) return (g >= 1) * dp[k - 1][1] * 1ll * n % Mod;
rep (i, 1, g) ans = (ans + 1ll * dp[k][i] * C(n - k - 1, i - 1) % Mod) % Mod;
return ans * 1ll * n % Mod;
}
};

[SRM625 Div1 Hard] Seatfriends的更多相关文章

  1. CF#345 (Div1)

    论蒟蒻如何被cf虐 以下是身败名裂后的题解菌=========== Div1 A.Watchmen 有n个点,每个点有一个坐标.求曼哈顿距离=欧几里得距离的点对数量. 只需要统计x或y一样的点对数量. ...

  2. 图论 SRM 674 Div1 VampireTree 250

    Problem Statement      You are a genealogist specializing in family trees of vampires. Vampire famil ...

  3. jq对象转为dom对象:$(".div1")[0] dom对象转为jq对象:$(dom对象)

    <!DOCTYPE html><html> <head> <meta charset="UTF-8"> <title>& ...

  4. 第一次div1做出3道题

    第一次div1做出3道题! 再接再厉! 哈利路亚!

  5. 第一次进div1了

    第一次进div1~好激动啊! 上帝依旧那么眷顾我!

  6. TopCoder 649 div1 & div2

    最近一场TC,做得是在是烂,不过最后challenge阶段用一个随机数据cha了一个明显错误的代码,最后免于暴跌rating,还涨了一点.TC题目质量还是很高的,非常锻炼思维,拓展做题的视野,老老实实 ...

  7. SRM DIV1 500pt DP

    SRM 501 DIV1 500pt SRM 502 DIV1 500pt SRM 508 DIV1 500pt SRM 509 DIV1 500pt SRM 511 DIV1 500pt SRM 5 ...

  8. codeforces #305 div1 done

    总算搞定了这一场比赛的题目,感觉收获蛮大 其中A,B,C都能通过自己的思考解决掉 D题思路好神,E题仔细想想也能想出来 以后坚持每两天或者一天做一场CF的div1的全套题目 除非有实在无法做出来的题目 ...

  9. [Codeforces Round #254 div1] C.DZY Loves Colors 【线段树】

    题目链接:CF Round #254 div1 C 题目分析 这道题目是要实现区间赋值的操作,同时还要根据区间中原先的值修改区间上的属性权值. 如果直接使用普通的线段树区间赋值的方法,当一个节点表示的 ...

随机推荐

  1. CEF3研究(三)

    一.Off-Screen Rendering 脱屏绘制 CEF的脱屏渲染并不创建源生的浏览器窗口,而是CEF提供主应用程序在无效区域和像素buffer里渲染,然后主应用程序通过鼠标.键盘和焦点事件通知 ...

  2. python爬虫(二)--了解deque

    队列-deque 有了上面一节的基础.当然你须要全然掌握上一节的全部方法,由于上一节的方法.在以下的教程中 会重复的用到. 假设你没有记住,请你返回上一节. http://blog.csdn.net/ ...

  3. Selenium 工具介绍、安装、配置及使用

    概念 selenium是一款跨平台开源工具,用于web的功能完整性測试,想做压力測试的能够换软件了 先看一眼官方站 http://docs.seleniumhq.org/download/ 首先要FQ ...

  4. 借助ltp语义分析提取特征,之后,文本生成

    """地点-哪里有做-业务-的(正规|靠谱)-公司?地点-做-业务-的(正规|靠谱)-公司(有哪些?|的联系方式是什么?|哪家口碑好值得信赖?)地点-做-业务-(怎么能省 ...

  5. sublime text3使用插件SublimeTextTrans设置透明度

    推荐一款在windows上设置sublime2和sublime3背景透明度的插件:SublimeTextTrans, 下载地址:https://github.com/vhanla/SublimeTex ...

  6. mongodb09----replicattion set--健壮性

    replication set复制集 replicattion set 多台服务器维护相同的数据副本,提高服务器的可用性.一台是服务器出问题了另外2台还可以接收干,secondary平时保持只读状态, ...

  7. 【转】 Android Studio --“Cannot resolve symbol” 解决办法

    Android Studio 无法识别同一个 package 里的其他类,将其显示为红色,但是 compile 没有问题.鼠标放上去后显示 “Cannot resolve symbol XXX”,重启 ...

  8. POJ3692 Kindergarten —— 二分图最大团

    题目链接:http://poj.org/problem?id=3692 Kindergarten Time Limit: 2000MS   Memory Limit: 65536K Total Sub ...

  9. 织梦仿站列表页pagelist分页显示竖排,如何修改成横排?

    织梦仿站列表页pagelist分页显示竖排,如何修改成横排? 织梦列表页的分页标签是采用pagelist来进行调用的,但是很多人在调用之后会出现一个列表竖着排列的问题(横排美观度好一些),还是非常不美 ...

  10. mac多线程下载神器

    本文参考:https://blog.csdn.net/orangleliu/article/details/46834429 神器:axel 安装(已经安装homebrew前提下,没有请参考:http ...