UVa Online Judge

  训练指南的题目。

  题意是,给出n个人,以及一些关系,要求对这n个人构成一个排列,其中父亲必须排在儿子的前面。问一共有多少种方式。

  做法是,对于每一个父节点,将它的儿子结点构成的子树看成无序状态,这样子对当前父节点整棵树计算一个排列数。如果把所有的这样的式子写出来,可以发现分子分母是可以相消的。假设点的总数是S,儿子的点的数目分别是A,B,C...,这样的话,对于这个结点,可以求得F(S)=F(A)+F(B)+F(C)+...。最后的结果是所有子树的F(S)*F(A)*F(B)*F(C)*...。最后消去以后就只剩下F(S)/(c(A)*c(B)*c(C)*...),其中c(X)是结点X的子树的结点个数。

代码如下:

 #include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <map>
#include <vector> using namespace std; const int N = ;
map<int, int> id;
int lsf[N], pn, prm[N >> ]; void getprm() {
id.clear();
lsf[] = lsf[] = ;
pn = ;
for (int i = ; i < N; i++) {
if (!lsf[i]) {
id[i] = pn;
prm[pn++] = i;
for (int j = i; j < N; j += i) lsf[j] = i;
}
}
} typedef long long LL;
const LL MOD = ;
int cnt[N], fcn[N];
bool vis[N], fa[N];
vector<int> rel[N]; LL multi(int a, int b) {
LL ret = , p = a;
while (b > ) {
if (b & ) ret *= p, ret %= MOD;
p *= p, p %= MOD;
b >>= ;
}
return ret;
} void dfs(int x) {
cnt[x] = ;
for (vector<int>::iterator vi = rel[x].begin(); vi != rel[x].end(); vi++) {
dfs(*vi);
cnt[x] += cnt[*vi];
}
} int main() {
getprm();
int n, m, T;
scanf("%d", &T);
while (T-- && ~scanf("%d%d", &n, &m)) {
memset(vis, , sizeof(vis));
memset(cnt, , sizeof(cnt));
memset(fcn, , sizeof(fcn));
memset(fa, , sizeof(fa));
int x, y;
for (int i = ; i <= n; i++) rel[i].clear();
while (m--) {
scanf("%d%d", &x, &y);
rel[y].push_back(x);
fa[x] = true;
}
for (int i = ; i <= n; i++) if (!fa[i]) rel[].push_back(i);
dfs();
for (int i = ; i <= n; i++) {
int t = i;
while (t > ) {
fcn[id[lsf[t]]]++;
t /= lsf[t];
}
t = cnt[i];
while (t > ) {
fcn[id[lsf[t]]]--;
t /= lsf[t];
}
//cout << "~~ " << i << ' ' << cnt[i] << endl;
//for (int i = 0; i < 10; i++) cout << prm[i] << ' ' << fcn[i] << endl;
}
//cout << pn << endl;
LL ans = ;
for (int i = ; i < pn; i++) {
ans *= multi(prm[i], fcn[i]);
ans %= MOD;
}
cout << ans << endl;
}
return ;
}

——written by Lyon  

