模乘(解决乘法取模爆long long) 二进制思想,变乘法为多次加法,具体思想跟着代码手算一遍就理解了,挺简单的 ll qmul(ll a,ll b,ll m) { ll ans=0; while(b){ if(b&1) ans=(ans+a)%m; a=(a+a)%m; b=b>>1; } return ans; } //快速幂中的ans*a和a*a可能会爆LL,使用模乘解决 ll qpow(ll a,ll b,ll m) { ll ans=1; while(b){ if(b&am

题目链接:http://hihocoder.com/problemset/problem/1388?sid=974337 题目大意:找出一个$k$,使得$\sum_{i=0}^{n-1}(A_{i}-B_{i+k \quad mod \quad n})^{2}$最小 把那个式子拆开..得到: $\sum _{i=0}^{n-1} A_{i}^{2}+\sum_{i=0}^{n-1}B_{i}^{2}-\sum_{i=0}^{n-1}2*A_{i}*B_{i+k}$ 前面的两个和式是定值,所以最小

HDU 4602 Partition f[i]表示和为i的方案数.已知f[i]=2i-1. dp[i]表示和为i,k有多少个.那么dp[i]=dp[1]+dp[2]+...+dp[i-1]+f[i-k]. 考虑与f有关的项: f[n-k]是答案的一部分,即2n-k-1是答案的一部分. 把与dp有关的项: 令s[i-1]=dp[1]+dp[2]+...+dp[i-1],那么s[n-1]是答案的一部分. s[i]=s[i-1]+dp[i],又dp[i]=s[i-1]+f[i-k].推出s[i]=2*

题目描述 众所周知,Hzwer学长是一名高富帅,他打算投入巨资发展一些小城市. Hzwer打算在城市中开N个宾馆,由于Hzwer非常壕,所以宾馆必须建在空中,但是这样就必须建立宾馆之间的连接通道.机智的Hzwer在宾馆中修建了N-1条隧道,也就是说,宾馆和隧道形成了一个树形结构. Hzwer有时候会花一天时间去视察某个城市,当来到一个城市之后,Hzwer会分析这些宾馆的顾客情况.对于每个顾客,Hzwer用三个数值描述他:(S, T, V)表示该顾客这天想要从宾馆S走到宾馆T,他的速度是V. Hz

1.首先需要掌握二进制数的一种特性,00=0,01=1,10=2,11=3.每一个二进制的值代表他前面的二进数的个数,比如11=3,他的前面就有三个二进制的数字,不过在本题中,题目数据是1-n,故把0抛弃掉,答案就变成了3-1=2个,但还要加上自身的话,就还是二进制的值.比如10=2,前面有01和10.11=3,前面有01 10 11.以此类推.2.那么问题就转化成了要求出小于等于给出十进制数字的最大二进制数,求出他的值即可.比如123,最大就是111,结果就是7.3. 问题转化成求最大的二进制

1.当调试时发现无法正常调用函数时,检查是否发生爆栈 对于每个栈仅有4MB的空间,开int只能开大约5*10^5. 大数组一定要开全局变量 2.当long long=int*int时会爆int,一定要注意 在进行long long计算时,最好全部开成long long ,或加上1LL 3.vector.size()的类型为unsigned int 当用其做减法时一定要注意边界条件,防止差小于0溢出 4.在求特殊最值时,一定要特殊定制初始化值,不能套用平常初始值,防止初始值不够. Ex:1<<3

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5446 C(n, m) % (p1*p2*p3*...*pk)的值 其实这个就是中国剩余定理最后算出结果后的最后一步求余 那C(n, m)相当于以前我们需要用中国剩余定理求的值 然而C(n, m)太大,我们只好先算出 C(n, m) % p1 = r1 C(n, m) % p2 = r2 C(n, m) % p3 = r3 . . . C(n, m) % pk = rk 用Lucas,这些r1,r2,

