DancingZ的博客

今天做了到题要用到慢速乘，然后发现自己的快速幂一直记得模板，没有弄懂原理= =快速幂：假设对于a^11，我们把11换成二进制11->1011那么a^11=a^8*a^2*a^1，转化成二进制之后，我们只用维护一个底数，就像十进制乘法那样。ll ksm(ll a,ll b,ll p){ ll ans=1,ret=a; for(;b;b>>=1){ if(b...

DescriptionBB很喜欢玩《文明》，因此它知道很多游戏的攻略。有一个攻略根据城市的N 种资源的数量，对地图上许多位置的城市都计算了一个价值分数并给出了如下公式：一座城市的价值分数=资源1的数量*资源1的价值+资源2 的数量*资源2的价值+资源3的数量*资源3的价值+……+资 源N 的数量*资源N 的价值。然而，SS却是首次接触这个游戏。所以它非常好奇每种资源的价值是如何得出的，于是它找...

扩展欧几里得算法，顾名思义，是求gcd(a,b)的拓展。对于ax+by=gcd(a,b)我们通过此算法求出一组满足条件的x,y那么，对于ax+by=c，若c%gcd(a,b)==0，则有且有无数组解，否则无解。//性质：对于未知数(x,y)，ax+by=k*gcd(a,b) int gcd(int a,int b){ return !b?a:gcd(b,a%b);}int...

Description　　有M个球，一开始每个球均有一个初始标号，标号范围为1～N且为整数，标号为i的球有ai个，并保证Σai = M。　　每次操作等概率取出一个球（即取出每个球的概率均为1/M），若这个球标号为k（k < N），则将它重新标号为k + 1；若这个球标号为N，则将其重标号为1。（取出球后并不将其丢弃）　　现在你需要求出，经过K次这样的操作后，每个标号的球的期望个数。I...

Description给定n序列An，从序列中选若干个数，使得乘积为完全平方数，有多少种方案。两个方案不同仅当选择的位置不同。方案数对1e9+7取模。Input第一行一个整数n,n<=100000第二行n个整数表示Ai,Ai<=70Output输出方案数对1e9+7取模后的值Sample Input41 1 1 1Sample Output1...

一、高斯消元求解线性方程组Description　　贾老二是个品学兼优的好学生，但由于智商问题，算术学得不是很好，尤其是在解方程这个方面。虽然他解决 2x=2 这样的方程游刃有余，但是对于 {x+y=3 x-y=1} 这样的方程组就束手无策了。于是他要你来帮忙。前提是一次方程组且保证在integer的范围内可以处理所有问题。Input　　第一行一个数字N（1≤N≤100）表示要求的...

刚好在复习dp，顺便把矩阵快速幂加速递推给复习了。首先矩阵乘法：f[i][j]=A[i][k]*B[k][j]我们观察斐波那契数列f[n]=f[n-1]+f[n-2]，我们构造矩阵乘法A:[ fn-2 fn-1]B:[0 1]    [1 1]C=A*B=[fn-1 fn-2+fn-1]=[fn-1 fn]所以对于C[fn-1 fn]，不难得出C=[f1 f2]*B^(...

 一、欧拉函数简单应用·定义：对于一个正整数n，它的欧拉函数的值即{p|p∈[1,n),p∈N+，p与n互质}的集合的大小，我们用φ(n)表示正整数n的欧拉函数φ(n)=n(1-1/p1)(1-1/p2)……(1-1/pk)，其中p1~pk为n的质因数附证明：https://blog.youkuaiyun.com/paxhujing/article/details/51353672·性...