ZOJ 3380 Patchouli's Spell Cards(DP,大数)

转载请注明出处,谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526       by---cxlove

题目:有m个位置,每个位置填入一个数,数的范围为1-n,问至少有i个位置的数字一样的概率为多少

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=3957 

完全就是仿造zeroclock的做法的。。。主要是练习JAVA

至少有i个位置一样,也就是i+1,i+2,……m个位置一样,不好处理

我们知道总的方案数,n^m也就是分母

可以dp处理有没有i个位置的数字一样的概率

dp[i][j]表示当前已经放了j个位置,用到了第i种颜色

即dp[i][j]=dp[i-1][j-k]*C(m-(j-k),k)    (k<=j&&k<l) 表示我们找出k个位置填入数字i

因为如果l>m则必然是无解的,l>m/2可以直接用组合计数算出,时间减少一半

数字很大,也没有取模,需要大数处理


import java.util.*;
import java.math.*;
public class Main
{
	static BigInteger[][] dp=new BigInteger[105][105];
	static BigInteger[][] c=new BigInteger[105][105];
	public static void main(String[] args)
	{
		Scanner cin = new Scanner(System.in);
		for(int i=0;i<105;i++)
		{
			c[i][0]=c[i][i]=BigInteger.ONE;
			for(int j=1;j<i;j++)
			{
				c[i][j]=c[i-1][j-1].add(c[i-1][j]);
			}
		}
		int m,n,l;
		while(cin.hasNext())
		{
			m=cin.nextInt();
			n=cin.nextInt();
			l=cin.nextInt();
			BigInteger total=BigInteger.valueOf(n).pow(m);
			if(l>m)
			{
				System.out.println("mukyu~");
				continue;
			}
			if(l>m/2)
			{      
                BigInteger ans = BigInteger.ZERO;  
                for (int i = l; i <= m; ++i)  
                    ans = ans.add(c[m][i].multiply(BigInteger.valueOf(n-1).pow(m-i)));  
                ans = ans.multiply(BigInteger.valueOf(n));  
                BigInteger gcd = ans.gcd(total);  
                System.out.println(ans.divide(gcd)+"/"+total.divide(gcd));  
                continue;
            }  
			for(int i=0;i<=n;i++)
			{
				for(int j=0;j<=m;j++)
				{
					dp[i][j]=BigInteger.ZERO;
				}
			}
			dp[0][0]=BigInteger.ONE;
			for(int i=1;i<=n;i++)
			{
				for(int j=1;j<=m;j++)
				{
					for(int k=0;k<l&&k<=j;k++)
					{
						dp[i][j]=dp[i][j].add(dp[i-1][j-k].multiply(c[m-(j-k)][k]));
					}
				}
			}
			BigInteger ans=new BigInteger("0");
			for(int i=1;i<=n;i++)
			{
				ans=ans.add(dp[i][m]);
			}
			ans=total.subtract(ans);
			BigInteger gcd=ans.gcd(total);
			System.out.println(ans.divide(gcd)+"/" + total.divide(gcd));
		}
	}
}


评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值