2017年08月_d12155214552

原创概率期望 LightOJ-1038

https://vjudge.net/contest/182687#problem/C题意：给出一个数字Ｄ我们可以选择１－Ｄ中可以被Ｄ整除的数字，然后用Ｄ出得到一个新的数字Ｄ１；然后在找所有Ｄ１的因子，用Ｄ１除，直到得到１；问除的次数的期望值；思路： d[i] 代表从ｉ除到０的期望步数；那么假设ｉ一共有ｃ个因子（包括１和本身） d[i] = ( d[1] + d[a2] + d[

2017-08-31 20:40:18 193

原创矩阵快速幂+高斯消元解递推方程-HDU6185

http://acm.hdu.edu.cn/showproblem.php?pid=6185题目描述Bob’s school has a big playground, boys and girls always play games here after school. To protect boys and girls from getting hurt when playing happi

2017-08-31 19:16:07 2162 1

https://vjudge.net/contest/182427#problem/BDescriptionGiven 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 lin

2017-08-30 22:08:50 621

原创数学-HDU4790

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4790 题目意思是：从区间[a,b]选取一个数字x，区间[c,d]选取一个数字y，使得(x+y)%p=m。问你几率有多大。最简分式输出（除gcd就可以）。总的事件总数就是(b-a+1)*(d-c+1)。关键是求满足(x+y)%p=m 的选择方法总数。如何求呢？其实直接求不好求，而且容易超

2017-08-28 20:52:15 295

原创数学模拟 HDU-4814

题目地址：https://vjudge.net/problem/HDU-4814题意很简单：给定一个十进制数，要求输出在(1+sqrt(5))/2进制中的数的表示。 φ + 1 = φ 2 .11(φ) = 100(φ).由于x的0次幂等于1 所以完全可以把n当成对应的x进制数的个位现在要做的就是把这个数转化成二进制题目给出的公式可以在两边同时乘x的若干次幂那么这两个式子就变成了 x

2017-08-27 21:05:46 214

转载数学学习计划

1.burnside定理，polya计数法这个大家可以看brudildi的《组合数学》，那本书的这一章写的很详细也很容易理解。最好能完全看懂了，理解了再去做题，不要只记个公式。 *简单题：（直接用套公式就可以了） pku2409 Let it Bead http://acm.pku.edu.cn/JudgeOnline/problem?id=2409

2017-08-26 10:56:59 573

原创各种博弈-HDU5754

题意：一个国际象棋棋盘，有四种棋子，从(n,m)走到(1,1)，走到(1,1)的人赢，先手赢输出B，后手赢输出G，平局输出D。题解：先把从(n,m)走到(1,1)看做是从(1,1)走到(n,m)。四种棋子的规则如下：1、王（King）：横、竖、斜都可以走，每次限走一格2、车（Rook）：横、竖均可走，不能斜走，格数不受限制，除王车易位的情况下，平时不能越子3、马（Knight）：每步棋先横走或竖走一

2017-08-26 10:23:50 164

转载积性函数性质-POJ2480

这个题需要用到欧拉函数的知识……欧拉函数的定义：对正整数n，欧拉函数是小于n的正整数中与n互质的数的数目（我们定义φ(1)=1）欧拉函数的的通式：φ(n)=n*(1-1/p1)(1-1/p2)(1-1/p3)……*(1-1/ps)(p1,p2,p3,……ps均是n的质因子)欧拉函数有这么几个比较重要的性质：性质1：如果n是质数p的k次幂，那么φ(n)=p^k-1*(p-1)性质2：欧拉函数是积性函数

2017-08-25 21:23:36 436

转载二分-POJ1905

大致题意：一根两端固定在两面墙上的杆受热弯曲后变弯曲求前后两个状态的杆的中点位置的距离解题思路：几何和二分的混合体如图，蓝色为杆弯曲前，长度为L 红色为杆弯曲后，长度为s h是所求依题意知 S=(1+n*C)*L又从图中得到三条关系式; （1）角度→弧度公式 θr = 1/2*s （2）三角函数公式 sinθ= 1/2*L/r （3）

2017-08-25 20:59:22 249

原创 SG函数打表-HDU1848

题意：取石子问题，一共有3堆石子，每次只能取斐波那契数个石子，先取完石子者胜利，问先手胜还是后手胜可选步数为斐波那契数，用SG函数计算，最终结果是所有SG值异或的结果#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>#define ll long longus

