Yukari's Birthday
Time Limit: 12000/6000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 3270 Accepted Submission(s): 697
Problem Description
Today is Yukari's n-th birthday. Ran and Chen hold a celebration party for her. Now comes the most important part, birthday cake! But it's a big challenge for them to place n candles on the top of the cake. As Yukari has lived for such a long long time, though she herself insists that she is a 17-year-old girl.
To make the birthday cake look more beautiful, Ran and Chen decide to place them like r ≥ 1 concentric circles. They place k i candles equidistantly on the i-th circle, where k ≥ 2, 1 ≤ i ≤ r. And it's optional to place at most one candle at the center of the cake. In case that there are a lot of different pairs of r and k satisfying these restrictions, they want to minimize r × k. If there is still a tie, minimize r.
To make the birthday cake look more beautiful, Ran and Chen decide to place them like r ≥ 1 concentric circles. They place k i candles equidistantly on the i-th circle, where k ≥ 2, 1 ≤ i ≤ r. And it's optional to place at most one candle at the center of the cake. In case that there are a lot of different pairs of r and k satisfying these restrictions, they want to minimize r × k. If there is still a tie, minimize r.
Input
There are about 10,000 test cases. Process to the end of file.
Each test consists of only an integer 18 ≤ n ≤ 10 12.
Each test consists of only an integer 18 ≤ n ≤ 10 12.
Output
For each test case, output r and k.
Sample Input
18 111 1111
Sample Output
1 17 2 10 3 10
Source
2012 Asia ChangChun Regional Contest
/*枚举r,然后用二分来求对应的k.
因为题目上k>=2,所以r肯定不会很大,二分就可以了。
需要注意的是二分时初始右边界不能太大,不然会导致计算过程超出long long*/
#include<stdio.h>
#include<math.h>
typedef __int64 ll;
ll fun(ll k,ll i,ll n)
{
ll s=1,j,m;
for(j=1;j<=i+1;j++)
{
s=s*k;
if(((s-k)/(k-1)) >n)
return 1;
}
m=(s-k)/(k-1);
if(m==n)
return 0;
else if(m==n-1)
return 0;
else if(m<n)
return -1;
else if(m>n)
return 1;
}
int main()
{
ll n,l,r,ans,mid,i,k,x,y;
while(scanf("%I64d",&n)!=EOF)
{
y=1; x=n-1;//是1的时候
for(i=2;i<=40;i++)
{
l=2;r=(ll)pow(n,1.0/i);//求出右边界的最大值
while(l<=r)
{
ll mid=(l+r)/2;
//printf("mid=%lld ans=%lld\n",mid,ans);
ans=fun(mid,i,n);
if(ans==0)
{
if(mid*i<x*y)
{
x=mid;
y=i;
}
break;
}
else if(ans<0)
l=mid+1;
else
r=mid-1;
}
}
printf("%I64d %I64d\n",y,x);
}
return 0;
}