淡墨æ末央

摘要： 题目：http://poj.org/problem?id=2096题意：一个软件有 s 个子系统，存在 n 种 bug。某人一天能找到一个 bug。 问，在这个软件中找齐 n 种 bug，并且每个子系统中至少包含一个 bug 的时间的期望值（单位：天）。 注意：bug 是无限多的，每个 bug 属于任何一种 bug 的概率都是 1/n；出现在每个系统是等可能的，为 1/s。思路：令f[i][j]表示已经找到了 i 种 bug，且 j 个子系统至少包含一个 bug，距离完成目标需要的时间的期望。 目标状态是 f[0][0]～ 转移方程： f[i][j] = i/n*j/s*f[i][... 阅读全文

摘要： 题目：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3551题意：开始有N-1个人和一个吸血鬼， 每天有两个生物见面，当人遇到吸血鬼时有p的概率变成吸血鬼，求全部变成吸血鬼所需要的时间的期望~思路： 设dp[i] 为还有 i 个人时，有一人变成吸血鬼的期望时间， p[i]为还有 i 个人时，有人变成吸血鬼的概率， 那么p[i]= p*i(N-i)/(N*(N-1)/2)~ dp[i]=1/p[i]; 由 E(X)=∑E(X=xi) 得 E[N]=∑dp[i]~ 1 #include 2 #include ... 阅读全文

摘要： 题目：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3329题意：有三个均匀的骰子，分别有k1,k2,k3个面，初始分数是0，当掷三个骰子的点数分别为a,b,c的时候，分数清零， 否则分数加上三个骰子的点数和，当分数>n的时候结束。求需要掷骰子的次数的期望~思路： 设 dp[i] 为分数为i时还需要的期望值， 那么dp[i]= ∑(dp[i+k]*pk ) + dp[0]*p0 +1; 设dp[i]= A[i] *dp[0] + B[i]; 代入上式得 ： dp[i] = ∑( A[i+k]*dp[0] +... 阅读全文

摘要： 题目：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=2271题意：一个N*N的矩阵( N = 2*k-1 , k 2 #include 3 #include 4 using namespace std; 5 const int dir[ 4 ][ 2 ] = {1,0,0,1,-1,0,0,-1}; 6 double dp[110][110][110]; 7 int N; 8 double F( int x, int y ) 9 {10 if( x=N )return -1;11 if( (x==0... 阅读全文

摘要： 本文转自：http://www.cnblogs.com/AbandonZHANG/archive/2012/07/21/2601889.html组合数学计数与统计2001 - 符文杰：《Pólya原理及其应用》2003 - 许智磊：《浅谈补集转化思想在统计问题中的应用》2007 - 周冬：《生成树的计数及其应用》2008 - 陈瑜希《Pólya计数法的应用》数位问题2009 - 高逸涵《数位计数问题解法研究》2009 - 刘聪《浅谈数位类统计问题》动态统计2004 - 薛矛：《解决动态统计问题的两把利刃》2007 - 余江伟：《如何解决动态统计问题》博弈2002 - 张一飞 阅读全文

摘要： 题目：http://poj.org/problem?id=2888题意：给定n（n 2 #include 3 #include 4 #include 5 using namespace std; 6 const int Mod=9973; 7 int N, M, K, T; 8 struct Mar 9 { 10 int m[15][15]; 11 inline void zero( ){ 12 memset(m, 0, sizeof m); 13 } 14 inline void one( ){ 15 z... 阅读全文

摘要： 链接:http://uva.onlinejudge.org/external/107/10733.pdf题意: N 种颜色可以涂成多少种立方体~思路: 使正六面体保持不变的运动群总共有:1.不变置换(1)(2)(3)(4)(5)(6), 共1个;2.沿对面中心轴旋转 90度, 270度 (1)(2345)(6), (1)(5432)(6) 同类共 6个;3.沿对面中心轴旋转 180度 (1)(24)(35)(6), 同类共 3个;4.沿对角线轴旋转 120度, 240度 (152)(346), (251)(643) 同类共 8个;5.沿对边中点轴旋转 180度 (16)(25)(43) 同类共 阅读全文

摘要： 链接：http://poj.org/problem?id=2154题意：给出两个整数 N 和 P，表示 N 个珠子，N种颜色，要求不同的项链数， 结果 %p ~思路： 利用polya定理解~定理内容：设是n个对象的一个置换群, 用m种颜色染图这n个对象，则不同的染色方案数为：其中，为的循环节数~本题只有旋转一种置换方式，那么共有 N 个置换， 每个置换的循环节为 gcd(N,i)~那么结果为∑(N^(gcd(N, i))) %P。 N为 1e9, 不能枚举 i , 但我们可以统计 gcd(N,i)==a 的有多少个~令L==N/a, i==a*t, 即 a==gcd(N, i)==gcd(L* 阅读全文

摘要： 暴力法： 1 int i;2 for (i=1;;i++) { if (((long long int)(n)*i-an)%M==0) break; }欧拉函数1 long long inv( long long n )2 {3 return pow( n, M - 2 )%M;4 }扩展欧几里得 1 //扩展欧几里德 2 void exp_gcd( LL a ,LL b,LL &x,LL &y) { 3 if( b == 0 ) { 4 x = 1; 5 y = 0; 6 } 7 else { 8 ... 阅读全文

摘要： 链接：http://poj.org/problem?id=1286http://poj.org/problem?id=2409 1 #include 2 #include 3 #include 4 #include 5 #include 6 using namespace std; 7 typedef long long LL; 8 LL P_M( LL a, LL b ) 9 {10 LL res=1, t=a;11 while (b){12 if(b&1)res*=t;13 t*=t;14 b>>=1;15 ... 阅读全文