题目来源:http://vjudge.net/problem/viewProblem.action?id=19592扩展欧几里得方程: ax+by=gcd(a,b)一定有解 把a=floor(x/k);b=ceil(x/k);floor,ceil分别为向下取余和向上取余。 ...
分类:
其他好文 时间:
2014-07-27 09:49:22
阅读次数:
297
Invoker
Problem Description
On of Vance's favourite hero is Invoker, Kael. As many people knows Kael can control the elements and combine them to invoke a powerful skill. Vance like Kael v...
分类:
其他好文 时间:
2014-07-26 02:16:36
阅读次数:
312
点击打开链接
Romantic
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2958 Accepted Submission(s): 1160
Problem Description
T...
分类:
其他好文 时间:
2014-07-26 02:12:16
阅读次数:
229
要死了,这个题竟然做了两天……各种奇葩的错误……HNU的12831也是这个题。题意: 给你两个等差数列,求这两个数列的公共元素的数量。 每个数列按照以下格式给出: N F D(分别表示每个数列的长度,首项,公差)。思路: 先用扩展欧几里得得到两个数列的一个交点,然后求出两个数列的第一个交点。然...
分类:
其他好文 时间:
2014-07-22 22:42:13
阅读次数:
209
一、题目POJ 1061 青蛙的约会【关于“欧几里得求最大公约数”和“扩展欧几里得算法”的题目】二、题目源程序#include using namespace std;#define LL long longLL gcd(LL a, LL b){ return b ? gcd(b, a%b)...
分类:
其他好文 时间:
2014-07-22 22:34:36
阅读次数:
253
扩展欧几里得算法就是求: ax + by = gcd(a, b)的一组整数解(x, y)一、非递归的实现:首先看a = 60, b = 22的情况:表格左边是欧几里得算法,右边等式计算ax + by = gcd(a, b)的解a = 2 × b + 1616 = a - 2bb = 1× ...
分类:
其他好文 时间:
2014-07-21 09:35:19
阅读次数:
290
将式子变形为ax-c=my可以看出原式有解当且仅当线性方程ax-my=c有解设g = gcd(a, m)则所有形如ax-my的数都是g的倍数因此如果g不整除c则原方程无解。下面假设g整除c:利用扩展欧几里得算法解出 au + mv =g 一个特解(u0, v0)所以可用整数c/g乘上上式au0*(c...
分类:
其他好文 时间:
2014-07-21 00:35:53
阅读次数:
297
刚学习的扩展欧几里得算法,刷个水题
求解 线性不定方程 和 模线性方程
求方程 ax+by=c 或 ax≡c (mod b) 的整数解
1、ax+by=gcd(a,b)的一个整数解:
void ex_gcd(int a,int b,int &d,int &x,int &y)//扩展欧几里得算法
{
if(!b){d=a;x=1;y=0;}
else {ex_gcd(...
分类:
其他好文 时间:
2014-07-19 18:26:00
阅读次数:
226
1 LL Ex_GCD(LL a,LL b,LL &x,LL& y) 2 { 3 if(b==0) 4 { 5 x=1; 6 y=0; 7 return a; 8 } 9 LL g=Ex_GCD(b,a%b,x,y);...
分类:
其他好文 时间:
2014-07-19 17:04:57
阅读次数:
249
Marbles
Input: standard input
Output: standard output
I have some (say, n) marbles (small glass balls) and I am going to buy some boxes to store them. The boxes are of two types:
Type 1: each ...
分类:
其他好文 时间:
2014-07-19 08:01:10
阅读次数:
264