数学相关
Saikaai
这个作者很懒,什么都没留下…
展开
专栏收录文章
- 默认排序
- 最新发布
- 最早发布
- 最多阅读
- 最少阅读
-
BZOJ 1076 [SCOI 2008] 奖励关 (概率与期望)
题目链接:BZOJ 1076 #include #include #include using 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 · 656 阅读 · 0 评论 -
bzoj 3112(zjoi 2013 防守战线) 线性规划+网络流
题目链接http://www.lydsy.com/JudgeOnline/problem.php?id=3112 话说这道题真心变态。。。终于见识到了浙江省选的难度了。。。 我们根据这道题的样例容易写出这样的不等式组: x2+x3>=1 x1+x2+x3+x4+x5>=4 x3+x4+x5>=2 所求为min{x1+5*x2+6*x3+3*x4+4*x5},其中x[i]表示在i位置修建原创 2015-02-10 11:59:59 · 2827 阅读 · 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=7 Input 输入仅一行,包含两个整数n, k。原创 2015-02-05 22:04:26 · 569 阅读 · 0 评论 -
poj1006(中国剩余定理)
Description Some 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 · 519 阅读 · 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 · 621 阅读 · 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 · 1171 阅读 · 0 评论 -
Poj 2096 Collecting Bugs (概率与期望)
题目链接:poj 2096 dp[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 · 625 阅读 · 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 0 1 2 0原创 2015-03-20 21:24:15 · 661 阅读 · 0 评论 -
BZOJ 1875 [SDOI 2009] HH去散步 (DP,矩阵乘法优化)
题目链接:BZOJ 1875 这道题的思路,主要是构建矩阵的思路很巧妙。我们普通的用矩阵乘法转移是用点来转移,但是这样不能去掉在一个地方逗留的情况。 一个很神奇的做法就是用边来构图(对于一条边i,除去一条边j满足i==(j^1)的情况,与其他的边都相连)转移t-1次, 然后同用与起点相连的边构造的一个矩阵(相当于系数矩阵)相乘。最后统计答案,只需将终点相连的边的答案加上即可。 #inc原创 2015-03-20 21:26:18 · 1532 阅读 · 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 · 860 阅读 · 0 评论 -
BZOJ 2440 完全平方数(莫比乌斯函数)
题目链接:BZOJ 2440 这道题需要惊人的联想能力和很大的脑洞= =。。。然后将它联系到莫比乌斯函数上。我也是看神犇的博客的。首先容易想到二分答案,然后判断当前的满足条件的数。具体怎样得出的请大家自己YY吧。。。 #include #include #include #include using namespace std; //利用莫比乌斯函数求值 #define maxn (5000原创 2015-02-27 10:15:17 · 621 阅读 · 0 评论 -
BZOJ 2301 Problem b (莫比乌斯反演)
题目链接:BZOJ 2301 这道题求a 为z的(x,y)的对数。然后我们就用莫比乌斯反演求出,具体推导什么的用公式编辑器太麻烦了。让我们一起膜拜PoPoQQQ的题解吧。。。华丽丽的传送门:神犇的博客 #include #include #include using namespace std; #define maxn (50000+10) int vis[maxn],mu[maxn],原创 2015-02-27 10:20:51 · 541 阅读 · 0 评论 -
hdu 4405 Aeroplane chess (概率与期望)
题目链接:hdu 4405 dp[i]表示在第i个格子,距离终点状态所需要的期望步数 方程转移: dp[i]=dp[fa[i]] ->若当前位置可以飞行 dp[i]+=(dp[i+j]+1)/6 (1 #include #include #include using namespace std; #define maxn (1000000+10) int N,M; int fa[m原创 2015-03-11 13:57:56 · 554 阅读 · 0 评论 -
poj 3233 Matrix Power Series(矩阵乘法)
Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak. Input The input contains exactly one test case. The first line of input contains three原创 2014-11-01 07:32:07 · 494 阅读 · 0 评论 -
Rational Resistance
迈克希望组装一些电子器件,他有无限多的1Ω的单位电阻,现在他想要用最少的单位电阻组装成一个电阻为a/b的电子器件,单位电阻之间可以串联也可以并联 [Intput] 输入一行两个正整数a,b(1 ≤ a, b ≤ 1018),表示要组装的值a/b,为最简分数,数据保证有解。 [Output] 输出一行一个答案。 [Sample ] Input原创 2014-10-23 22:12:24 · 781 阅读 · 0 评论 -
BZOJ 3527 [ZJOI 3014] 力 (FFT)
题目链接:BZOJ 3527 因为这道题在BZOJ上没有贴,所以我还是附上数据吧。 input: 5 4006373.885184 15375036.435759 1717456.469144 8514941.004912 1410681.345880 output: -16838672.693 3439.793 7509018.566 4595686.886 10原创 2015-03-01 19:35:17 · 663 阅读 · 0 评论 -
bzoj 3450(数学)
#include #include #include using 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 · 679 阅读 · 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 · 463 阅读 · 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 · 499 阅读 · 0 评论 -
POJ 3844 Divisible Subsequences
[Description] 给定一个由正整数构成的序列,并给出一个数 d,求序列中存在多少个连续子序列,使得 子序列之和能被 d 整除。 [Input] 第一行一个数 T,表示有 T 组测试数据。 每组测试数据第一行两个数,d,n 表示除数和序列的长度 接下来一行 n 个不大于 10^9 的数,表示这个序列。 [Output] 对于每组测试数据输出一个数,表示找到符合条件的子串的数原创 2014-10-18 07:43:05 · 750 阅读 · 0 评论 -
hdu 4336 Card Collector (概率与期望)
题目链接:hdu 4336 dp[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 · 486 阅读 · 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 · 975 阅读 · 0 评论
分享