uva 10673

來源:互聯網
上載者:User

ax+by=gcd(a,b); 為方便 g=gcd(a,b), ax+by=g

這個證明有些複雜就不寫了,而如何構造一個可行解(x1,y1)其實也在證明過程中

在得到一個可行解後就可以得到無數組解,他們是(x1-k*(b/g) , y1+k*(a/g)) , (其中g=gcd(a,b),k是整數)

而對於方程ax+by=c,只要c是g倍數那麼就有整數解,否則沒有

看完原題,p是x/k的下整,q是x/k的上整,然後p*m+q*n=x,這個方程其實就是ax+by=c的形式,而且這個方程一定有整數解

因為d=gcd(p,q),若p=q,d=p,若p!=q即|p-q|=1,則d=1

所以無論如何x一定是d的倍數,這個方程一定有整數解

#include<iostream>#include<cstdio>#include<string.h>#include<math.h>#include<stdlib.h>using namespace std;int x,y,d;void exgcd(int a,int b){    if(b==0)    {        x=1;        y=0;        d=a;    }    else    {        exgcd(b,a%b);        int t=x;        x=y;        y=t-a/b*x;    }}int main(){    int t,s,k,a,b;    double xx,kk;    scanf("%d",&t);    while(t--)    {        scanf("%d%d",&s,&k);        xx=s;        kk=k;        a=floor(xx/kk);        b=ceil(xx/kk);        exgcd(a,b);        x=s/d*x;        y=s/d*y;        printf("%d %d\n",x,y);    }    return 0;}


相關文章

聯繫我們

該頁面正文內容均來源於網絡整理,並不代表阿里雲官方的觀點,該頁面所提到的產品和服務也與阿里云無關,如果該頁面內容對您造成了困擾,歡迎寫郵件給我們,收到郵件我們將在5個工作日內處理。

如果您發現本社區中有涉嫌抄襲的內容,歡迎發送郵件至: info-contact@alibabacloud.com 進行舉報並提供相關證據,工作人員會在 5 個工作天內聯絡您,一經查實,本站將立刻刪除涉嫌侵權內容。

A Free Trial That Lets You Build Big!

Start building with 50+ products and up to 12 months usage for Elastic Compute Service

  • Sales Support

    1 on 1 presale consultation

  • After-Sales Support

    24/7 Technical Support 6 Free Tickets per Quarter Faster Response

  • Alibaba Cloud offers highly flexible support services tailored to meet your exact needs.