码迷,mamicode.com
首页 > 其他好文 > 详细

Modified LCS (扩展欧几里德)细写了对这个算法思路的理解

时间:2014-05-23 00:26:20      阅读:377      评论:0      收藏:0      [点我收藏+]

标签:扩展欧几里德   gcd   

题目:Modified LCS

为过此题去仔细研究了下扩展欧几里德算法,领悟了一些精华。

模板为:

void gcd(ll a, ll b, ll& d, ll& x, ll& y) 
{ 
    if(!b) {d = a; x = 1; y = 0;} 
    else{ gcd(b, a%b, d, y, x); y -= x*(a/b);} 
} 
这里算出来的x,y是对于方程ax+by=gcd(a,b)而言的一组解。
为什么叫扩展欧几里德说明肯定用了欧几里德算法的原理即:gcd(a,b)=gcd(b,a%b);
ax+by=gcd(a,b)=gcd(b,a%b)=。。。可以一直使用欧几里德算法直到b=0,即就是刚才算法模板的第一句话,
if(!b) {d = a; x = 1; y = 0;} 
当b=0时,a = gcd(a,b),ax+by=gcd(a,b)这个方程就变为ax=a了,即x=1,y等于几其实是随便的,但是必须要有数字为了反推出之前的解,这里就为了简洁就等于0。PS:可以尝试不填0说不定靠点运气能提高效率(玩笑)。

然后得到了最后一组解(1,0)如何推出第一个方程的解呢?

显然可知:gcd(b,a%b) = bx+a%by。注意这里的 x和y跟第一个方程的x和y是不一样的。
连等式得:ax1 + by1 = bx2 + a%by2 = bx2 + (a - a/b*b)y2这里是利用取模的性质.即a%b = a - a/b*b。
交换位置得: ax1 + by1 = ay2 + b(x2 - a/b * y2);
因为是构造解所以根据恒等式设 x1=y2 ,y1 = (x2 - a/b * y2);
注意模板里第二句
else{ gcd(b, a%b, d, y, x); y -= x*(a/b);} 

x和y是交换位置的。就是因为 x1=y2,即上一层的 y是当前的x,当前gcd的y应该是上一层的x减去a/b×y。就是y1 = (x2 - a/b * y2);
但是在这
gcd(b, a%b, d, y, x);
之后。你上一层的x的值变成了当前y的值,上一层y的值变成了当前x的值。当前的x值就是这个,但是y值却不是上一层的x,需要计算,所以就有了
 y -= x*(a/b);
如果还有不懂可以拿笔模拟一下,先不停的递归直到终点即b==0时,然后给x,y赋值1和0,再一步一步回推执行 y -= x*(a/b);


好了。可以讲那道题了。题意就是2个等差序列,告诉你个数,起点,公差,问你有几个数是2个序列共有的?

此题就是要解F1 + D1x = F2 + D2y在0=<x<=N1-1 && 0=<y<=N2-1最多的整数解。(F为起点,D为公差,N为个数)
移项可得 D1x - D2y = F2 - F1一个标准的可用扩展欧几里德的方程。
唯一不同的是方程b的地方要用-D2,先用模板解出一个解x0,y0。

先设g = gcd(d1,d2)即最小公约数。

明确下思路,此题是要求最小且大于等于0的x和y,然后这个2个点就是起点了,其后一样的数就是起点加上2个序列公差的最小公倍数×k了。最小公倍数为d1*d2/g。
第一个序列相同数的间隔为最小公倍数除以d1,即p1 = d1*d2/g/d1 = d2/g,同理第二个序列相同数的间隔为p2 = d1/g。
故第一个序列可能的相同数个数为(N1 - 1 - x)/p1,第二个为(N1 - 1 - y)/p2。
答案为2个值的较小者。

刚才说的间隔就可以用来求方程的其他解,即x = x0 + k*p1,y = y0 + k*p2。
还要分2种情况来求x,y。
其实是2个方向,如果x0,y0有一个小于0,说明k要大于0,一直分别加p1,p2直到2者都大于等于0为止。
如果都大于0,则说明可能太大了,要求最小且大于等于0的x和y嘛,说明k要小于0。然后一直分别减p1,p2直到2者不能再减了,再减就有一项小于0为止。
具体思路已经非常清晰了,如果没有看懂希望多看几遍应该就能想通。

AC代码:
#include<cstdio> 
#include<ctype.h> 
#include<algorithm> 
#include<iostream> 
#include<cstring> 
#include<vector> 
#include<stack> 
#include<cmath> 
#include<queue> 
#include<set> 
#include<ctime> 
using namespace std; 
#define ll long long 
  
void gcd(ll a, ll b, ll& d, ll& x, ll& y) 
{ 
    if(!b) {d = a; x = 1; y = 0;} 
    else{ gcd(b, a%b, d, y, x); y -= x*(a/b);} 
} 
  
int main() 
{ 
    ll n1,f1,d1,n2,f2,d2; 
    int t; 
    scanf("%d",&t); 
    while(t--) 
    { 
        scanf("%lld%lld%lld%lld%lld%lld",&n1,&f1,&d1,&n2,&f2,&d2); 
        ll x,y,g; 
        gcd(d1,-d2,g,x,y); 
        ll f = f2-f1; 
        x *= f/g;//小心溢出 
        y *= f/g; 
        g = abs(g); 
        ll a = d1/g; 
        ll b = d2/g; 
        if(x < 0 || y < 0) 
        { 
            while(x < 0 || y < 0) 
            { 
                x += b; 
                y += a; 
            } 
        } 
        else
        { 
            while(x >= 0 && y >= 0) 
            { 
                x -= b; 
                y -= a; 
            } 
            x += b; 
            y += a; 
        } 
        ll num1 = (n1-1-x)/b; 
        ll num2 = (n2-1-y)/a; 
        printf("%lld\n",min(num1,num2)+1); 
    } 
    return 0; 
} 
  
/************************************************************** 
    Problem: 1446 
    User: HNU_TEAM_3 
    Language: C++ 
    Result: Accepted 
    Time:0 ms 
    Memory:1484 kb 
****************************************************************/


Modified LCS (扩展欧几里德)细写了对这个算法思路的理解,布布扣,bubuko.com

Modified LCS (扩展欧几里德)细写了对这个算法思路的理解

标签:扩展欧几里德   gcd   

原文地址:http://blog.csdn.net/glqac/article/details/26515207

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!