【数学】
文章平均质量分 74
chm517
这个作者很懒,什么都没留下…
展开
-
Hdu 4611 数论
Hdu 4611 数论(输出答案注意%I64d而不是%lld)x≡a(mod A), x≡b(mod B) ……*1 nA+a=x, mB+b=x; ……*2 nA+a=mB+b ……*3 |nA-mB|=|a-b| ……*4 诡异,跟x没有关系!!! 其实关系是x∈[nA,(n+1)A-1] &&x∈[nA,(n+1)A原创 2013-07-26 18:59:27 · 502 阅读 · 0 评论 -
14年省赛K题:www
#define N 100005int vis[N],num[N];const int mod=1000000007;int o;int n,m,g,A,B;int a[105],b[105];LL ans[105];vector aa;map ma;int dp[105][8000];LL pow(LL a,LL n,LL mod){ LL res=1; whi原创 2014-10-05 19:54:42 · 693 阅读 · 0 评论 -
##数论模版##
/**************以下为常用算法************///把n的约数的莫比乌斯函数值用map的形式返回O(sqrt(n))map<int,int > moebius (int n){ map<int,int > res; vector<int > pri; for (int i=2;i*i<=n;i++) ...原创 2014-09-16 13:48:38 · 474 阅读 · 0 评论 -
uva 11139 格点计数问题
#include#include#include#include#includeusing namespace std;typedef long long LL;int n;#define N 125int gc[N][N],sum1[N][N],sum2[N][N];LL dp[N][N];LL fff(int x,int y){ if (x>y) swap(x原创 2014-09-27 21:43:09 · 1032 阅读 · 0 评论 -
##位运算若干应用##
//把右起第一个0变成1 | (100101111->100111111) | x or (x+1)//把右边连续的0变成1 | (11011000->11011111) | x or (x-1)//取右边连续的1 | (100101111->1111) | (x xor (x+1)) shr 1//取右起的第一个1 |原创 2014-09-09 22:12:40 · 361 阅读 · 0 评论 -
POJ 2115 ax≡c(mod b)求最小非负x
//ax≡c(mod b),已知a,b,c,求最小非负x解(记为ans)//Solution://存在唯一(x,y)ax+by=gcd(a,b);//令gcd(a,b)=d,调用exgcd求x,y,d;//if (c%d!=0) 无解!//a1=a/d;//b1=b/d;//c1=c/d;//ans=(x*c1%b1+b1)%b1;原创 2014-09-15 18:19:34 · 517 阅读 · 0 评论 -
Hdu 1573 线性同余方程组
#include #include #include #include #include #include #include #include #include #define pb push_back#define mp make_pair#define eps 1e-9#define zero(x) (fabs(x)<eps)#define pi acos(-1.0)原创 2014-09-12 14:45:57 · 348 阅读 · 0 评论 -
hdu 4919 数论+大数
1、扩栈语句2、将递归转为非递归求解(标称这么做的)//f[n]=4f[k]+6k (n==2k+1)// =2f[k]+2[k-1]+4k-4 (n==2k)BigInteger contribution(const BigInteger &n, const BigInteger &coefficient) { BigInteger result; i原创 2014-08-06 14:53:38 · 528 阅读 · 0 评论 -
##高斯消元##
参考 http://www.cnblogs.com/kuangbin/archive/2012/09/01/2667044.html//模板1:整数方程求整数解#include#include#include#include#includeusing namespace std;const int MAXN=50;int a[MAXN][MAXN];//增广矩阵int x[M原创 2014-03-04 21:00:01 · 494 阅读 · 0 评论 -
LA 2995
大白书 P12 #include#include#include#define for(i,t) for (int i=0;i<(t);++i)using namespace std;#define N 10char v[6][N][N],b[N][N][N];int n;//z//| x//| ///|///-----yvoid get(int k,int i,i原创 2014-01-26 21:50:38 · 447 阅读 · 0 评论 -
【转】线性筛素数,线性筛欧拉函数,求前N个数的约数个数
转自http://blog.csdn.net/ji414341055/article/details/5771066转载 2014-07-16 09:32:23 · 360 阅读 · 0 评论 -
##单纯形法##
/* * Description: * 线性规划 单纯形法实现 * Usage: * MAXN:最大变量个数 * MAXE:最大不等式个数 * n=变量个数,m=不等式个数 * 目标函数为 sum{x[i]*A[0][i];0<=i<n}+A[0][n] * 约束条件为 * 1) x[i]>=0 for 0<=i<n * 2) sum{x[i]*A[原创 2015-04-15 17:15:35 · 626 阅读 · 0 评论