[BestCoder] Round #7

1001

http://acm.hdu.edu.cn/showproblem.php?pid=4985

输出循环节

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <stdlib.h>
#include <cmath>
#include <iomanip>
#include <vector>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <cctype>
#define rd(x) scanf("%d",&x)
#define rd2(x,y)  scanf("%d%d",&x,&y)
#define rd3(x,y,z) scanf("%d%d%d",&x,&y,&z)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
int dx[8]={0,0,-1,1,1,1,-1,-1};
int dy[8]={1,-1,0,0,-1,1,-1,1};
const int maxn=1e5+2;
bool vis[maxn];
int p[maxn];
int n;

int main()
{
    while(rd(n)!=EOF)
    {
        for(int i=1;i<=n;i++)
            rd(p[i]);
        memset(vis,0,sizeof(vis));
        int position;
        bool first;
        for(int i=1;i<=n;i++)
        {
            position=i;
            first=1;
            while(!vis[position])
            {
                vis[position]=1;
                if(first)
                {
                    cout<<"("<<position;
                    first=0;
                }
                else
                    cout<<" "<<position;
                position=p[position];
                if(vis[position])
                    cout<<")";

            }
        }
        cout<<endl;
    }
    return 0;
}

1002

http://acm.hdu.edu.cn/showproblem.php?pid=4986

题意为在n个盒子里随机放着一个钥匙,且该钥匙和某一个盒子相对应，用1个魔法可以打开一个盒子，问把这些盒子全部打开，使用魔法个数的期望.

设p[i]为打开i个盒子的期望，那么另外一个盒要么是单独的一个循环节，那么插入到前面中去。p[i]=  (p[i-1]+1)/i + p[i-1]  *(i-1)/i  得p[i]=p[i-1]+ 1/i

也就是p[n]=1+ 1/2 +1/3 + .....1/n

这个数组是发散的，所以没有求和公式，只有一个近似的求解方法：
   1+1/2+1/3+......+1/n ≈ lnn+C(C=0.57722......一个无理数,称作欧拉初始,专为调和级数所用)当n很大时，有：1+1/2+1/3+1/4+1/5+1/6+...1/n = 0.57721566490153286060651209 + ln(n)//C++里面用log(n)，pascal里面用ln(n)  
0.57721566490153286060651209叫做欧拉常数  

#include <iostream>
#include <stdio.h>
#include <cmath>
#include <map>
#include <stack>
#include <queue>
#include <cctype>
#define rd(x) scanf("%d",&x)
#define rd2(x,y)  scanf("%d%d",&x,&y)
#define rd3(x,y,z) scanf("%d%d%d",&x,&y,&z)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const double C=0.57721566490153286060651209;//欧拉常数
const int maxn=100000;
double p[maxn];
int n;

void init()
{
    p[0]=0;
    for(int i=1;i<=maxn;i++)
        p[i]=p[i-1]+(1.0/i);
}

int main()
{
    init();
    while(rd(n)!=EOF)
    {
        if(n<=maxn)
            printf("%.4lf\n",p[n]);
        else
            printf("%.4lf\n",log(n)+C);

    }
    return 0;
}