ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

「AGC020F」 Arcs on a Circle

2021-02-23 09:02:44  阅读:270  来源: 互联网

标签:Arcs int 部分 --- ge 圆弧 Circle AGC020F 小数


「AGC020F」 Arcs on a Circle

Link

这个题非常 Amazing 啊。果然AtCoder全是智商题

首先你可以注意到数据范围真的是小得离谱,让你想要爆搜。

然后你发现不可做,那考虑状压。

首先你发现这是一个环很烦,所以我们随便找一个端点断环为链。

问题转换为求能覆盖整个圆的不同的覆盖方式。

但是显然这样的方案有无数种,我们需要考虑优化。

注意到两段圆弧是否相交,仅与它们的起始位置的整数部分的值和小数部分有关。更进一步地,他们小数部分的相对大小决定了其是否相交。

若有 \(p_i\le p_j\),则两圆弧相交 \(\iff p_i+l_i\ge p_j\)。

拆分为整数和小数部分,有 \([p_i]+\{p_i\}+l_i\ge [p_j]+\{p_j\}\)。

若 \([p_i]+l_i\neq [p_j]\),则小数部分无影响。

若 \([p_i]+l_i= [p_j]\),则 \([p_i]+\{p_i\}+l_i\ge [p_j]+\{p_j\}\iff \{p_i\}\ge \{p_j\}\)。

不妨设相对位置互不相同。

于是我们可以考虑枚举所有圆弧小数部分的相对位置,然后状压DP即可。

/*---Author:HenryHuang---*/
/*---Never Settle---*/
#include<bits/stdc++.h>
using namespace std;
int s;
int a[10],tot;
long long ans;
int n,c;
int f[50][500];
int solve(){
	memset(f,0,sizeof f); f[0][a[n]*(n+1)]=1;
	for(int i=1;i<(n+1)*c;++i){
		if(i%(n+1)){
			int t=i%(n+1)-1,len=min(i+(n+1)*a[t],(n+1)*c);
			for(int j=0;j<=s;++j){
				if((j&(1<<t))==0){
					for(int k=i;k<=(n+1)*c;++k)
						f[j|(1<<t)][max(k,len)]+=f[j][k];
				}
			}
		}
	}
	return f[s][(n+1)*c];
}
int main(){
	ios::sync_with_stdio(0);
	cin.tie(0),cout.tie(0);
	cin>>n>>c;
	s=(1<<(n-1))-1;--n;
	for(int i=0;i<=n;++i) cin>>a[i];
	sort(a,a+n+1);
	do{
		ans+=solve();
		++tot;	
	}while(next_permutation(a,a+n));
	cout<<setprecision(14)<<fixed<<(long double)ans/tot/pow(c,n)<<'\n';
	return 0;
}

标签:Arcs,int,部分,---,ge,圆弧,Circle,AGC020F,小数
来源: https://www.cnblogs.com/HenryHuang-Never-Settle/p/14434040.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有