2018.11.21: 1.[BZOJ 4868][SHOI 2017] 从后往前枚举最后位置即可,如果$A<B$,用尽可能多的$A$替换$B$操作 Tip:很大的$C$可能爆$longlong$,注意特判掉与$C$相乘超过上限的数! 2.[BZOJ 4872][SHOI 2017] 首先每个灯最多按一次,接下来发现正确按取集合唯一,只不过顺序任意 设$dp[i]$表示从正确按取$i$个到$i-1$个的期望次数,则有: $dp[i]=\frac{i}{n}+\frac{n-i}{n}*(dp[i

bzoj1009:kmp想法+递推+矩阵快速幂.很好的想法,考虑用长串去kmp匹配短串,dp[i][j]表示匹配指针分别指在i.j位置时候,前i位母字符串一共有多少种可能性,那么dp[i][j]=Σdp[i-1][k]*p[k][j] p[k][j]就是状态k后加一位数字转移到状态j一种多少种可能(通过自己匹配自己来预处理p),那么可以看做一个矩阵dp[i-1]和p[][]不断相乘,矩阵快速幂即可 bzoj2120:莫对算法.最基础的莫队,O(1)转移的,按照block[l]为第一关键字,r为第

联赛除去今天刚好只有一个星期了,最后一个星期也很关键,要吃好睡好保持心情愉悦.当然也免不了最后的复习计划. 首先是模板,之前还有很多模板没有复习到,这些东西是一定要落实到位的. 每天往后面写一点...一定要写完...都是很基础的板子了... [x] 11.3 [模板]KMP字符串匹配 [ ] 11.3 [模板]左偏树(可并堆) [x] 11.3 [模板]最近公共祖先(LCA) [x] 11.3 [模板]最长公共子序列 [x] 11.4 [模板]网络最大流 [x] 11.4 [模板]最小费用最大流

题意: 给你一个1e9-1e14的质数P,让你找出这个质数的前一个质数Q,然后计算Q!mod P 题解: 1e14的数据范围pass掉一切素数筛法,考虑Miller-Rabin算法. 米勒拉宾算法是一种判断素数的随机化算法,由于其随机性,它不能保证总是正确的,但其对于一个素数,总会返回素数的结果,对于一个合数,才有极小概率返回素数的结果(假阳性). 米勒拉宾算法对于单个素数的判断时间复杂度为$O(log^3n)$.(因为1e14相乘会爆longlong,模乘要写成龟速乘,因此要多一个log) 1

bzoj 4386: [POI2015]Wycieczki 这题什么素质,爆long long就算了,连int128都爆……最后还是用long double卡过的……而且可能是我本身自带大常数吧,T了好长时间…… 先说一下超级汇点的计数吧,先说结论: 1.将所有点(此题中只有一级点)向一个超级汇点0连边,将矩阵乘n次,相应的f[i][j]即为从i到j的走n步方案数,f[i][0]为i到0走n步的方案数,若在给他乘一个ans矩阵(ans在前),则f[0][0]-n(点数)为所有长度等于n(指数)的

给每个国家建一个链表,这样分治过程中的复杂度就和序列长度线形相关了,无脑套整体二分就可以. (最坑的地方是如果所有位置都是一个国家,那么它的样本个数会爆longlong!!被这个坑了一次,大于p[i]的时候break就行了). #include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<vector> #define N 300005 #de

4386: [POI2015]Wycieczki Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 197  Solved: 49[Submit][Status][Discuss] Description 给定一张n个点m条边的带权有向图,每条边的边权只可能是1,2,3中的一种.将所有可能的路径按路径长度排序,请输出第k小的路径的长度,注意路径不一定是简单路径,即可以重复走同一个点. Input 第一行包含三个整数n,m,k(1<=n<=40,1&

http://www.lydsy.com/JudgeOnline/problem.php?id=2154 题意:求$\sum_{i=1}^{n} \sum_{j=1}^{m} lcm(i, j)$, $n,m<=1e7$ #include <bits/stdc++.h> using namespace std; typedef long long ll; const int N=1e7+10, MD=20101009; int p[N], pcnt, n, m; bool np[N];

题意:求A^B的所有因子之和 很容易知道,先把分解得到,那么得到,那么 的所有因子和的表达式如下 第一种做法是分治求等比数列的和  用递归二分求等比数列1+pi+pi^2+pi^3+...+pi^n: (1)若n为奇数,一共有偶数项,则:      1 + p + p^2 + p^3 +...+ p^n = (1+p^(n/2+1)) + p * (1+p^(n/2+1)) +...+ p^(n/2) * (1+p^(n/2+1))      = (1 + p + p^2 +...+ p^(n/

题不是很难,但是啊,人很傻啊...机子也很鬼畜啊... 3231: [Sdoi2008]递归数列 Time Limit: 1 Sec Memory Limit: 256 MB Submit: 569 Solved: 241 [Submit][Status][Discuss] Description 一个由自然数组成的数列按下式定义: 对于i <= k:ai = bi 对于i > k: ai = c1ai-1 + c2ai-2 + - + ckai-k 其中bj和 cj (1<=j<

题意:有m个石头围成一圈,编号分别为0到m-1,现在有n只青蛙,都在0号石头上,第i只青蛙会从当前编号为p的石头跳到编号为(p+ai)%m的石头上.被青蛙经过的石头都会被占领,求这m块石头中所有被占领过的石头的编号和. 题解:对于第i只青蛙,它所能跳到的最小的位置是gcd(ai, m) 设最小位置为z,需要跳x圈,跳了y步,可得方程:x*m+z=ai*y 即:x*m-ai*y = z  由扩展欧几里得定理可知,z的最小整数解为gcd(m,ai) 因为对于单独的每一只青蛙计算结果会重复计算,所以利

P<=10一开始是吓死我了 后来想到这就是一个经典的决策单调性解决1d1d动态规划的题目 像决策单调性完全可以打表找规律,这里有一篇严谨的证明https://www.byvoid.com/blog/noi-2009-poet 关于1d1d动归的优化可以看<1d1d动态规划优化初步> 注意可能会爆longlong,所以用extended计算 type node=record l,r,x:longint; end; ..] of node; f:..] of extended; s:..]

2016-06-01 16:28:25 题目链接: Xn数列 (Codevs No.1281) 题目大意: 给定一种递推式为 Xn=(A*Xn-1+C)%M 的数列,求特定的某一项%G 解法: 矩阵乘法 不会的去看看高中矩阵的那本选修,起码知道都是啥意思,好理解得多 矩阵构造:           向量构造: A C                   X0 0  1                    1 需要注意的地方: 1.超大整数乘法,写个快速乘,防止爆longlong 2.函数的代值

Acperience 题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5734 Description Deep neural networks (DNN) have shown significant improvements in several application domains including computer vision and speech recognition. In computer vision, a particul