LightOJ - 1306 Solutions to an Equation
Posted 橙子用户
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了LightOJ - 1306 Solutions to an Equation相关的知识,希望对你有一定的参考价值。
题目大意:给你一个方程ax+by+c=0(x1<=x<=x2&&y1<=y<=y2)
求其整数解个数
具体思路:扩展欧几里得求一个特解,然后求一下x和y的最少变化量,
然后求出当x1<=x<=x2,与其对应的y的范围,和y1<=y<=y2做一个并,在算一下有几个解就好了
花絮:一开始我发现x=x0+xx*t(x0为特解,xx为最少变化量)y=y0-yy*t,并试图求出t的范围,wa了几发之后发现有许多+1-1的奇怪地方需要处理,于是果断换了一个做法
结果却还是wa,通过一些特殊手段发现好像是我为了把方程变为ax+by=c,我把读入的c取反,但前面的特判写的比较早~~~
扩展欧几里得是什么勒?
欧几里得算法可以求最大公约数
int gcd(int a,int b){if(b==0)return a;else return gcd(b,a%b);}
然后扩展欧几里得可以求ax+by=gcd(a,b)的一组解
怎么求的呢?
因为gcd(a,b)=gcd(b,a%b)
所以ax+by=gcd(a,b)可以变为
b*x1+(a%b)*y1=gcd(b,a%b)=ax+by
然后一顿公式变形
b*x1+(a%b)*y1=b*x1+(a-[a/b]*b)*y1([]为下取整)
=b*x1+(a-[a/b]*b)*y1=a*y1+b*(x1-[a/b]*y1)=ax+by
所以x=y1 y=x1-[a/b]*y1
于是可以写成这样
void exgcd(int a,int b,int& d,int& x,int & y) { if(b==0)d=a,x=1,y=0; else exgcd(b,a%b,d,y,x),y-=(a/b)*x; }
最后是AC代码
#include<cstdio> #include<algorithm> #include<cmath> #define int long long using namespace std; int T; int gcd(int a,int b){if(b==0)return a;else return gcd(b,a%b);} void exgcd(int a,int b,int& d,int& x,int & y) { if(b==0)d=a,x=1,y=0; else exgcd(b,a%b,d,y,x),y-=(a/b)*x; } main() { scanf("%lld",&T); int t; for (t=1;t<=T;t++) { int a,b,c,X1,X2,Y1,Y2; scanf("%lld%lld%lld%lld%lld%lld%lld",&a,&b,&c,&X1,&X2,&Y1,&Y2);c=-c; printf("Case %lld: ",t); if(a==0&&b==0) { if(c==0) printf("%lld\\n",(Y2-Y1+1)*(X2-X1+1));else puts("0"); continue; } else if(a==0) { if(c%b!=0) puts("0"); else { int z=c/b; if(z>=Y1&&z<=Y2) printf("%lld\\n",(X2-X1+1));else puts("0"); } continue; } else if(b==0) { int z=c/a; if(c%a!=0) puts("0"); else if(z>=X1&&z<=X2) printf("%lld\\n",(Y2-Y1+1)); else puts("0"); continue; } int g=gcd(a,b); if(c%g) { puts("0"); continue; } int x,y; exgcd(a,b,g,x,y); x=x*(c/g),y=y*(c/g); int xx=b/g,yy=a/g,Ymn,Ymx; X1=((x-X1)%xx+abs(xx))%abs(xx)+X1; X2=((x-X2)%xx-abs(xx))%abs(xx)+X2; if(X1>X2){printf("0\\n");continue;} Ymn=(-yy)*((X1-x)/xx)+y; Ymx=(-yy)*((X2-x)/xx)+y; if(Ymn>Ymx)swap(Ymn,Ymx); Y1=max(Y1,Ymn); Y2=min(Y2,Ymx); Y1=((y-Y1)%yy+abs(yy))%abs(yy)+Y1; Y2=((y-Y2)%yy-abs(yy))%abs(yy)+Y2; if(Y1>Y2)puts("0");else printf("%lld\\n",(Y2-Y1)/abs(yy)+1); } return 0; }
以上是关于LightOJ - 1306 Solutions to an Equation的主要内容,如果未能解决你的问题,请参考以下文章
LightOJ - 1306 Solutions to an Equation
(状压) Marriage Ceremonies (lightOJ 1011)