Examining the Rooms - HDU 3625 第一类斯特林数

Examining the Rooms

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 1205    Accepted Submission(s): 727


Problem Description
A murder happened in the hotel. As the best detective in the town, you should examine all the N rooms of the hotel immediately. However, all the doors of the rooms are locked, and the keys are just locked in the rooms, what a trap! You know that there is exactly one key in each room, and all the possible distributions are of equal possibility. For example, if N = 3, there are 6 possible distributions, the possibility of each is 1/6. For convenience, we number the rooms from 1 to N, and the key for Room 1 is numbered Key 1, the key for Room 2 is Key 2, etc.
To examine all the rooms, you have to destroy some doors by force. But you don’t want to destroy too many, so you take the following strategy: At first, you have no keys in hand, so you randomly destroy a locked door, get into the room, examine it and fetch the key in it. Then maybe you can open another room with the new key, examine it and get the second key. Repeat this until you can’t open any new rooms. If there are still rooms un-examined, you have to randomly pick another unopened door to destroy by force, then repeat the procedure above, until all the rooms are examined.
Now you are only allowed to destroy at most K doors by force. What’s more, there lives a Very Important Person in Room 1. You are not allowed to destroy the doors of Room 1, that is, the only way to examine Room 1 is opening it with the corresponding key. You want to know what is the possibility of that you can examine all the rooms finally.
 

Input
The first line of the input contains an integer T (T ≤ 200), indicating the number of test cases. Then T cases follow. Each case contains a line with two numbers N and K. (1 < N ≤ 20, 1 ≤ K < N)
 

Output
Output one line for each case, indicating the corresponding possibility. Four digits after decimal point are preserved by rounding.
 

Sample Input
3 3 1 3 2 4 2
 

Sample Output
0.3333 0.6667 0.6250

题意:有n个房间,其中每个房间都有1把某一个房间钥匙,你最多强行打开m个房间,并且不能强行打开第一个房间。问你能够打开所有房间的概率是多少。

思路:因为每个房间都有一把钥匙,所以可以把其变成几个环,对于每一个环,只要强行打开一种的一个房间即可打开这个环中所有的房间。这里就用到了第一类斯特林数,num[i][j]表示i个房间形成j个环的情况数目。分母为n个房间形成1-n个环的情况数,分子为n个房间形成1-m个环的情况数,减去其中1号房间单独一个环的情况数。

AC代码如下:

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
int T,t,n,m;
ll num[25][25];
void init()
{
    int i,j,k;
    num[1][1]=1;
    for(i=2;i<=20;i++)
       for(j=1;j<=i;j++)
          num[i][j]=num[i-1][j]*(i-1)+num[i-1][j-1];
}
int main()
{
    int i,j,k;
    ll S=1,a,b;
    init();
    scanf("%d",&T);
    for(t=1;t<=T;t++)
    {
        a=b=0;
        scanf("%d%d",&n,&m);
        for(i=1;i<=n;i++)
           b+=num[n][i];
        for(i=1;i<=m;i++)
           a+=num[n][i];
        for(i=1;i<=m;i++)
           a-=num[n-1][i-1];
        printf("%.4f\n",1.0*a/b);
    }
}



[user@localhost ~]$ cd /home/user/Homework/web [user@localhost web]$ gdb httpd.o GNU gdb (GDB) Red Hat Enterprise Linux 7.6.1-115.el7 Copyright (C) 2013 Free Software Foundation, Inc. License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html> This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law. Type "show copying" and "show warranty" for details. This GDB was configured as "x86_64-redhat-linux-gnu". For bug reporting instructions, please see: <http://www.gnu.org/software/gdb/bugs/>... Reading symbols from /home/user/Homework/web/httpd.o...(no debugging symbols found)...done. (gdb) help List of classes of commands: aliases -- Aliases of other commands breakpoints -- Making program stop at certain points data -- Examining data files -- Specifying and examining files internals -- Maintenance commands obscure -- Obscure features running -- Running the program stack -- Examining the stack status -- Status inquiries support -- Support facilities tracepoints -- Tracing of program execution without stopping the program user-defined -- User-defined commands Type "help" followed by a class name for a list of commands in that class. ---Type <return> to continue, or q <return> to quit--- Type "help all" for the list of all commands. Type "help" followed by command name for full documentation. Type "apropos word" to search for commands related to "word". Command name abbreviations are allowed if unambiguous. (gdb) stack Undefined command: "stack". Try "help". (gdb) running Undefined command: "running". Try "help". (gdb) help List of classes of commands: aliases -- Aliases of other commands breakpoints -- Making program stop at certain points data -- Examining data files -- Specifying and examining files internals -- Maintenance commands obscure -- Obscure features running -- Running the program stack -- Examining the stack status -- Status inquiries support -- Support facilities tracepoints -- Tracing of program execution without stopping the program user-defined -- User-defined commands Type "help" followed by a class name for a list of commands in that class. ---Type <return> to continue, or q <return> to quit---q Quit (gdb) stack help Undefined command: "stack". Try "help". (gdb) step 1 The program is not being run. (gdb) running Undefined command: "running". Try "help". (gdb) exit Undefined command: "exit". Try "help". (gdb) r Starting program: /home/user/Homework/web/httpd.o /bin/bash: /home/user/Homework/web/httpd.o: Permission denied /bin/bash: line 0: exec: /home/user/Homework/web/httpd.o: cannot execute: Permission denied During startup program exited with code 126. (gdb) gdb httpd Undefined command: "gdb". Try "help". (gdb) 如何退出gdb调试
08-22
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值