●BZOJ 1416 [NOI2006]神奇的口袋
题链:
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]神奇的口袋的更多相关文章
- BZOJ 1416: [NOI2006]神奇的口袋( 高精度 )
把x1~xn当成是1~n, 答案是不会变的. 然后直接模拟就行了...... P.S 双倍经验... BZOJ1416 && BZOJ1498 -------------------- ...
- Bzoj1498&1416: [NOI2006]神奇的口袋
什么鬼双倍经验题??? Sol 考虑在第\(k\)次摸到\(y\)的概率 如果上次摸到\(y\),目前有\(sum\)个球,\(y\)有\(a[y]\)个,那么概率就是\(\frac{a[y]+d}{ ...
- 【BZOJ1416/1498】【NOI2006】神奇的口袋(数论,概率)
[BZOJ1416/1498][NOI2006]神奇的口袋(数论,概率) 题面 BZOJ1416 BZOJ1498 洛谷 题面都是图片形式是什么鬼.. 题解 考虑以下性质 1.\(x[1],x[2]. ...
- BZOJ 1006 [HNOI2008] 神奇的国度(简单弦图的染色)
题目大意 K 国是一个热衷三角形的国度,连人的交往也只喜欢三角原则.他们认为三角关系:即 AB 相互认识,BC 相互认识,CA 相互认识,是简洁高效的.为了巩固三角关系,K 国禁止四边关系,五边关系等 ...
- 神奇的口袋(dp)
有一个神奇的口袋,总的容积是40,用这个口袋可以变出一 些物品,这些物品的总体积必须是40. John现在有n(1≤n ≤ 20)个想要得到的物品,每个物品 的体积分别是a1,a2……an.John可 ...
- [codeup] 2044 神奇的口袋
题目描述 有一个神奇的口袋,总的容积是40,用这个口袋可以变出一些物品,这些物品的总体积必须是40.John现在有n个想要得到的物品,每个物品的体积分别是a1,a2--an.John可以从这些物品中选 ...
- dp 神奇的口袋
有一个神奇的口袋,总的容积是40,用这个口袋可以变出一 些物品,这些物品的总体积必须是40. John现在有n(1≤n ≤ 20)个想要得到的物品,每个物品 的体积分别是a1,a2--an.Joh ...
- 九度OJ 1114:神奇的口袋 (DFS、DP)
时间限制:1 秒 内存限制:32 兆 特殊判题:否 提交:948 解决:554 题目描述: 有一个神奇的口袋,总的容积是40,用这个口袋可以变出一些物品,这些物品的总体积必须是40.John现在有n个 ...
- 九度oj 题目1114:神奇的口袋
题目描述: 有一个神奇的口袋,总的容积是40,用这个口袋可以变出一些物品,这些物品的总体积必须是40.John现在有n个想要得到的物品,每个物品的体积分别是a1,a2……an.John可以从这些物品中 ...
随机推荐
- Beta 第一天
一.今日任务 重新熟悉整体项目 对整个项目在未来的beta冲刺中进程有一个合理的规划 由于我们送出的是一个负责前端的成员,引入的也是一个负责前端工作的女生,(女生做起美工比起男生更加得心应手吧)所以我 ...
- 学号:201621123032 《Java程序设计》第14周学习总结
1:本周学习总结 2:使用数据库技术改造你的系统 2.1:简述如何使用数据库技术改造你的系统.要建立什么表?截图你的表设计. 建立一个图书馆的表 建立读者用户个人的借书信息表---但是目前没有办法做到 ...
- 简单的C语言编译器--语义制导翻译
语法分析是最难写的,而这部分确实最伤脑的.大量的语义动作分析差点把我逼疯. 简而言之,这部分的作用就是在每次归约之后,都进行一些语义动作,最终让我们得到测试程序的三地址码,即中间代码. 1. ...
- iOS极光推送SDK的使用流程
一.极光推送简介 极光推送是一个端到端的推送服务,使得服务器端消息能够及时地推送到终端用户手机上,整合了iOS.Android和WP平台的统一推送服务.使用起来方便简单,已于集成,解决了原生远程推送繁 ...
- video与audio的使用
HTML5 DOM 为 <audio> 和 <video> 元素提供了方法.属性和事件. 这些方法.属性和事件允许您使用 JavaScript 来操作 <audio> ...
- Oracle银行存取钱系统
Oracle银行存取钱系统 /* 银行系统 要求: 1.创建一个用户信息表(userinfo).一个交易信息表(deal) 2.用户信息表字段:用户编号.用户名称.密码.余额 交易信息表字段:编号.交 ...
- javascript 中的类型
javascript 中的类型 js 是一门弱语言,各式各样的错误多种多样,特别是确定返回值有问题的时候,你会用什么来进行表示错误? 我一般有三个选择: null '' error {} 第一个选择 ...
- BizTalk 2016 配置 RosettaNet遇到的坑
本文只针对已经安装好BizTalk 2016 需要在安装RosettaNet加速器的伙伴. IIS配置 权限问题 错误信息 Failed to get IIS metabase property. E ...
- 用‘+=’拼接字符串,打印时总会出现一个undefined
var str; for(var i = 0; i < 5; i++){ str += String(i); } console.log(str); 他喵的,打印的结果竟然是"unde ...
- Docker学习笔记 - Docker客户端和服务端
学习内容: Docker客户端和服务端的通讯方式:client和自定义程序 Docker客户端和服务端的连接方式:socket 演示Docker客户端和服务端之间用remote-api通讯:nc ...