bzoj2817[ZJOI2012]波浪
题目链接:
http://www.lydsy.com/JudgeOnline/problem.php?id=2817
波浪
【问题描述】
阿米巴和小强是好朋友。 阿米巴和小强在大海旁边看海水的波涛。小强第一次面对如此汹涌的海潮, 他兴奋地叫个不停。而阿米巴则很淡定,他回想起曾经的那些日子,事业的起伏, 情感的挫折……总之今天的风浪和曾经经历的那些风雨比起来,简直什么都不 算。 于是,这对好朋友不可避免地产生了分歧。为了论证自己的观点,小强建立 了一个模型。他海面抽象成一个 1 到N的排列P[1…N]。定义波动强度等于 相邻 两项的差的绝对值的和 ,即: L = | P2 – P1 | + | P3 – P2 | + … + | PN – PN-1 | 给你一个N和M,问:随机一个 1…N的排列,它的波动强度 不小于 M的概率 有多大? 答案请保留 小数点后 K 位输出,四舍五入 。
【输入格式】
输入文件 wavel.in 的第一行包含三个整数 N, M 和 K,分别表示排列的长度, 波动强度,输出位数。
【输出格式】
输出文件 wavel.out 包含一个小数点后 K 位的实数。
【样例输入】
3 3 3
【样例输出】
0.667
【样例说明】
N = 3 的排列有 6 个:123,132,213,231,312,321;他们的波动强度分 别为 2,3,3,3,3,2。所以,波动强度不小于 3 的概率是 4/6,即 0.667。 你也可以通过下面的代码来验证这个概率: int a[3]={0,1,2},s=0,n=3; for (int i=0;i<1000000;i++){ random_shuffle(a,a+n); int t=0; for (int j=0;j<n-1;j++) t+=abs(a[j+1]-a[j]);
if (t>=3) s++; } printf("%.3f\n",s/1000000.0);
【数据规模】
对于 30%的数据,N ≤ 10。 对于另外 30%的数据,K ≤ 3。 对于另外 30%的数据,K ≤ 8。 对于另外 10%的数据,N ≤ 50。 对于 100%的数据,N ≤ 100,K ≤ 30,0 ≤ M ≤ 2147483647。
题解
啊。。
我从未写过如此**的题目啊T^T
这种题估计也就fhq出的出来。。
从此我对小强和阿米巴产生了心理阴影。。
网上好像没多少题解啊。。那我再发个良心题解吧。。
好的现在开始正文。。
首先考虑从小到大在序列中插入每个数字。
比如。。
涂黑的格子表示已经放了数字的位置,然后。。我们现在要再放一个数字。
显然易见。。涂黑的格子的数字都比当前数字小
然后有这几种情况
红色格子表示当前填的数字
先看第一张图,新添的格子产生了新的一段!
而且,这一段两端的数还没填。。说明,当前位置贡献的权值要-2*i(假设i为当前填的数字),因为两端的数都比这个大,又因为你有(k+1)段可以插入新的一段,所以转移系数为(k+1)。。
再看第二张图,i合并了原来的两段,说明i两端的数字都填好了,这代表着当前位置贡献的权值为2*i,因为两端的数字都比这个小,又因为你有(k-1)段可以合并两个段,多已转移系数为(k-1)。。
发现了没?我们可以DP了!
设f[i][j][k][t]表示当前插入到第i个数字,权值和为j,有k段,t。。t表示整个序列两端的情况。比如t=0表示这个序列两端都没有放数字,1表示放了一个,2表示放了两个。
然后,我们可以对第一张图列出方程:f[i][j-2*i][k+1][0]+=f[i-1][j][k][0]*(k+1)
同理,可以对第二段建立转移方程:f[i][j+2*i][k-1][0]+=f[i-1][j][k][0]*(k-1)
第三张图就是f[i][j][k][0]+=f[i-1][j][k][0]*2*k
然后第四张图,出现了边界的情况,也依然可以DP。
1、当前放在边界上的数字和别的段合成了新的一段,那么i的贡献就是i,由于有两个端点,所以转移系数是2,即方程为:f[i][j+i][k][1]+=f[i-1][j][k][0]*2
2、当前放在边界上的数字创造了新的一段,那么i的贡献就是-i,由于有两个端点,所以转移系数是2,即方程为:f[i][j-i][k+1][1]+=f[i-1][j][k][0]*2
-----------------------------------------------------------------------------以上是t=0的情况。。t=1的情况也类似,只不过。。由于你能放的位置会减少。。所以转移系数也会改变。。具体看代码。----------------------------------------------------
对了,当t=2时,就不用考虑边界情况了!
好的大体内容讲完了!
但是
但是
但是
精度30位会炸啊!!!!!
这时神器来了:__float128
优点:精度高
缺点:速度慢速度慢速度慢速度慢速度慢速度慢速度慢速度慢速度慢空间大空间大空间大空间大空间大空间大空间大空间大
又由于出题人****卡内存卡时限。。
你会发现这个代码T了。。
#include <bits/stdc++.h>
#define MAX 4500
#define N 110
using namespace std;
typedef long long ll;
double f[][MAX*+][N][];
int i,j,k,n,m,x,y,t,q,la,no;
void print2(__float128 x,int k){
for (i=;i<=n;i++)x=x/(__float128)i;
ll a=x;x-=a;int s[],l,j=;s[]=;
while (a)s[++s[]]=a%,a/=;if (!s[])s[++s[]]=;l=s[];
for (i=;i<=k+;i++){x*=;a=(int)x;x-=a;s[++s[]]=a;}
if (s[s[]]>=)j=;int p=s[];
while (j&&p){s[p-]++;j=s[p-]/;s[p-]%=;p--;}
if (j){for (i=s[];i>;i--)s[i]=s[i-];s[]=;}
putchar(''+s[]);if (k)putchar('.');for (i=;i<=k+;i++)putchar(s[i]+'');
return;
}
inline void solve(int n,int m,int q){
f[][MAX-*][][]=;f[][MAX-][][]=;f[][MAX][][]=;no=;la=;
for (i=;i<=n;i++){
no=-no,la=-la;
memset(f[no],,sizeof f[no]);
for (j=;j<=MAX*;j++)
for (k=;k<=n-;k++){
if (f[la][j][k][]){
if (j>=*i)f[no][j-*i][k+][]+=f[la][j][k][]*(k+);
if (j+*i<=MAX*)f[no][j+*i][k-][]+=f[la][j][k][]*(k-);
f[no][j][k][]+=f[la][j][k][]*k*;
if (j+i<=MAX*)f[no][j+i][k][]+=f[la][j][k][]*;
if (j>=i)f[no][j-i][k+][]+=f[la][j][k][]*;
}
if (f[la][j][k][]){
if (j>=*i)f[no][j-*i][k+][]+=f[la][j][k][]*k;
if (j+*i<=MAX*)f[no][j+*i][k-][]+=f[la][j][k][]*(k-);
f[no][j][k][]+=f[la][j][k][]*(*k-);
if (j+i<=MAX*)f[no][j+i][k][]+=f[la][j][k][];
if (j>=i)f[no][j-i][k+][]+=f[la][j][k][];
}
if (f[la][j][k][]){
if (j+i*<=MAX*)f[no][j+i*][k-][]+=f[la][j][k][]*(k-);
if (j>=*i)f[no][j-i*][k+][]+=f[la][j][k][]*(k-);
f[no][j][k][]+=*f[la][j][k][]*(k-);
}
}
}
__float128 ans=;
for (i=m+MAX;i<=MAX*;i++)ans+=(__float128)(f[no][i][][]*1.000000000000);
print2(ans,q);
}
int main(){
scanf("%d%d%d",&n,&m,&q);
solve(n,m,q);
return ;
}
TLE啊。。
555~
本机测试22s。。其中两个点9s。。
说明
bzoj太慢了!!!
bzoj太慢了!!!
bzoj太慢了!!!
bzoj太慢了!!!
bzoj太慢了!!!
bzoj太慢了!!!
bzoj太慢了!!!
bzoj太慢了!!!
这时,膜了一下别人的代码,发现好妙啊
原来可以将k<=8的用double做,k>8的用__float128做!
写个namespace。。。inline什么的。。就过了!!跑得飞快啊!
AC代码
#include <bits/stdc++.h>
#define N 110
#define MAX 4500
using namespace std;
typedef long long ll;
int i,j,k,n,m,x,y,t,q,la,no;
namespace dob {typedef double db;db f[][][][];}
namespace fl {typedef __float128 db;db f[][][][];}
template <class T> inline
void print2(T x,int k){
for (i=;i<=n;i++)x=x/(T)i;
ll a=x;x-=a;int s[],l,j=;s[]=;
while (a)s[++s[]]=a%,a/=;if (!s[])s[++s[]]=;l=s[];
for (i=;i<=k+;i++){x*=;a=(int)x;x-=a;s[++s[]]=a;}
if (s[s[]]>=)j=;int p=s[];
while (j&&p){s[p-]++;j=s[p-]/;s[p-]%=;p--;}
if (j){for (i=s[];i>;i--)s[i]=s[i-];s[]=;}
putchar(''+s[]);if (k)putchar('.');for (i=;i<=k+;i++)putchar(s[i]+'');
return;
}
template <class T> inline
void solve(T f[][][][]){
f[][MAX-*][][]=;f[][MAX-][][]=;f[][MAX][][]=;no=;la=;
for (int i=;i<=n;i++){
no=-no,la=-la;
memset(f[no],,sizeof f[no]);
for (int j=;j<=MAX*;j++)
for (int k=;k<=n-;k++){
if (f[la][j][k][]){
if (j>=*i)f[no][j-*i][k+][]+=f[la][j][k][]*(k+);
if (j+*i<=MAX*)f[no][j+*i][k-][]+=f[la][j][k][]*(k-);
f[no][j][k][]+=f[la][j][k][]*k*;
if (j+i<=MAX*)f[no][j+i][k][]+=f[la][j][k][]*;
if (j>=i)f[no][j-i][k+][]+=f[la][j][k][]*;
}
if (f[la][j][k][]){
if (j>=*i)f[no][j-*i][k+][]+=f[la][j][k][]*k;
if (j+*i<=MAX*)f[no][j+*i][k-][]+=f[la][j][k][]*(k-);
f[no][j][k][]+=f[la][j][k][]*(*k-);
if (j+i<=MAX*)f[no][j+i][k][]+=f[la][j][k][];
if (j>=i)f[no][j-i][k+][]+=f[la][j][k][];
}
if (f[la][j][k][]){
if (j+i*<=MAX*)f[no][j+i*][k-][]+=f[la][j][k][]*(k-);
if (j>=*i)f[no][j-i*][k+][]+=f[la][j][k][]*(k-);
f[no][j][k][]+=*f[la][j][k][]*(k-);
}
}
}
T ans=;
for (int i=m+MAX;i<=MAX*;i++)ans+=(T)f[no][i][][];
print2(ans,q);
}
int main(){
scanf("%d%d%d",&n,&m,&q);if (q<=)solve(dob::f);else solve(fl::f);
return ;
}
AC啦~
本机测试5s!!
飞起~
然后就是几个注意点。。
1:用滚动数组。。
2:要输出优化。。我曾分类讨论了K为1~30的所有情况。。。
3、做好卡评测的准备。。
4、为bzoj买一台好点的机子
完结?撒花!
bzoj2817[ZJOI2012]波浪的更多相关文章
- 洛谷2612&&bzoj2817 [ZJOI2012]波浪
洛谷2612&&bzoj2817 [ZJOI2012]波浪 原题链接 题解 因为有abs不太好搞,考虑拆掉abs. 生成排列的方法之一:n个空位,从1到n一次插入一个空位. 这样搞的话 ...
- 【BZOJ2817】[ZJOI2012]波浪(动态规划)
[BZOJ2817][ZJOI2012]波浪(动态规划) 题面 BZOJ 洛谷 题解 首先这个差值最大也就是\(n^2\)级别的. 那么这样子就可以压进状态啦. 我们把这个操作看成一个个加数的操作,按 ...
- [ZJOI2012]波浪
Description: L = | P2 – P1 | + | P3 – P2 | + - + | PN – PN-1 | 给你一个N和M,问:随机一个1-N的排列,它的波动强度(L)不小于M的概率 ...
- Luogu2612 ZJOI2012 波浪 DP
传送门 花掉了自己用来搞学科的时间做了这道题-- 一道类似的题:Here 考虑拆开绝对值计算贡献.那么我们对于\(1\)到\(N\)的排列,从小到大地将插入它们插入排列中. 假设我们现在计算到了数\( ...
- [ZJOI2012]波浪弱化版(带技巧的DP)
题面 \(solution:\) 这道确实挺难的,情况特别多,而且考场上都没想到如何设置状态.感觉怎么设状态不能很好的表示当前情况并转移,考后发现是对全排列的构造方式不熟而导致的,而这一题的状态也是根 ...
- Luogu P2612 [ZJOI2012]波浪
题目 我们考虑从\(1\)到\(n\)把每个数放到序列里面去,以消掉绝对值. 在最后的序列中,如果\(i\)的某一边是序列的边界,那么\(i\)会产生\(0\)的贡献.如果\(i\)的某一边是一个比\ ...
- 题解 洛谷 P2612 【[ZJOI2012]波浪】DP+高精
题目描述 题目传送门 分析 因为有绝对值不好处理,所以我们强制从小到大填数 设 \(f[i][j][p][o]\) 为当前填到了第 \(i\) 个数,波动强度为 \(j\),有 \(p\) 个连续段并 ...
- Noip前的大抱佛脚----赛前任务
赛前任务 tags:任务清单 前言 现在xzy太弱了,而且他最近越来越弱了,天天被爆踩,天天被爆踩 题单不会在作业部落发布,所以可(yi)能(ding)会不及时更新 省选前的练习莫名其妙地成为了Noi ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
随机推荐
- 20155209 林虹宇 Exp3 免杀原理与实践
Exp3 免杀原理与实践 使用msf生成后门程序的检测 将上周msf生成的后门文件放在virscan.org中进行扫描 结果很危险 使用msf编码一次进行扫描 使用msf编码10次进行扫描 结果同样很 ...
- 《网络对抗》Exp5 MSF基础应用
20155336<网络对抗>Exp5 MSF基础应用 一.基础知识回答 用自己的话解释什么是exploit,payload,encode exploit:渗透攻击的模块合集,将真正要负责攻 ...
- Android开发——Android进程保活招式大全
)前台进程(Foreground process),即用户当前操作所必需的进程,通常数量不多.举例如下: //拥有用户正在交互的 Activity(已调用 onResume()) //拥有某个 Ser ...
- Luogu P1198 [JSOI2008]最大数
我会用高级(???)的单调栈来打这道题吗? 线段树即可水过. 假设这个数列刚开始所有数都是0,然后我们每次只要进行一个点的修改和区间求和即可. 这不就是 线段树大法. 只要用一个len记录一下当前数列 ...
- Remote 桌面的win2003 servre端设定
Microsoft Windows [版本 5.2.3790](C) 版权所有 1985-2003 Microsoft Corp. C:\Documents and Settings\Administ ...
- Linux中tty、pty、pts的概念区别 转载
基本概念: > tty(终端设备的统称): tty一词源于Teletypes,或teletypewriters,原来指的是电传打字机,是通过串行线用打印机键盘通过阅读和发送信息的东西,后来这东西 ...
- Linux 学习日记 1
这是我第一次系统地学习Linux,希望通过这个学习日记收获一些东西把-- @_@ Grub - 启动管理器 在启动时让用户选择要启动的系统.(但是windows比较霸道--重装windows后会将 ...
- Markdown打造高逼格博客
这里首先假设读者你已经掌握了Markdown与GitHub的基本用法 如果不会, 请先自行百度或Google, 我目前还没写Markdown与GitHub的教程 看云只是一个推荐, 可以认为协助生成格 ...
- Asp.Net_Mvc3.5语法_<%%>的用法
一. <%%>这种格式实际上就是和asp的用法一样的,只是asp中里面是vbscript或 者javascript代码,而在asp.net中用的是.net平台下支持的语言.特别 注意:服务 ...
- NodeMCU学习(一) : 开始之前的准备
安装Aduino开发环境 在官网中下载Arduino开发环境,或者在网盘中下载: 网盘地址: https://pan.baidu.com/s/1OjMhYgKOYW69YC2dEwFgyw: 提取码: ...