风城

HDOJ4549问题递推比较最后的形式为： 其中F（n）为斐波那契数列先用矩阵快速幂求斐波那契数的通项F（n）与F（n-1），由于他们在指数位置上，算出来非常大，所以需要用到欧拉降幂公式，其中（mod）=mod-1，1e9+7是个素数。 最后，用快速幂把F（n）结果算出来即可 AC代码：/***************************************************

计算几何的圆的模板题+二分判断题目中给了n个圆，要求：在这n个圆中取某个圆的圆心，然后找到最小的半径，使得这个新的圆与这n个圆的相交面积，会不小于这n个圆的面积的一半说起来，很绕口，但是因为n不大，是可以枚举的！对于每个圆心，我们都求一个最小的半径，然后n个值中间取最小的就是答案那么，如何求得这个最小呢？很简单，化计算性问题为判定性问题我们二分半径，然后验证这个半径会不会满足题中面积的条件代码如下：

HDU1576A/B求A/B（mod9973）,并且已经知道gcd(B,9973) = 1，所以 A/B(mod9973)=A*B^（-1）(mod9973)，重在求乘法逆元：由于（B，9973）=1,利用费马小定理最简单，逆元即为B^(mod-2) AC代码：/*************************************************************** By

HDU5515 首先这是一道阅读理解题，建议大家一定要好好地读题，不忽略每一个细节。 推荐大家读读这个博客题目翻译 思路比较简单，考虑fighter与speeder相遇在#2到#3的线段和相遇在#3和#4的线段的两种情况，分别卡列方程卡限制条件 一元二次方程最好能够用二分求解，精度好并且能够直接二分解x的范围 下面的解法直接对方程进行求解，过程中注意精度:AC代码，略丑陋:#inclu

HDU3697对课程结束时间由小到大进行排序，之后直接暴力枚举。 AC代码:#include <bits/stdc++.h>using namespace std;const int maxn=300+5;struct node{ int start,end;};node a[maxn];bool vis[maxn];bool cmp(node a,node b){

HDU3694 前面说了一大推都没什么用，最后是求四边形的一个fermat点。 多边形的fermat点是欧拉平面上的一个点到多边形每个点的距离之和最小的点 四边形的fermat点在四个顶点与对角线交点值5个点中较小的那个，直接套模板计算取最小值即可 在判断多边形的对角线时小小用了一下极角排序 AC代码：#include <iostream>#include <cstdio>#inclu

题目链接 统计在大区间[A,B]中与N互质的数的个数，其中A，B都在10的15次方级。暴力肯定不行，正难则反，思考求[A,B]中不与N互质的数的个数，那么该数必然是N的其中一个因子的倍数，所以首先对N的所有因子进行打表，然后根据每一个因子求在[A,B]区间中其倍数的个数。 在不同的因子进行计数的时候可能产生重复，所以进行容斥。 最后，小技巧将区间转化成[1,B]-[1,A-1]求更加方便。

HDU2276 题目是一个常见的一路路灯开开关关的问题，要求1代表开灯，0代表关灯，整个路灯围成一个圆，并且当该路灯的左边的路灯是开着的，即1时，状态翻转，其余情况状态不变。 分析状态变化: 得到转台转移矩阵后，矩阵快速幂套路即可 AC代码：/*************************************************************** By: Xingxi

HDU2256 直接去求double会损失精度，肯定不行，所以要构造递推关系去求解 最后注意一下负数取模的处理 AC代码:/************************************************ * By: Xingxing * * Date: 2016-09-26

Beauty ContestTime Limit: 3000MS Memory Limit: 65536KTotal Submissions: 35550 Accepted: 11012DescriptionBessie, Farmer John's prize cow, has just won first place

HDU2669 首先题目描述和题目本身卵关系没有 就是求解不定方程的解的解，注意x的值必须大于零。 在模版e_gcd求得特解之后，需要知道该不定方程的通解: 利用它最小的正整数x AC代码:/*************************************************************** By: Xingxing* Date: 2016-09-27*

HDU4609 这道三个智障的题目把我做成了智障，先膜拜bin神，他的题解已经很详细 Bin神的HDU4609题解 这里稍作一些解释: num[i]num[i]数组记录的不同长度的strick的个数，num[i]∗num[i]num[i]*num[i]后代表任意可重复地那两根构成的长度的种类数 ,题目要求是不可重复的，所以要进行第一次去重之后考虑num[i]num[i]的前缀和数组sum[

HDU5894 问题是求n个不同座位围成一个圆，m个孩子至少间隔k个座位坐，一共有多少种方法。 首先，第一个孩子做，有n种方法。 第一个孩子坐下后，剩下的能够坐的座位有n-1-k*m,将间隔的椅子拿走 剩下的m-1个孩子们一共还有，又因为这是一个圆排列，所以要除以孩子的总数m，最终的结果为中方法。上面仔细思考是用意确定的，接下来是题目的关键，大组合数取模 第一种，利用逆元：因为mod=1e9

HDOJ2842 Chinese Rings 自己手动推一下前几项，找到递推规律 得到 从而得到状态转移矩阵 利用矩阵快速幂求解即可 AC代码：/************************************************************************* > File Name: hdu2842.cpp > Author: X

题目链接：HDU5895 注意到题目的数据非常之大，并且递归之后求的指数幂非常之大，所以要先知道欧拉降幂公式: ,相当于想成了指数的循环节。 根据递推关系得到： 建立矩阵: 用矩阵快速幂求解出g(n*y)，之后再利用快速幂取模求出最后的结果由于数据很大，过程中把数据开大，及时取模。AC代码：#include <cstdio>#include <cstring>#include<

HDU2604Queuing 题目大意求n阶序列，每一位可以填f或者m，求不存在fff或者fmf的字串的个数，根据题目建立递推关系式: 之后用矩阵快速幂套路掉 AC代码:#include <iostream>#include <cstdio>#include <cstring>using namespace std;const int N=4;int l,mod

HDU2604Queuing 问题是求，n没有给定范围，应该是无限大，但是题目的精度卡在1e-5是一个世纪难题，但是该问题存在极限PI26\frac{PI^2}{6},这里有它详细的证明，很NB，巴塞尔问题极限证明打表确定n的上限：/************************************************************************* > Fil

HDU1061 求N^N（mod10）的值，最简单就是mod10乘 在这里试了一下快速幂的模版，注意数据开long long 或者__int64才能存下/************************************************************************* > File Name: HDU1061.cpp > Author:Xingxi

套矩阵快速幂模版，中间步骤及时mod即可/************************************************************************* > File Name: HDU1575.CPP > Author: Xingxing > Created Time: 2016/9/19 23:30:31 **************

矩阵快速幂模版题，注意k<10的初值即可。#include <iostream>#include <cstdio>#include <cmath>#include <cstring>const int N=10;int k,mod,a[10];using namespace std;struct mat{ int m[N][N]; mat() {} mat un

题目链接 题目不难看懂，首先把公式变形: ∑i=0n−1Ai+B(i+k)modn2=∑i=0n−1A2i+∑i=0n−1B2(i+k)modn−2∑i=0n−1Ai∗B(i+k)modn\sum_{i=0}^{n-1} {A_i+B_{(i+k)modn}}^2 =\sum_{i=0}^{n-1}A_i^2+\sum_{i=0}^{n-1}B_{(i+k)modn}^2-2\sum_{i=0