Note: This article is available in Chinese only. 本文暂无英文版本。
View original
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=43 题意：其实就是给了两个式子。。。(N+1)^h=a,N^h=b,a,b已知，然后求关于N的两个式子.。。 思路：数学上这个方程貌似不可解。。？ 所以只能枚举一下==。。。注意精度问题把。。。
然后用换底公式求对数的时候要向上取整。
还有b为1的时候是特殊数据。
1/* ***********************************************
2Author :111qqz
3Created Time :2016年01月28日 星期四 16时46分38秒
4File Name :code/uva/107.cpp
5************************************************ */
6
7#include <cstdio>
8#include <cstring>
9#include <iostream>
10#include <algorithm>
11#include <vector>
12#include <queue>
13#include <set>
14#include <map>
15#include <string>
16#include <cmath>
17#include <cstdlib>
18#include <ctime>
19#define fst first
20#define sec second
21#define lson l,m,rt<<1
22#define rson m+1,r,rt<<1|1
23#define ms(a,x) memset(a,x,sizeof(a))
24typedef long long LL;
25#define pi pair < int ,int >
26#define MP make_pair
27
28using namespace std;
29const double eps = 1E-12;
30const int dx4[4]={1,0,0,-1};
31const int dy4[4]={0,-1,1,0};
32const int inf = 0x3f3f3f3f;
33int a,b;
34int main()
35{
36 #ifndef ONLINE_JUDGE
37 freopen("code/in.txt","r",stdin);
38 #endif
39
40 while (scanf("%d %d",&a,&b)!=EOF)
41 {
42 if (a==0&&b==0) break;
43 int n ,h;
44 int ans1,ans2=0;
45 //b为1是特殊数据，除数为0了。。。。
46 if (b==1)
47 {
48 h = ceil(log(a)/log(2));
49 cout<<h<<" "<<2*a-1<<endl;
50 continue;
51 }
52
53 for ( int i = 1 ; ;i++)
54 {
55 if (fabs(log(a)*log(i)-log(b)*log(i+1))<eps)
56 {
57 n = i;
58 break;
59 }
60 }
61 h = (ceil)(log(b)/log(n)); //应该向上取整。。。可以看做一般规律。。换地公式求对数的时候。。想想为什么。
62 // cout<<"h:"<<h<<endl;
63 ans1 = (int)((1-b)/(1-n));
64 ans2=0;
65 int curh = 1;
66 for ( int i = 0 ; i <= h ;i++)
67 {
68 // cout<<"curh:"<<curh<<" b:"<<b<<endl;
69 ans2+=curh*b;
70 curh*=(n+1);
71 b/=n;
72 }
73 // cout<<ans1<<" "<<ans2<<endl;
74 printf("%d %d\n",ans1,ans2);
75
76
77 }
78
79 #ifndef ONLINE_JUDGE
80 fclose(stdin);
81 #endif
82 return 0;
83}