数学相关
Saikaai
这个作者很懒,什么都没留下…
展开
-
BZOJ 1076 [SCOI 2008] 奖励关 (概率与期望)
题目链接:BZOJ 1076#include#include#includeusing namespace std;int K,N;int a[50],v[50];double dp[110][32800];inline int read(){ int x=0,f=1;char ch=getchar(); while(ch'9'){if(ch=='-')f=-1;ch=g原创 2015-03-12 12:42:37 · 595 阅读 · 0 评论 -
bzoj 3112(zjoi 2013 防守战线) 线性规划+网络流
题目链接http://www.lydsy.com/JudgeOnline/problem.php?id=3112话说这道题真心变态。。。终于见识到了浙江省选的难度了。。。我们根据这道题的样例容易写出这样的不等式组:x2+x3>=1x1+x2+x3+x4+x5>=4x3+x4+x5>=2所求为min{x1+5*x2+6*x3+3*x4+4*x5},其中x[i]表示在i位置修建原创 2015-02-10 11:59:59 · 2695 阅读 · 0 评论 -
bzoj 1257(数学)
Description给出正整数n和k,计算j(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其中k mod i表示k除以i的余数。例如j(5, 3)=3 mod 1 + 3 mod 2 + 3 mod 3 + 3 mod 4 + 3 mod 5=0+1+0+3+3=7Input输入仅一行,包含两个整数n, k。原创 2015-02-05 22:04:26 · 509 阅读 · 0 评论 -
poj1006(中国剩余定理)
DescriptionSome people believe that there are three cycles in a person's life that start the day he or she is born. These three cycles are the physical, emotional, and intellectual cycles, and the原创 2015-02-05 21:58:15 · 464 阅读 · 0 评论 -
codevs 1482 路线统计
【问题描述】给出一个n个点有向图,求从s点到f点恰好经过时间t的路径总数。不能,可以重复的走各点。【输入数据】第一行包含一个整数n, 所有点是从0到n-1编号.接下来n行,每行包含n个字符.第i行第j个字符表示i到j需要的时间,字符只可能是’1’到’5’,或者是’.’表示i不能到达j, 保证主对角线都是’.’。接下来一行3个整数s, f, t。【输出数据】输原创 2014-10-30 16:50:25 · 555 阅读 · 0 评论 -
概率和期望
这几天集训有几道概率和期望的题,才知道到我之前的概率和期望的题简直就像小学生数数一样QAQ。所以我打算看一下概率和期望。下面是一些总结,至于那几道非人哉的题之后弄清楚了再写吧。不定期更新。1.连续型变量的期望:概率密度积分什么的。2.离散型变量的期望:E(x)=sigma(x[i]*p[i])3.平方的期望:E(x^2)=E(x)^2+D(x)证明:E(x^2)原创 2015-03-25 19:12:47 · 980 阅读 · 0 评论 -
Poj 2096 Collecting Bugs (概率与期望)
题目链接:poj 2096dp[i][j]:descrip in this step we have found i kinds of bugs and these bugs belong to j kinds of system.dp[i][j]=(i/j)*(n*s)*dp[i][j]+(i/j)*(1-j/s)*dp[i][j+1]+(1-i/n)*(j/s)*dp[i+1][j]+原创 2015-03-19 22:04:49 · 571 阅读 · 0 评论 -
Poj 3734 Blocks(DP,矩阵乘法优化)
题目链接:poj 3734这道题用矩阵乘法优化DP。考虑到直接转移的话,N太大,会TLE。由于转移的方案数,转移的状态很少,所以可以将转移的方案用矩阵来表示。设dp[i][k]表示涂到第i个格子,状态为k的方案数。其中状态k的定义为:k=1:红色和绿色都为奇数,k=2:红奇绿偶,k=3:红偶绿奇,k=4:红偶绿偶DP方程自己写。表示为矩阵为:2 1 1 01 2 0原创 2015-03-20 21:24:15 · 598 阅读 · 0 评论 -
BZOJ 1875 [SDOI 2009] HH去散步 (DP,矩阵乘法优化)
题目链接:BZOJ 1875这道题的思路,主要是构建矩阵的思路很巧妙。我们普通的用矩阵乘法转移是用点来转移,但是这样不能去掉在一个地方逗留的情况。一个很神奇的做法就是用边来构图(对于一条边i,除去一条边j满足i==(j^1)的情况,与其他的边都相连)转移t-1次,然后同用与起点相连的边构造的一个矩阵(相当于系数矩阵)相乘。最后统计答案,只需将终点相连的边的答案加上即可。#inc原创 2015-03-20 21:26:18 · 1449 阅读 · 1 评论 -
BZOJ 3560 DZY Loves Math V
题目链接:BZOJ 3560首先,可以根据phi函数为积性函数,可以分解质因数,计算每个质因数对答案的贡献,最后乘起来即可。那么现在问题就来了,怎样计算质因数对答案的贡献?我们可以对每个数ai分解质因数,对于它的一个质因数p,记录p出现的次数bi,那么这个质数p对答案的贡献,是对于每个数ai,p的j(1华丽丽的传送门。(不会用公式编辑器的哭了QAQ。)#include#include#原创 2015-04-01 20:45:49 · 801 阅读 · 0 评论 -
BZOJ 2440 完全平方数(莫比乌斯函数)
题目链接:BZOJ 2440这道题需要惊人的联想能力和很大的脑洞= =。。。然后将它联系到莫比乌斯函数上。我也是看神犇的博客的。首先容易想到二分答案,然后判断当前的满足条件的数。具体怎样得出的请大家自己YY吧。。。#include#include#include#includeusing namespace std;//利用莫比乌斯函数求值 #define maxn (5000原创 2015-02-27 10:15:17 · 552 阅读 · 0 评论 -
BZOJ 2301 Problem b (莫比乌斯反演)
题目链接:BZOJ 2301这道题求a为z的(x,y)的对数。然后我们就用莫比乌斯反演求出,具体推导什么的用公式编辑器太麻烦了。让我们一起膜拜PoPoQQQ的题解吧。。。华丽丽的传送门:神犇的博客#include#include#includeusing namespace std;#define maxn (50000+10)int vis[maxn],mu[maxn],原创 2015-02-27 10:20:51 · 483 阅读 · 0 评论 -
hdu 4405 Aeroplane chess (概率与期望)
题目链接:hdu 4405dp[i]表示在第i个格子,距离终点状态所需要的期望步数方程转移:dp[i]=dp[fa[i]] ->若当前位置可以飞行dp[i]+=(dp[i+j]+1)/6 (1#include#include#includeusing namespace std;#define maxn (1000000+10)int N,M;int fa[m原创 2015-03-11 13:57:56 · 502 阅读 · 0 评论 -
poj 3233 Matrix Power Series(矩阵乘法)
DescriptionGiven a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.InputThe input contains exactly one test case. The first line of input contains three原创 2014-11-01 07:32:07 · 449 阅读 · 0 评论 -
Rational Resistance
迈克希望组装一些电子器件,他有无限多的1Ω的单位电阻,现在他想要用最少的单位电阻组装成一个电阻为a/b的电子器件,单位电阻之间可以串联也可以并联[Intput]输入一行两个正整数a,b(1 ≤ a, b ≤ 1018),表示要组装的值a/b,为最简分数,数据保证有解。 [Output]输出一行一个答案。 [Sample ]Input原创 2014-10-23 22:12:24 · 730 阅读 · 0 评论 -
BZOJ 3527 [ZJOI 3014] 力 (FFT)
题目链接:BZOJ 3527因为这道题在BZOJ上没有贴,所以我还是附上数据吧。input:54006373.88518415375036.4357591717456.4691448514941.0049121410681.345880output:-16838672.6933439.7937509018.5664595686.88610原创 2015-03-01 19:35:17 · 601 阅读 · 0 评论 -
bzoj 3450(数学)
#include#include#includeusing namespace std;#define maxn (300000+10)int n;char a[maxn];double f[maxn],d[maxn];int main(){ scanf("%d\n",&n); scanf("%s",a); for(int i=1;i<=n;i++){ if(a[i原创 2015-02-05 22:10:08 · 617 阅读 · 0 评论 -
关于·拓展欧几里得
扩展欧几里德算法基本算法:对于不完全为 0 的非负整数 a,b,gcd(a,b)表示 a,b 的最大公约数,必然存在整数对 x,y ,使得 gcd(a,b)=ax+by。证明:设 a>b。 1,显然当 b=0,gcd(a,b)=a。此时 x=1,y=0; 2,ab!=0 时 设 ax1+by1=gcd(a,b); bx2+(a mod b)y2=gcd(b,a mo原创 2014-11-05 19:10:03 · 416 阅读 · 0 评论 -
codeforces Tower of Hanoi
[Description]汉诺塔游戏,在我们熟知的汉诺塔游戏的规则下增加一条,从第i柱移动一个盘子到j柱需要消耗tij的时间,你的任务是求n个盘子从1柱移动到3柱的最少时间。[Intput]首先是一个3*3的矩阵表示tij (1 ≤ tij ≤ 10000;i ≠ j),tii==0,接下来一行一个整数n(1 ≤ n ≤ 40),表示初始时1柱上的圆盘数。[Output]原创 2014-10-24 22:10:58 · 452 阅读 · 0 评论 -
POJ 3844 Divisible Subsequences
[Description]给定一个由正整数构成的序列,并给出一个数 d,求序列中存在多少个连续子序列,使得子序列之和能被 d 整除。[Input]第一行一个数 T,表示有 T 组测试数据。每组测试数据第一行两个数,d,n 表示除数和序列的长度接下来一行 n 个不大于 10^9 的数,表示这个序列。[Output]对于每组测试数据输出一个数,表示找到符合条件的子串的数原创 2014-10-18 07:43:05 · 676 阅读 · 0 评论 -
hdu 4336 Card Collector (概率与期望)
题目链接:hdu 4336dp[i]表示在i这个状态下,收集齐所有卡片所需要的买的袋数,转移方程为:dp[i]=(dp[i]+1)*(1-tot)+(dp[i|a[j]]+1)*p[j] (j为i状态中不包含的袋数,tot为所有符合条件的j的p[j]的和),即i可以转移到拓展一袋的状态,也可以转移到没有拓展的状态。由于方程左右两边都有要求的量dp[i]那么把等式右边展开,移项可以得到最后的原创 2015-03-11 19:10:22 · 427 阅读 · 0 评论 -
各种数论知识
博客好久没有更新了。。。感觉自己这段时间集训就像是在学数学一样ORZ。高一的那些触们会那么多高数是要闹哪样QAQ。下面是一些自己这段时间的学到的一些数论知识。我太弱,不要鄙视内容。。。原根1.设M为正整数,a为整数,若a % M 的阶(群论的定义下为元素的个数)为phi(M),则a为模M的一个原根。2.设g为P(P为素数)的原根,则满足,g ^ (P -1) = 1 (mo原创 2015-04-01 19:21:48 · 857 阅读 · 0 评论