1. 程式人生 > >●BZOJ 1416 [NOI2006]神奇的口袋

●BZOJ 1416 [NOI2006]神奇的口袋

bzoj pos mem targe big 隨機 max prime 神奇

題鏈:

http://www.lydsy.com/JudgeOnline/problem.php?id=1416
題解:

Pòlya瓦罐模型:
給定罐子裏每種顏色的球的個數A[i],按題目要求隨機操作若幹次之後,摸到i號球的概率仍然為A[i]/sum(A[i])。

所以可以把題目看成若幹次這種遊戲,得出每次取球時的概率,然後更新對應的A[i]和sum(A[i])。
把每次的概率相乘即是答案。

要用到高精度,為了避免除法,可以先把分子分母質因數分解(最大質因子不超過20000),
把分子分母相同的因子相互抵消後,再用高精度乘起來。


代碼:

#include<bits/stdc++.h>
#define BIT 10000
#define MAXN 1005
using namespace std;
struct Bigint{
	int A[2*MAXN],len;
	Bigint(){memset(A,0,sizeof(A)); len=1;}
	void operator = (int rtm){
		*this=Bigint();
		do{
			A[len++]=rtm%BIT; rtm/=BIT;
		}while(rtm);
		len--;
	}
	Bigint operator * (const Bigint &rtm) const{
		static Bigint now; now=Bigint();
		for(int i=1;i<=len;i++)
			for(int j=1;j<=rtm.len;j++){
				now.A[i+j-1]+=A[i]*rtm.A[j];
				now.A[i+j]+=now.A[i+j-1]/BIT;
				now.A[i+j-1]%=BIT;
			}
		now.len=len+rtm.len;
		while(!now.A[now.len]) now.len--;
		return now;
	}
	Bigint operator ^ (int b){
		static Bigint now,base;
		base=*this; now=1;
		while(b){
			if(b&1) now=now*base;
			b>>=1; base=base*base;
		}
		return now;
	}
	void Print(){
		printf("%d",A[len]);
		for(int i=len-1;i;i--)
			printf("%04d",A[i]);
	}
};
int T,N,D,S,pnt;
int A[MAXN],prime[20005],cnt[2][20005];
void prework(){
	static bool np[20005];
	for(int i=2;i<=20000;i++) if(!np[i]){
		prime[++pnt]=i; 
		for(int j=i*i;j<=20;j+=i) np[j]=1;
	}
}
void Decomposition(int k,int val){
	static int mincnt;
	if(val&&cnt[k][0]!=0) cnt[k][0]=1;
	else return (void)(cnt[k][0]=0);
	for(int i=1;i<=pnt&&val!=1;i++){
		while(val%prime[i]==0)
			cnt[k][i]++,val/=prime[i];
		mincnt=min(cnt[k][i],cnt[k^1][i]);
		cnt[k][i]-=mincnt; cnt[k^1][i]-=mincnt;
	}
}
Bigint combination(int k){
	Bigint ret,tmp; ret=cnt[k][0];
	for(int i=0;i<=pnt;i++) if(cnt[k][i]){
		tmp=prime[i];
		tmp=tmp^cnt[k][i];
		ret=ret*tmp;
	}
	return ret;
}
int main(){
	prework(); prime[0]=1;
	cnt[0][0]=cnt[1][0]=-1;
	scanf("%d%d%d",&T,&N,&D);
	for(int i=1;i<=T;i++)
		scanf("%d",&A[i]),S+=A[i];
	for(int i=1,x,y;i<=N;i++){
		scanf("%d%d",&x,&y);
		Decomposition(0,A[y]);
		Decomposition(1,S);
		A[y]+=D; S+=D;
	}
	Bigint a=combination(0);
	Bigint b=combination(1);
	if(!cnt[0][0]) b=1;
	a.Print(); putchar(‘/‘);
	b.Print(); putchar(‘\n‘);
	return 0;
}

  

●BZOJ 1416 [NOI2006]神奇的口袋