题目链接: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. 下一代的中间件必须是支持docker规范的

    下一代的中间件必须是支持docker规范的,这是中间件技术走向标准规范化的必经之路. 什么是 Docker? 答案是:Docker 是下一代的云计算模式.Docker 是下一代云计算的主流趋势. Do ...

  2. 定制 ArcEngine 要素编辑工具

    来自:http://blog.sina.com.cn/s/blog_4d780fc10101d2d5.html 先初步了解到大概用到的下面的接口和类: IEngineEditor IEngineEdi ...

  3. Python pandas学习笔记

    参考文献:<Python金融大数据分析> #导入模块 import pandas as pd #生成dataframe df = pd.DataFrame([10,20,30,40], c ...

  4. Visual Studio VS如何修改代码字体

    工具-选项-环境-字体和颜色

  5. 怎样使用SSH连接OpenStack上的云主机

    转载请注明出处.否则将追究法律责任http://blog.csdn.net/xingjiarong/article/details/47021815 在上一篇博客中我介绍了怎样在OpenStack中创 ...

  6. 获取当前时间 YYYY-MM-DD

    1.函数封装 /** * 获取当前时间 * 格式YYYY-MM-DD */ Vue.prototype.getNowFormatDate = function() { var date = new D ...

  7. Python的调用程序

    任务 调用系统命令ping 判断局域网内有哪些主机存活 假设你用c语言写了一个算法,需要对该算法进行测试.测试的数据集几百个.这时可以使用过GCC生成test.exe,再使用python批量调用该ex ...

  8. [办公应用]让WORD自动显示到四级目录

    一般情况下,word的目录默认显示到三级目录.如果需要显示到四级目录,你会怎么操作呢? 只要按下图所示,单击“插入”-“引用”-“索引和目录” 然后单击“目录”选项卡,将“显示级别”处的3改为4即可. ...

  9. mysql + Fluently NHibernate + WebAPI + Autofac

    MySQL.Fluently NHibernate.WebAPI.Autofac,对我来说每一个都是麻烦疙瘩,现在它们为了一个共同的项目而凑合到一起了.一路磕磕碰碰,现在貌似有了一点眉目. 作为一个步 ...

  10. spark安装和登陆配置

    1.下载安装包: http://www.igniterealtime.org/downloads/index.jsp 2.点击安装后打开登陆界面: 2.点击“高级”,设置相关配置: 3.点击“登陆”后 ...