uva 11174 Stand in a Line (排列组合)的更多相关文章

  1. uva 11174 Stand in a Line

    // uva 11174 Stand in a Line // // 题目大意: // // 村子有n个村民,有多少种方法,使村民排成一条线 // 使得没有人站在他父亲的前面. // // 解题思路: ...

  2. UVA 11174 Stand in a Line 树上计数

    UVA 11174 考虑每个人(t)的所有子女,在全排列中,t可以和他的任意子女交换位置构成新的排列,所以全排列n!/所有人的子女数连乘   即是答案 当然由于有MOD 要求逆. #include & ...

  3. UVA 11174 Stand in a Line (组合+除法的求模)

    题意:村子里有n个人,给出父亲和儿子的关系,有多少种方式可以把他们排成一列,使得没人会排在他父亲的前面 思路:设f[i]表示以i为根的子树有f[i]种排法,节点i的各个子树的根节点,即它的儿子为c1, ...

  4. UVA 11174 Stand in a Line 树dp+算

    主题链接:点击打开链接 题意:白书的P103. 加个虚根就能够了...然后就是一个多重集排列. import java.io.PrintWriter; import java.util.ArrayLi ...

  5. UVA 11174 Stand in a Line,UVA 1436 Counting heaps —— (组合数的好题)

    这两个题的模型是有n个人,有若干的关系表示谁是谁的父亲,让他们进行排队,且父亲必须排在儿子前面(不一定相邻).求排列数. 我们假设s[i]是i这个节点,他们一家子的总个数(或者换句话说,等于他的子孙数 ...

  6. 【递推】【推导】【乘法逆元】UVA - 11174 - Stand in a Line

    http://blog.csdn.net/u011915301/article/details/43883039 依旧是<训练指南>上的一道例题.书上讲的比较抽象,下面就把解法具体一下.因 ...

  7. 数学:UVAoj 11174 Stand in a Line

    Problem J Stand in a Line Input: Standard Input Output: Standard Output All the people in the bytela ...

  8. UVa 12712 && UVaLive 6653 Pattern Locker (排列组合)

    题意:给定 一个n * n 的宫格,就是图案解锁,然后问你在区间 [l, r] 内的所有的个数进行组合,有多少种. 析:本来以为是数位DP,后来仔细一想是排列组合,因为怎么组合都行,不用考虑实际要考虑 ...

  9. UVa Problem 10132 File Fragmentation (文件还原) 排列组合+暴力

    题目说每个相同文件(01串)都被撕裂成两部分,要求拼凑成原来的样子,如果有多种可能输出一种. 我标题写着排列组合,其实不是什么高深的数学题,只要把最长的那几个和最短的那几个凑一起,然后去用其他几个验证 ...

随机推荐

  1. java.lang.NoSuchMethodError: org.apache.commons.io.FileUtils.copyInputStreamToFile(Ljava/io/InputStream;Ljava/io/File;)V

    昨天用的好好的,今天就不行了 也懒得搞 具体原因就是引用的问题 找了个能用的pom 直接贴过去完事儿

  2. JavaScript远程调用

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  3. spring boot 异常处理(转)

    spring boot在异常的处理中,默认实现了一个EmbeddedServletContainerCustomizer并定义了一个错误页面到"/error"中,在ErrorMvc ...

  4. 扩展 Microsoft.Owin.Security

    微软在 OWIN 框架中对 OAuth 认证的支持非常好, 使用现有的 OWIN 中间件可以做到: 使用 Microsoft.Owin.Security.OAuth 搭建自己的 OAuth2 服务端, ...

  5. 调用Lua脚本print(xxx)报attempt to call a nil value (global 'print')错误

    在自己程序里调用Lua脚本print(xxx) 报出attempt to call a nil value (global 'print')错误 解决方法: luaopen_base(L); 或者 l ...

  6. onethink二级导航调用

    <ul class="nav-main">//添加tree参数 <think:nav name="nav" tree="true&q ...

  7. php_imagick是怎么实现复古效果的呢?

    php_imagick程序示例 1.创建一个缩略图并显示出来 <?phpheader('Content-type: image/jpeg');$image = new Imagick('imag ...

  8. 日期格式之——new Date()的用法

    获取时间: 1 var myDate = new Date();//获取系统当前时间 获取特定格式的时间: 1 myDate.getYear(); //获取当前年份(2位) 2 myDate.getF ...

  9. 利用JDBC连接Oracle数据库(转)

    http://blog.csdn.net/wahaha1_/article/details/8512438 JDBC是Sun公司制定的一个可以用Java语言连接数据库的技术. 一.JDBC基础知识 J ...

  10. python 模块的执行环境