12 2015 档案
摘要：求出所有线段的交点，然后利用叉乘求四边形面积即可。//// main.cpp// poj1408//// Created by 陈加寿 on 15/12/31.// Copyright (c) 2015年 chenhuan001. All rights reserved.//#include...
阅读全文
摘要：#include <iostream> #include <cmath> #include <vector> #include <stdio.h> #include <string.h> #include <stdlib.h> #include <algorithm> using namespace
阅读全文
摘要：做完了才发现，好像没有人和我的做法一样的，不过我怎么都觉得我的做法还是挺容易想的。我的做法是：把周围的方框按顺时针编号，然后对于每一条边，如果点出现在边的一侧，则把另一侧所有的点加1，这样最后统计最小值+1即可。离散化一下 O(n)//// main.cpp// poj1066//// Cre...
阅读全文
摘要：1. d维空间中一个规则几何体n刀最多能切成几块f(d,n)。f(d,n) = f(d,n-1) + f(d-1,n-1) 这不就是组合数递推公式。。。于是 f(1,n)=1+nf(2,n)=1+n+n*(n-1)/2f(3,n)=1+n+n*(n-1)/2+n*(n-1)*(n-2)/6 = (...
阅读全文
摘要：精度比acos , asin 什么的高些。ParametersyValue representing the proportion of the y-coordinate.xValue representing the proportion of the x-coordinate.If both a...
阅读全文
摘要：计算几何头疼的地方一般在于代码量大和精度问题，代码量问题只要平时注意积累模板一般就不成问题了。精度问题则不好说，有时候一个精度问题就可能成为一道题的瓶颈，简直“画龙点睛”。这些年的题目基本是朝着越来越不卡精度的方向发展了，但是也不乏一些%^&%题#$%$^，另外有些常识不管题目卡不卡，都是应该知道的...
阅读全文
摘要：#include #include #include #include #define MAX_N 100using namespace std; /////////////////////////////////////////////////////////////////////常量区c...
阅读全文
摘要：1.余弦定理 cos A=(b²+c²-a²)/2bc 2.海伦公式 假设在平面内，有一个三角形，边长分别为a、b、c，三角形的面积S可由以下公式求得： S = sqrt(s1* (s1-a) * (s1-b)*(s1-c))，其中s1表示的是三角形的周长的一半。 S = sqrt(s1* (s1-
阅读全文
摘要：/****************************** 大组合取模之:110^9，log(mod)*log(b),否则log(b) ***************/long long Mod_Mul(long long a,long long b,long long mod){ lon...
阅读全文
摘要：typedef long long ll;/********************************** 大组合数取模之lucas定理模板,1=A,M>=1//输出：返回x的范围是[1,M-1]ll GetNi(ll A,ll M){ ll rex=0,rey=0; ll td=...
阅读全文
摘要：基本的模板题，统计分子分母中p出现的次数，然后求逆元取模.//// main.cpp// fzu2020//// Created by 陈加寿 on 15/12/27.// Copyright (c) 2015年 chenhuan001. All rights reserved.//#inc...
阅读全文
摘要：Lucas 定理（证明）A、B是非负整数，p是质数。AB写成p进制：A=a[n]a[n-1]...a[0]，B=b[n]b[n-1]...b[0]。则组合数C(A,B)与C(a[n],b[n])*C(a[n-1],b[n-1])*...*C(a[0],b[0])mod p 相同即：Lucas(n,m...
阅读全文
摘要：1.错排公式n各有序的元素应有n！种不同的排列。如若一个排列式的所有的元素都不在原来的位置上，则称这个排列为错排。任给一个n，求出1,2,……,n的错排个数Dn共有多少个。递归关系式为：D(n)=(n-1)(D(n-1)+D(n-2))D(1)=0,D(2)=12.组合数递推公式c(n,m) = c...
阅读全文
摘要：/************************************* 求解x^a=b(mod c) x在[0,c-1]上解的个数模板 输入：1e9>=a,b>=1,1e9>=c>=3. 返回：调用xaeqbmodc(a,b,c),返回解的个数 复杂度: 找原根的复杂度很低，所以总的复杂度为O...
阅读全文
摘要：/************** 快速幂模板 调用:Quk_Mul(a,b,mod) 返回:a^b%mod 复杂度:当mod>10^9，log(mod)*log(b),否则log(b) ***************/ long long Mod_Mul(long long a,long long b,long long mod) { long long msum=0; while...
阅读全文
摘要：刚开始看题，想了一会想到了一种容斥的做法。复杂度O( n(3/2) )但是因为题目上说有3000组测试数据，然后吓尿。完全不敢写。 然后想别的方法。唉，最近精神有点问题，昨天从打完bc开始想到1点多，没想到什么好的方法，然后躺床上睡不着，迷迷糊糊又好像挺清醒的，大概想到了用莫比乌斯反演的一种解法，初...
阅读全文
摘要：SumdivTime Limit:1000MSMemory Limit:30000KTotal Submissions:17039Accepted:4280DescriptionConsider two natural...
阅读全文
摘要：1.扩展欧几里得求逆元typedef long long ll;//ax + by = gcd(a,b)//传入固定值a,b.放回 d=gcd(a,b), x , yvoid extendgcd(ll a,ll b,ll &d,ll &x,ll &y){ if(b==0){d=a;x=1;y=...
阅读全文
摘要：1.无成段更新 2.成段更新 (区间替换) 3.（区间增减）o(╯□╰)o 两行代码的差别。
阅读全文
摘要：ZYB's Biology简单题，直接判断替换，然后判断即可。#include #include #include #include #include #include using namespace std;char s[110],t[119];int chgs(char c){ if(c=...
阅读全文
摘要：这次题目挺有意思的，可惜D题差几分钟没交对。E题的博弈还是很基础的SG函数，让我复习了下博弈。A - Uncowed Forces基本的输入输出，一个公式输出来即可，但是要注意的是，写之前先化简变成全为整数的情况。#include #include using namespace std;int m...
阅读全文
摘要：转自:Sprague-Grundy Function-SG函数--博弈论(3)公平游戏的Sprague-Grundy定理公平游戏是一种双人游戏，在游戏中双方都有完整的信息，没有牵涉，任何状态的合法操作对双方来说都是相同的。一个公平游戏可以抽象地用一个有向无环图来表示，这个图中每个点都对应这一个状态，...
阅读全文
摘要：整数的阶：设a和n是互素的正整数，使得a^x=1（mod n）成立的最小的正整数x称为a模n的阶//求阶模板：A^x=1(mod M),调用GetJie(A,M)//输入：10^10>A,M>1//输出：无解返回-1，有解返回最小正整数x//复杂度：O(M^(0.5))long long gcd(l...
阅读全文
摘要：定理1:对正整数n，欧拉函数是小于等于n的数中与n互质的数的数目（1的欧拉函数为1） 求法（容斥）： ,其中p1, p2……pn为x的所有质因数，x是不为0的整数。φ(1)=1（唯一和1互质的数(小于等于1)就是1本身） 定理2: 也就一个正整数n的所有因子的欧拉函数之和等于n。 证明： 先来看一个
阅读全文
摘要：/************************************* ---高次同余方程模板BabyStep-GiantStep--- 输入：对于方程A^x=B(mod C),调用BabyStep(A,B,C),(0>=1; a=(a*a)%mod; } retu...
阅读全文
摘要：随机拆分，简直机智。关于过程可以看http://wenku.baidu.com/link?url=JPlP8watmyGVDdjgiLpcytC0lazh4Leg3s53WIx1_Pp_Y6DJTC8QkZZqmiDIxvgFePUzFJ1KF1G5xVVAoUZpxdw9GN-S46eVeiJ6Q...
阅读全文
摘要：根据费马小定理:对于素数n，a(0>=1; a = (a+a)%mod; } return msum;}long long Quk_Mul(long long a,long long b,long long mod){ long long qsum=1; whi...
阅读全文
摘要：裸的pell方程。 然后加个快速幂.No more tricks, Mr NanguoTime Limit: 3000/1000 MS (Java/Others)Memory Limit: 65535/32768 K (Java/Others)Total Submission(s): 320Acce...
阅读全文
摘要：佩尔方程x*x-d*y*y=1，当d不为完全平方数时，有无数个解，并且知道一个解可以推其他解。 如果d为完全平方数时，可知佩尔方程无解。假设(x0,y0)是最小正整数解。则：xn=xn-1*x0+d*yn-1*y0yn=xn-1*y0+yn-1*x0证明只需代入。 如果忘记公式可以自己用(x0*x0...
阅读全文
