Throwing Dice LightOJ - 1064 || (勉强能用的)分数类
方法:
设ans[i][j]表示i个骰子点数恰好为j的概率。那么ans[1][1]到ans[1][6]都为1/6。
显然,$ans[i][j]=sum\{ans[i-1][j-k]\}(1<=k<=6,j>k)$
n和x上限很小,直接处理出所有点数恰好为某个值的结果,然后再做一遍类似前缀和的东西处理出所有点数大于等于某个值的结果。这里答案需要分数,于是乱写了一个分数结构体。
分数运算(没有优化的):
初始化空的分数:分子为0,分母为1
加法:先通分,然后加分子,然后约分
乘法:先互相约分,然后分子分母分别相乘
错误次数:1次
错误原因:由于把前缀和的范围限制与原数据的范围限制搞混,66行>=0写成>=i,导致WA
#include<cstdio>
typedef long long LL;
LL gcd(LL a,LL b)
{
LL t;
while(b!=)
{
t=a;
a=b;
b=t%b;
}
return a;
}
struct X
{
LL a,b;
X()
{
a=;b=;
}
X(int x,int y)
{
a=x;b=y;
}
X operator+(const X& c)
{
X ans;
LL tmp=gcd(b,c.b);
ans.b=b/tmp*c.b;
ans.a=ans.b/b*a+ans.b/c.b*c.a;
tmp=gcd(ans.a,ans.b);
ans.a/=tmp;
ans.b/=tmp;
return ans;
}
X operator*(const X& c)
{
LL tmp1=gcd(a,c.b),tmp2=gcd(b,c.a);
X ans;
ans.a=a/tmp1*c.a/tmp2;
ans.b=b/tmp2*c.b/tmp1;
return ans;
}
void print()
{
if(b==)
printf("%lld",a);
else
printf("%lld/%lld",a,b);
}
}ans[][],ans2[][];
LL T,TT,n,x;
int main()
{
int i,j,k;
for(i=;i<=;i++)
ans[][i]=(X){,};
for(i=;i<=;i++)
for(j=i;j<=*i;j++)
for(k=;k<=;k++)
if(j-k>)
ans[i][j]=ans[i][j]+ans[i-][j-k]*ans[][];//ans[1][1]就是1/6
for(i=;i<=;i++)
{
ans2[i][*i]=ans[i][*i];
for(j=*i-;j>=;j--)
ans2[i][j]=ans2[i][j+]+ans[i][j];
}
scanf("%lld",&T);
for(TT=;TT<=T;TT++)
{
scanf("%lld%lld",&n,&x);
printf("Case %lld: ",TT);
ans2[n][x].print();
puts("");
}
return ;
}
Throwing Dice LightOJ - 1064 || (勉强能用的)分数类的更多相关文章
- Throwing Dice(概率dp)
C - Throwing Dice Time Limit:2000MS Memory Limit:32768KB 64bit IO Format:%lld & %llu Lig ...
- lightoj 1064 Throwing Dice
题意:给你n个骰子,求n个骰子的和不小于x的概率. 刚开始想每给一组数就计算一次~~太笨了- -,看了别人的代码,用dp,而且是一次就初始化完成,每次取对应的数据就行了.WA了好多次啊,首先不明白的就 ...
- Light OJ 1064 - Throwing Dice
题目大意: 给你n个骰子, 问点数大于等于x的概率是多少? #include<cstdio> #include<cstring> #include<iostream> ...
- LightOJ1064 Throwing Dice(DP)
第一眼以为是概率DP,我还不会.不过看题目那么短就读读,其实这应该还不是概率DP,只是个水水的DP.. dp[n][s]表示掷n次骰子点数和为s的情况数 dp[0][0]=1 dp[i][j]=∑dp ...
- lightoj刷题日记
提高自己的实力, 也为了证明, 开始板刷lightoj,每天题量>=1: 题目的类型会在这边说明,具体见分页博客: SUM=54; 1000 Greetings from LightOJ [简单 ...
- 10409 - Die Game
Problem G: Die Game Life is not easy. Sometimes it is beyond your control. Now, as contestants of AC ...
- UVA 657 The die is cast
The die is cast InterGames is a high-tech startup company that specializes in developing technolo ...
- The Die Is Cast(poj 1481简单的双dfs)
http://poj.org/problem?id=1481 The Die Is Cast Time Limit: 1000MS Memory Limit: 10000K Total Submi ...
- Boost 1.61.0 Library Documentation
http://www.boost.org/doc/libs/1_61_0/ Boost 1.61.0 Library Documentation Accumulators Framework for ...
随机推荐
- Xammp修改端口
How can I get XAMPP working on port 80 under Windows 10? By default, Windows 10 starts Microsoft IIS ...
- Arcgis Engine(ae)接口详解(7):地图(map)操作
IMap map = null; //跟map同一层次的activeView对象,他们都是“地图”的对象,map管理地图内容,activeView管理显示内容 IActiveView activeVi ...
- iOS优化策略
之前写的一系列文章或者小经验一直没有时间去整理放在博客上,今天整理出来,之前是写在作业部落,语法是markdown,点击链接浏览,仅供参考,希望对你有帮助. https://www.zybuluo.c ...
- python 目录下的__init__.py
1 一个目录要成为一个package必须有__init__.py文件 The __init__.py files are required to make Python treat the direc ...
- Flume-ng-sdk源码分析
Flume 实战(2)--Flume-ng-sdk源码分析 - mumuxinfei - 博客园 http://www.cnblogs.com/mumuxinfei/p/3823266.html
- Java,如何获取文件的MD5值
MessageDigest类封装得很不错,简单易用 不多说,直接上代码 import java.io.FileInputStream;import java.security.MessageDiges ...
- C# 自定义控件及引用自动义控件
1.http://www.cnblogs.com/hjxzjp/p/7823292.html 优先考虑从现有的控件中进行派生,并添加所需要的功能. 在解决方案资源管理器窗口中设置:引用----&g ...
- 关于JAVA中的前期绑定 后期绑定(动态绑定)
前期绑定,在程序执行前根据编译时类型绑定,调用开销较小,如C语言只有前期绑定这种方法调用. 后期绑定,是指在运行时根据对象的类型进行绑定,又叫动态绑定或运行时绑定.实现后期绑定,需要某种机制支持,以便 ...
- design.js
//模块式开发 var myNamespace = (function () { var myPrivateVar = 0; var myPrivateMethod = function (foo) ...
- java 连接飞信API
通过java连接飞信api给自己的好友(包括自己)发送飞信内容.如果对方的手机号非你的飞信好友则不能发送.1. [代码]飞信发送类 package per.artisan.fetion; ...