首页 > 其他分享 >TopCoder SRM478C RandomApple 题解

TopCoder SRM478C RandomApple 题解

时间:2024-02-17 20:22:06浏览次数:34  
标签:箱子 SRM478C int 题解 pb 苹果 siz TopCoder define

题意:有\(k\)种苹果和\(n\)个箱子,每个箱子中有一些苹果,先等概率选取\(n\)个箱子组成集合的非空子集,再从选出的苹果中随机选一个,问每种苹果被选中的概率是多少

箱子\(i\)有\(a_{i,j}\)个第\(j\)种苹果,第\(i\)个箱子的总苹果数\(siz_i=\sum\limits_{j=1}^ka_{i,j}\),苹果总数\(sum=\sum\limits_{i=1}^n siz_i\)

每个箱子中的苹果被选中的概率是相同的,所以先考虑选箱子,设\(f_{i,j}\)表示在前\(i\)个箱子中选\(j\)个苹果的方案数(只能一箱一箱选),那么\(f_{i,j}=f_{i-1,j}+f_{i-1,j-siz_i}\)

再设\(g_{i,j}\)表示在所有不是\(i\)的箱子中选\(j\)个苹果的方案数,那么\(g_{i,j}=f_{n,j}-g_{i,j-siz_i}\)

于是,选中了箱子\(i\)的情况对箱子\(i\)中每个苹果被选中的概率贡献为\(\sum\limits_{j=0}^{sum-siz_i}\frac{g_{i,j}}{2^n-1}\times\frac{1}{j+siz_i}\)

解释一下

\(f\) 的递推显然。

\(g\) 的递推看似很神奇,其实移项后可以发现 \(g_{i,j}+g_{i,j-siz_i}=f_{n,j}\),也就是 \(f_{n,j}\) 的递推,当不选第 \(i\) 个箱子时方案数为 \(g_{i,j}\),反之即为选第 \(i\) 个箱子里的 \(siz_i\) 个苹果,那么方案数为其他箱子选 \(j-siz_i\) 个苹果的方案数,也就是 \(g_{i,j}\)。

计算答案时,枚举到 \(i\),是已经假设箱子 \(i\) 被选中后的;枚举 \(j\) 表示其他箱子里的苹果数,从其他箱子里选出 \(j\) 个苹果的方案数为 \(g_{i,j}\),对于选出的任意子集,其概率为 \(\frac{g_{i,j}}{2^n-1}\),在选出来的所有苹果中,每一个苹果被选中的概率为 \(\frac{1}{j+siz_i}\)。

不知道为什么本题不压维过不了。

#include <bits/stdc++.h>
using namespace std;
#define srand srand(time(NULL))
#define random(x) rand() % (x)
#define il inline
#define ptc putchar
#define reg register
#define mp make_pair
#define pb push_back
#define R(i, l, r) for (int i = l; i <= r; ++i)
#define debug puts("--------------------------------------------")
typedef long long ll;
typedef pair<int, int> PII;
namespace kunkun
{
    template <typename T>
    il void read(T &x)
    {
       x = 0; T f = 1; char ch;
       while (!isdigit(ch = getchar())) f -= (ch == '-') << 1;
       while (isdigit(ch)) x = (x << 1) + (x << 3) + (ch & 15), ch = getchar(); x *= f;
    }
    template <typename T, typename ...L>
    il void read(T &x, L &...y) {read(x); read(y...);}
    template <typename T>
    il void write(T x)
    {
        if (x < 0) ptc('-'), x = -x;
        if (x > 9) write(x / 10);
        ptc(x % 10 + '0');
    }
    template <typename T, typename ...L>
    il void write(T &x, L &...y) {write(x), ptc(' '); write(y...);}
}
using namespace kunkun;
const int N = 55;
// O(50*50*50*200)
int n, k, a[N][N], siz[N], sum;
ll f[500005], g[500005];
class RandomApple
{
	public:
		vector <double> theProbability(vector <string> hundred, vector <string> ten, vector <string> one)
		{
			n = hundred.size();
			k = hundred[0].size();
			for (int i = 1; i <= n; ++i)
			{
				for (int j = 1; j <= k; ++j)
				{
					a[i][j] = (hundred[i - 1][j - 1] - '0') * 100 + (ten[i - 1][j - 1] - '0') * 10 + one[i - 1][j - 1] - '0';
					siz[i] += a[i][j];
				}
				sum += siz[i];
			}
//			cout << "sum: " << sum << endl;
			vector <double> ans;
			R(i, 1, k) ans.pb(0);
			// 设f[i][j]表示在前i个箱子中选j个苹果的方案数
			f[0] = 1;
//			f[1][siz[1]] = 1;
			for (int i = 1; i <= n; ++i)
			{
				for (int j = sum; j >= siz[i]; --j) f[j] += f[j - siz[i]];
//				for (int j = 0; j <= sum; ++j)
//				{
//					f[i][j] = f[i - 1][j];
//					if (j >= siz[i]) f[i][j] += f[i - 1][j - siz[i]];
//				}
			}
			for (int i = 1; i <= n; ++i)
			{
				long double tmp = 0;
				for (int j = 0; j <= sum; ++j)
				{
//					cout << i << ' ' << j << endl;
					g[j] = f[j];
					if (j >= siz[i]) g[j] -= g[j - siz[i]];
					if (j + siz[i] <= sum) tmp += g[j] * 1.0 / (j + siz[i]);
				}
				tmp /= ((long double)pow(2, n) - 1);
				for (int j = 1; j <= k; ++j)
				{
					ans[j - 1] += tmp * a[i][j];
					// 选到箱子i中任意一个苹果的概率为tmp 
				}
//				cout << "???\n";
			}
//			debug;
			return ans;
		}
} ;
/*
signed main()
{
	vector <string> h, t, o;
	h.pb("01010110"), h.pb("00011000"), h.pb("00001000"), h.pb("10001010"), h.pb("10111110");
	t.pb("22218214"), t.pb("32244284"), t.pb("68402430"), t.pb("18140323"), t.pb("29043145");
	o.pb("87688689"), o.pb("36101317"), o.pb("69474068"), o.pb("29337374"), o.pb("87255881");
//	h.pb("10");
//	t.pb("00");
//	o.pb("00");
//	h.pb("0000"), h.pb("0000"), h.pb("0000");
//	t.pb("2284"), t.pb("0966"), t.pb("9334");
//	o.pb("1090"), o.pb("3942"), o.pb("4336");
//	h.pb("00"), h.pb("00");
//	t.pb("00"), t.pb("00");
//	o.pb("21"), o.pb("11");
//	debug;
	vector <double> ans = solver.theProbability(h, t, o);
	for (auto x : ans) printf("%.10f ", x);
    return 0;
}
//*/