2017-08-25 19:18:23 349

原创生成函数（背包）-HDU1171

https://vjudge.net/contest/181055#problem/C生成函数的解法要注意最后的判断，生成函数的用时更长。背包是标准解法。#include <iostream>#include <cstdio>#include <cstring>#include <cmath>#include <algorithm>#include <vector>#include

2017-08-23 14:28:26 501

转载生成函数

原文：http://blog.csdn.net/xiaofei_it/article/details/17042651母函数，又称生成函数，是ACM竞赛中经常使用的一种解题算法，常用来解决组合方面的题目。本文讲解母函数，但不讲解该算法的基础理论。读者随便找一本组合数学教材便可找到相应的内容，或者直接在网上搜索一下。母函数通常解决类似如下的问题：给5张1元，4张2元，3张5元，要得到15元，有多少种组

2017-08-23 13:11:58 166

原创生成函数-HDU1398

https://vjudge.net/contest/181055#problem/A给你一个数让你求出1~17中每个数的平方任意组合有几种情况满足条件。有多少种满足方案输出a[n] 最小的取不到的面值输出i，此时a[i]为0且i最小#include <iostream>#include <cstdio>#include <cstring>#include <cmath>#inclu

2017-08-23 12:34:10 161

原创容斥＋莫比乌斯反演＋分块优化-BZOJ2301

对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y) = k，gcd(x,y)函数为x和y的最大公约数。 Input 第一行一个整数n，接下来n行每行五个整数，分别表示a、b、c、d、k Output 共n行，每行一个整数表示满足要求的数对(x,y)的个数 Sample Input 22 5 1 5 11 5 1 5 2 Sample Out

2017-08-21 21:18:50 221

转载莫比乌斯反演-HDU1695

https://vj.xtuacm.cf/contest/view.action?cid=57#problem/OHDU1695的题目大意是这样的，给你 a , b , c , d , k 五个值（题目说明了你可以认为 a=c=1） x 属于 [1,b] ，y属于[1,d] 让你求有多少对这样的（x,y）满足gcd（x,y）==k。给你的时间是 3000 MS。 0 < a <= b

2017-08-21 19:54:00 252

原创莫比乌斯函数-BZOJ2440

