uva 10169 Urn-ball Probabilities !

概率的公式很容易得出，一般求这种至少1个的概率p就是先求┐p。

小数点后有多少个0可以考虑，a * (1/a) = 1 把a和(1/a)都化为标准型即为 a化为 b.cdef*10^k (其中1 <= b <10 bcdefghi代表每一位数)  。(1/a) 化为 0.ghi*10^k1 。b.cdef* 0.ghi = 1 所以 k+k2=0 所以1/a 的0的个数和a 的位数-1是相等的。

所以考虑当前a有多少位， a = 10^log 10(a) 所以只要计算得到log 10(a) 的整数部分即可（也就是10的多少次方）。

#include <iostream>
#include<stdio.h>
#include<cmath>
#define maxn 1000000
using namespace std;

long long ans[maxn];
double a[maxn];
long long b[maxn];

void biao()
{
   double p1=1;
   double p2 = 0;
   for(int i = 1;i <= maxn;i++)
   {
       p1=p1*((double)i*(i+1)-1)/((double)i*(i+1));
       p2 += log10((double)i*i+i);
       a[i] = 1 - p1;
       b[i] = (int)(p2 - fmod(p2,1)+0.5);
   }
}

int main()
{
    biao();
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        printf("%.6lf %lld\n",a[n],b[n]);
    }
    return 0;
}