给出一个关于变量x,y的不定方程,显然这个方程可能有多个整数解。
如果有解,使得最小的一组整数解是什么。你只需要输出的最小值。
输入描述:
第一行三个空格隔开的整数a,b,c(0 ≤ a,b,c≤ 105)。
第二行两个空格隔开的整数p1,p2(1 ≤ p1,p2 ≤ 105)。
第三行两个空格隔开的整数q1,q2(1 ≤ q1,q2 ≤ 105)。
输出描述:
如果方程无整数解,输出“Kuon”。
如果有整数解,输出的最小值。
示例1
输入
2 2 1
1 1
1 1
输出
Kuon
示例2
输入
1 2 3
1 1
1 1
输出
4
思路:
其实有拓展欧几里得的影子
用GCD来判断不定方程是否有解,如果有解,将不定方程和带入方程进行化简;
过程如下:
;
将y带入方程可得:
由抛物线的性质可知:方程的最小值在对称轴附近
由上式化简可知,抛物线的对称轴为:
然后在这条抛物线两边来找满足的整数解就可以了(即满足),然后求两边满足要求的最小解就可以了
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <math.h>
using namespace std;
ll gcd(ll a,ll b)
{return b?gcd(b,a%b):a;}
int main(int argc, char const *argv[])
{
ll a,b,c;
ll p1,p2;
ll q1,q2;
cin>>a>>b>>c;
cin>>p1>>p2;
cin>>q1>>q2;
ll _=gcd(a,b);
if(c%_)
cout<<"Kuon"<<endl;
else
{
ll x=((2*a*c*q2)/(b*b)+a*q1/b-p1)/(2*(p2+(a*a*q2)/(b*b)));
ll x1=x,x2=x;
while((c-a*x1)%b)
x1--;
while((c-a*x2)%b)
x2++;
ll a1=(c-a*x1)/b;
ll a2=(c-a*x2)/b;
ll res1=(p1*x1+p2*x1*x1+q2*a1*a1+q1*a1);
ll res2=(p1*x2+p2*x2*x2+q2*a2*a2+q1*a2);
cout<<min(res1,res2)<<endl;
}
return 0;
}
其实是down大神的这是链接