ブルーブリッジカップ-買えない数
タイトル:買えない数
明ちゃんはキャンディ屋を開いた.彼は工夫を凝らした:フルーツキャンディを4粒1パックと7粒1パックの2種類に包んだ.キャンディは分解して売ってはいけません.
子供が砂糖を買いに来たとき、彼はこの2つの包装で組み合わせた.もちろん、砂糖を10粒買うなど、組み合わせられないキャンディもあります.
コンピュータでテストしてみてください.このような包装の場合、最大で買えない数は17です.17より大きい数字はすべて4と7で組み合わせることができます.
本題の要件は,2つのパッケージの数が既知である場合に,最大で組み合わせられない数字を求めることである.
入力:
2つの正の整数は、各包装中の糖の粒数(いずれも1000未満)を表す.
要求出力:
最大で買えない糖の数を表す正の整数
無解を考慮する必要はない
例:
ユーザー入力:
4 7
プログラムは出力するべきです:
17
次に例を示します.
ユーザー入力:
3 5
プログラムは出力するべきです:
7
分析:
この問題の一般的な解法は、この2,3日拡張されたユークリッドアルゴリズムを勉強したばかりなので、このアルゴリズムを借りて性能を高めることができると思います.コードを修正したり削除したりします.
まず、理論的根拠はこの文章から来ています.http://www.cppblog.com/yuanyuelang/articles/95378.html
この文章には、a(x+qb)+b(y-qa)=c;qは任意の整数
コードでは、mはここのaに対応し、nはここのbに対応する.
この式から,x+qb>=0とすることができる.y-qa>=0 (1)
ここで、x=(c/d)*x‘y=(c/d)*y’
ここでd=1であるため、x=c*x',y=c*y'
x′とy′については拡張ユークリッドアルゴリズムで求めることができ,
式(1)は、次のように変換できます.
c*x'+qb>=0; (2)
c*y'-qa>=0; (3)
以上の2つの式を満たすと,未知変数にはcとqがあり,長い間試算されていたが,ここではcの値をさらに直接得ることができなかったため,一時的に1つ1つ検証cが条件を満たすか否かを後退するしかなかった.
y'>0の場合、x'は必ず<0であり、(2)(3)と結合し、(aをm、bをnに置換する)
入手可能:
-m*x'*q<=|x'*y'*c|<=n*y'*q (4)
ここでcの値を直接テストすることができて、cがどんな値を取って、qを解くことができません(qは整数です)
例えばm=3,n=4(x’=-1,y’=1)
(4)によると、
3 q<=c<=4 q(ここではcの値を直接導き出せばテストする必要はないが、cの値をどのように直接導き出すかはしばらく考えられない)
cが何の値を取ってqが解けないとき、cはm*nから下へ減少します.
c=12,11,10,9,8,7,6のときqはいずれも解がある
c=5のときqは無解である
だからc=5は対応する答えです
以上がy'>0の場合、y'<0の場合、対応する(4)式は、
-n*y'*q<=|x'*y'*c|<=m*x'*q (5)
例えば、m=4、n=7(x'=2、y'=-1)
(5)式により、
7q<=2c<=8q
c=28~18でqはすべて解があります
c=17のときqは解けない
したがってc=17は結果である
次に、対応するjavaコードを示します.
入力:
22 23
出力:
461
明ちゃんはキャンディ屋を開いた.彼は工夫を凝らした:フルーツキャンディを4粒1パックと7粒1パックの2種類に包んだ.キャンディは分解して売ってはいけません.
子供が砂糖を買いに来たとき、彼はこの2つの包装で組み合わせた.もちろん、砂糖を10粒買うなど、組み合わせられないキャンディもあります.
コンピュータでテストしてみてください.このような包装の場合、最大で買えない数は17です.17より大きい数字はすべて4と7で組み合わせることができます.
本題の要件は,2つのパッケージの数が既知である場合に,最大で組み合わせられない数字を求めることである.
入力:
2つの正の整数は、各包装中の糖の粒数(いずれも1000未満)を表す.
要求出力:
最大で買えない糖の数を表す正の整数
無解を考慮する必要はない
例:
ユーザー入力:
4 7
プログラムは出力するべきです:
17
次に例を示します.
ユーザー入力:
3 5
プログラムは出力するべきです:
7
分析:
この問題の一般的な解法は、この2,3日拡張されたユークリッドアルゴリズムを勉強したばかりなので、このアルゴリズムを借りて性能を高めることができると思います.コードを修正したり削除したりします.
まず、理論的根拠はこの文章から来ています.http://www.cppblog.com/yuanyuelang/articles/95378.html
この文章には、a(x+qb)+b(y-qa)=c;qは任意の整数
コードでは、mはここのaに対応し、nはここのbに対応する.
この式から,x+qb>=0とすることができる.y-qa>=0 (1)
ここで、x=(c/d)*x‘y=(c/d)*y’
ここでd=1であるため、x=c*x',y=c*y'
x′とy′については拡張ユークリッドアルゴリズムで求めることができ,
式(1)は、次のように変換できます.
c*x'+qb>=0; (2)
c*y'-qa>=0; (3)
以上の2つの式を満たすと,未知変数にはcとqがあり,長い間試算されていたが,ここではcの値をさらに直接得ることができなかったため,一時的に1つ1つ検証cが条件を満たすか否かを後退するしかなかった.
y'>0の場合、x'は必ず<0であり、(2)(3)と結合し、(aをm、bをnに置換する)
入手可能:
-m*x'*q<=|x'*y'*c|<=n*y'*q (4)
ここでcの値を直接テストすることができて、cがどんな値を取って、qを解くことができません(qは整数です)
例えばm=3,n=4(x’=-1,y’=1)
(4)によると、
3 q<=c<=4 q(ここではcの値を直接導き出せばテストする必要はないが、cの値をどのように直接導き出すかはしばらく考えられない)
cが何の値を取ってqが解けないとき、cはm*nから下へ減少します.
c=12,11,10,9,8,7,6のときqはいずれも解がある
c=5のときqは無解である
だからc=5は対応する答えです
以上がy'>0の場合、y'<0の場合、対応する(4)式は、
-n*y'*q<=|x'*y'*c|<=m*x'*q (5)
例えば、m=4、n=7(x'=2、y'=-1)
(5)式により、
7q<=2c<=8q
c=28~18でqはすべて解があります
c=17のときqは解けない
したがってc=17は結果である
次に、対応するjavaコードを示します.
import java.util.Scanner;
public class ys_cbk_08 {
//
public static void main(String[] args) {
Scanner scanner = new Scanner(System.in);
// 1,
int m = scanner.nextInt(), n = scanner.nextInt();
// 1. m*x+n*y=gcd(m,n) x,y
// , ,x,y
int gcd = gcd(m, n);
if (gcd != 1) {
System.out.println(" ");
return;
}
if (y < 0) {// y y’
int a = -n * y;
int b = m * x;
int c = 0;
for (int i = n * m; i >= 1; i--) {
c = -x * y * i; // c
if ((c / a) * b < c) {
System.out.println(i);
break;
}
}
} else {
int a = -m * x;
int b = n * y;
int c = 0;
for (int i = n * m; i >= 1; i--) {
c = -x * y * i;
if ((c / a) * b < c) {
System.out.println(i);
break;
}
}
}
}
private static int x;
private static int y;
/**
* x,y
* @param m
* @param n
* @return
*/
public static int gcd(int m, int n) {
// 1,
// if(n==0){
// x=1;
// y=0;
// return m;
// }
// int d=gcd(n,m%n);
// int tmp=x;
// x=y;
// y=tmp-m/n*y;
// return d;
// :http://www.cnblogs.com/frog112111/archive/2012/08/19/2646012.html
// 2,
int x1, y1, x0, y0;
x0 = 1;
y0 = 0;
x1 = 0;
y1 = 1;
x = 0;
y = 1;
int r = m % n;
int q = (m - r) / n;
while (r != 0) {
x = x0 - q * x1;
y = y0 - q * y1;
x0 = x1;
y0 = y1;
x1 = x;
y1 = y;
m = n;
n = r;
r = m % n;
q = (m - r) / n;
}
return n;
}
}
入力:
22 23
出力:
461