其实这是我在某个莫比乌斯反演的PPT里看到的，但是这个题不是反演只是个莫比乌斯函数的应用。具体做法是二分答案。只需要一个小小的check函数来判断当前二分到的答案是否比k大或小即可。手动模拟了一下发现某个规律对于一个数t，t以内的数里的非完全平方数倍数的个数 num=1的倍数的数量−一个质数平方数(9,25,49…)的倍数的数量+两个质数的积平方数(36,100,225…

2017-08-21 13:48:16 285

转载阶梯博弈-HDU3389

其实这道题不难，不要理解错了就行，就像我一开是把题目中的A，B理解成盒子中卡片的数量了，就直接做不出来了，所以由此可见，认真读题是多么的重要！！！ 1 3 4号盒子是不能够再转移卡片到其他盒子中去了的，其他盒子中的卡片经过若干步的转移最终也一定会转移到1 3 4号盒子中去。具体来说，n % 6 == 0 或 2 或 5的盒子，经过奇数步转移到1 3 4中去，其他的则须经过偶数步才能转移过去。下面来

2017-08-18 10:54:08 242

转载阶梯博弈

今天在POJ做了一道博弈题..进而了解到了阶梯博弈…下面阐述一下我对于阶梯博弈的理解..首先是对阶梯博弈的阐述…博弈在一列阶梯上进行…每个阶梯上放着自然数个点..两个人进行阶梯博弈…每一步则是将一个集体上的若干个点( >=1 )移到前面去..最后没有点可以移动的人输.. 如这就是一个阶梯博弈的初始状态 2 1 3 2 4 … 只能把后面的点往前面放…如何来分析这个问题呢…其实阶梯博弈经过转换可以

2017-08-17 22:04:23 254

原创阶梯博弈-HDU4315

https://vj.xtuacm.cf/contest/view.action?cid=125#rank参考：http://blog.csdn.net/shahdza/article/details/7858032题意：有n个人爬山，山顶坐标为0，其他人按升序给出，不同的坐标只能容纳一个人（山顶不限），Alice和Bob轮流选择一个人让他移动任意步，但不能越过前面的人，且不能和前面一个人在相同的位

2017-08-17 22:00:45 286

原创组合数+二项式反演（容斥）-UVALive - 7040

题意：t组数据，每组给定n,m,k。有n个格子，m种颜色，要求把每个格子涂上颜色且正好适用k种颜色且相邻的格子颜色不同，求一共有多少种方案，结果对1e9+7取余。二项式反演（重点）设g(i)表示正好有i个颜色涂格子那么全部的C(n, i)*g(i)加起来正好就是所有涂格子的情况，总共n!种情况n! = Σ C(n, i)*g(i) (i从0到n)那么f(n) = n!，所以f(x) = x!那么我们

2017-08-17 13:10:27 379

原创博弈规律-HDU4642

题目: http://acm.hdu.edu.cn/showproblem.php?pid=4642题意: 给你一个n*m的矩阵,上面有正面或者反面的棋子, 每一次只能选(x,y)正面的棋子进行翻转, 区域(x,y)-(n,m) 中的棋子都会被翻转. A,B一块玩, A先,最后要使得所有的棋子都是反面的. 求出谁会赢得比赛.思路: 博弈. 因为每一次翻转都会翻到格子(n,m), 最后要使得所

2017-08-16 21:36:22 196

原创博弈推理-UVA11892

做博弈题自己先尝试推一推，先找一些简单的例子试一试题目大意：给定n堆石子的个数，两人轮流选择石子堆取石子，直到不能取为失败，附加条件，如果前一次操作，即队手的操作，没有将选中石子堆中的石子取完，那么当前操作者必须在该堆中取石子。解题思路：只要有一个石子堆的个数大于2，那么先手就获得必胜态，可控。对于全是1的情况判断奇偶性。简单推一下前几组，可以发现只有(1,1),(1,1,1,1),(1,1,1,1

2017-08-16 20:05:37 292

原创快速幂或找规律求a的b次方的最后一位-HDU1097

https://vj.xtuacm.cf/contest/view.action?cid=125#problem/EA hard puzzle Problem Description lcy gives a hard puzzle to feng5166,lwg,JGShining and Ignatius: gave a and b,how to know the a^b.everybody

2017-08-16 15:43:14 1239

原创威佐夫博弈输出第一次如何取-HDU2177

https://vj.xtuacm.cf/contest/view.action?cid=125#problem/D暴力枚举第一个分支，从两堆物品中同时取出相同数量的物品第二个分支，只从一堆物品中取物品（从多的那一堆中取，为什么只从多的那一堆中取捏?Because 从少的那一堆中取与从多的那一堆中取中会有一部分重合的地方，重合的就是从少的那一堆取的情况。）#include <iostream>

2017-08-16 15:14:03 267

原创博弈-HDU2897

邂逅明下 https://vj.xtuacm.cf/contest/view.action?cid=125#problem/BProblem Description 当日遇到月，于是有了明。当我遇到了你，便成了侣。那天，日月相会，我见到了你。而且，大地失去了光辉，你我是否成侣？这注定是个凄美的故事。（以上是废话）小t和所有世俗的人们一样，期待那百年难遇的日食。驻足街头看天，看日月渐渐走近

2017-08-16 14:23:03 183

原创高斯消元求解模线性方程-poj2947

https://vj.xtuacm.cf/contest/view.action?cid=115#problem/P这题目难点在于矩阵构建和取模strcmp(s, “MON”)==0 如果等于0两个字符串相同#include<stdio.h>#include<algorithm>#include<iostream>#include<string.h>#include<math.h>us

2017-08-10 15:11:00 267

原创高斯消元模板（kuangbin版）

//基本版 #include<stdio.h> #include<algorithm> #include<iostream> #include<string.h> #include<math.h> using namespace std; const int MAXN=50; int a[MAXN][MAXN];//增广矩阵 int x[MAXN];//解集 bo

2017-08-10 14:43:14 211

原创矩阵构造+矩阵快速幂-HDU5950

https://vj.xtuacm.cf/contest/view.action?cid=115#problem/I题目的难点在于矩阵的构造题解：典型的矩阵快速幂的运用。关键是i^4怎么维护？我们可以当成求第i+1项，那么i^4就变成了（i+1）^4。那么这时我们可以用二项式定理从i^4、i^3、i^2、i^1、i^0的组合中得到（i+1）^4。也就是说总共需要维护：f[i+1]、f[i]、（i+1

2017-08-09 10:30:18 241

原创高斯消元求方案数-HDU3364

https://vjudge.net/problem/HDU-3364#include<iostream>#include<stdlib.h>#include<stdio.h>#include<string>#include<vector>#include<deque>#include<queue>#include<algorithm>#include<set>#include<m

2017-08-08 15:56:31 200

原创高斯消元求解-HDU3359

https://vj.xtuacm.cf/contest/view.action?cid=115#problem/M这个题目用的模板是直接求解方程式题解及代码：这是一个double精度的高斯消元，方程建立直接暴力枚举每个点与当前点的距离就可以了，注意一下题目中给出的输入先是列数，后是行数#include <iostream>#include <cstdio>#include <cmath>

2017-08-08 15:35:57 203

原创高斯消元求方案数 POJ-1830

https://vj.xtuacm.cf/contest/view.action?cid=115#problem/N 题目难点在于建立方程组，求解直接用模板高斯消元和HDU的3364一样是开关问题 http://blog.csdn.net/luricheng/article/details/52506388 这里只是多了个小变形如果初始状态是010 目标状态是111 那其实就是等价于为

2017-08-08 15:32:49 249

原创 HDU 2814 斐波那契循环节+欧拉函数降幂

https://vj.xtuacm.cf/contest/view.action?cid=115#status//-/0/ 求F(a^b)^(F(a^b)^(n-1))%c a,b,n是longlong范围内的正整数 c是一个小于300的余数思路：不难发现c小于300是有意而为的这里有两个定理： 1.在模c的状态下，斐波那契具有循环节性质，即F[a^b]=F[a^b%len] len

2017-08-07 10:46:55 360

原创二分-UVA10341

https://vj.xtuacm.cf/contest/view.action?cid=115#problem/J 二分简单题#include <iostream>#include <cstdio>#include <cmath>#include <cstring>#include <algorithm>#define ll long longusing namespace std;

2017-08-05 10:05:17 179

原创三分-HDU3714

https://vj.xtuacm.cf/contest/view.action?cid=115#problem/D三分简单题#include <iostream>#include <cstdio>#include <cmath>#include <cstring>#include <algorithm>#define ll long longconst int maxn=10000;

2017-08-05 10:02:27 194

原创唯一分解定理+容斥-HDU4497

https://vj.xtuacm.cf/contest/view.action?cid=59#problem/L 这个题目的关键在于map的用法，通过map的灵活运用减少代码量和思维难度题意：求同时满足gcd(x,y,z)==g&&lcm(x,y,z)==L的(x,y,z)组合数.题解： 1.很直观的：若l%g！=0||g>l时，答案为0； 2.若果满足l%g==0,用唯一分解定理把个g

2017-08-03 19:11:14 234

原创 1-N中素数的个数（N非常大）-HDU5901

https://vj.xtuacm.cf/contest/view.action?cid=59#problem/D 题意：求区间[1,N]的质数的个数(1≤N≤10^11)。思路：模板题，上两个模板模板代码一：复杂度大概O（n^(3/4)） #include <bits/stdc++.h> #define ll long long using namespace

2017-08-03 11:48:41 2544

原创 lucas定理+费马小定理方法求逆元-HDU3944

https://vj.xtuacm.cf/contest/view.action?cid=59#problem/E 这题目用的线性方法求逆元要记住题目的意思是说，在杨辉三角上找一条路径，每次必须往下走一行，使得这个路径上所有数字的总和最小，输出取模值。这里对于组合数有一个变换。C(n-1,k-1)+C(n-1,k)=C(n,k)同时，仔细观察杨辉三角就知道，对于每行，从左至中间都是递增的，同时又根

2017-08-03 10:30:40 309

原创二次筛法找区间两点特别大（超过10亿）的素数-POJ2689

https://vj.xtuacm.cf/contest/view.action?cid=59#problem/C 如果区间两点小于3亿可以用bitset 这题区间两点可能大于20亿，只能用两次筛法，先找出0到47000之间的素数，再对目标区间用筛法。两次筛法是模板参考： http://blog.csdn.net/userluoxuan/article/details/38404185

2017-08-02 20:48:46 450

组合数学第五版

空空如也