标签:箱子,SRM478C,int,题解,pb,苹果,siz,TopCoder,define
From: https://www.cnblogs.com/cyyhcyyh/p/18018312

相关文章

  • 无限酒店 题解
    题目链接由于间隔不变,对于下面\(3\)个操作,只需记录起始位置与间隔即可。对于无数人到达酒店:所有位置的起始点\(\times2\),间隔\(\times2\),新的团队起始点为\(1\),间隔为\(2\)。对于\(k\)个人到达酒店:所有点的起始点\(+k\),间隔不变,新的团队起始点为\(0\),间隔为\(1\)......
  • CF1365G Secure Password 题解
    Description本题是交互题。有一个固定的数组\(A\),同时通过数组\(A\)构造出数组\(P\),具体来讲,\(P_i\)是\(A\)中除\(A_i\)外的所有元素的按位或。你需要在最多\(13\)次询问中得到最后的\(P\)数组。\(2\leqn\leq1000\)。Solution首先有一个\(2\logn\)的是注......
  • 传纸条 题解
    题目描述小渊和小轩是好朋友也是同班同学,他们在一起总有谈不完的话题。一次素质拓展活动中,班上同学安排做成一个m行n列的矩阵,而小渊和小轩被安排在矩阵对角线的两端,因此,他们就无法直接交谈了。幸运的是,他们可以通过传纸条来进行交流。纸条要经由许多同学传到对方手里,小渊坐在矩阵......
  • 2024.2.17模拟赛T1题解
    先考虑\(q=(1...n)\)的情况:发现如果设\(divcnt(p)\)表示将\(p\)划分为极小值域连续段的个数,那满足\(divcnt(p)\gem\)的排列都是合法的。那现在要求出有多少个排列符合条件可以先算出长度为\(i\),\(divnct\)为\(1\)的排列个数,这可以用dp解决然后再背包一下,就能求......
  • P2042 [NOI2005] 维护数列 题解
    题目链接:维护数列比较不好码的题,首先确保自己会一种文艺平衡树的书写,这点因人而异,比较推荐初学者学\(fhq\)平衡树,坑比较少,比较好码,写平衡树合并、持久化类的题时,也比较好写。注意到空间需求比较大,还涉及删除,我们的删除用各种树类数据结构中最常用的回收标记用于新的节点开辟。......
  • CF1929E 题解
    题意:给定一棵\(n\)个节点数和\(k\)条路径\((a_i,b_i)\),求至少将多少条边染色,使得给定路径都至少有一条染色的边。\(n\le10^5,k\le20\)。思路:好题。显然状压\(dp\),\(dp[S]\)表示至少染多少条边使得\(S\)中的路径都满足条件。正常思路是枚举子集转移,考虑\(T\s......
  • 整数划分 题解
    题目描述如何把一个正整数N(N长度<20)划分为M(M>1)个部分,使这M个部分的乘积最大。N、M从键盘输入,输出最大值及一种划分方式。输入格式第一行一个正整数T(T<=10000),表示有T组数据。接下来T行每行两个正整数N,M。输出格式对于每组数据第一行输出最大值。第二行输出划分方案,将N按......
  • 情书密码 题解
    题目描述有消息称:绝恋找到了自己的MissRight,正准备自己的表白。绝恋已经写好了情书,但为了避免其它人截获,他对情书进行加密。为了打探绝恋的私密,你冒着生命危险终于搞到了这封情书。原以为可以轻易将情书解密,结果竟然发现聪明的绝恋并没有直接写出加密用的密码,而是在那粉红色的......
  • CF484B题解
    朴素的办法是枚举每两个数然后更新取模的结果。时间复杂度为\(O(n^2)\)不能通过。这个朴素的过程可以看作枚举一个数然后找对其取模最大值的过程。我们可以枚举一个数,然后再枚举以它的倍数为两端的区间,找其中取模的最大值。找最大值的过程可以二分或双指针。值域很小,也可以......
  • 牛客小白87题解A-G
    A小苯的石子游戏本题贪心模拟即可B小苯的排序疑惑考虑到最简单的操作把n-1个数排好,最后一个看能否有序。即:a[1]<=min(a[2],a[3],a[4]..,a[n])||a[n]>=max(a[1],a[2],a[3]....,a[n-1])满足条件,反之易得不可能C&D小苯的IDE括号问题本题考察题意理解,简单版本我们可以只关注逻......