• sss


    <更新提示>


    <正文>

    T1 Gift

    Description

    ​ 人生赢家老王在网上认识了一个妹纸,然后妹纸的生日到了,为了表示自己的心 意,他决定送她礼物。可是她喜爱的东西特别多,然而他的钱数有限,因此他想 知道当他花一定钱数后剩余钱数无法再购买任何一件剩余物品(每种物品他最多 买一个)时有多少种方案,两种方案不同,当且仅当两种方案中至少有一件品不 同,可是由于他忙着准备泡下一个妹纸(chi),因此麻烦聪明的你帮帮忙。

    Input Format

    ​ 输入第一行 n 和 m, n 表示妹纸喜欢的礼物数目, m 表示现有的钱数,第二行 n 个数,表示 n 个物品的价格。

    Output Format

    ​ 输出一行一个数表示方案数目,答案对 1000000007 取模。

    Sample Input

    6 25 
    8 9 8 7 16 5
    

    Sample Output

    15
    

    Hint

    30%的数据: 0<=n<=100 0<=m<=500
    100%的数据:0<=n<=1000 0<=m<=1000
    注意:所有物品价格均小于 m

    解析

    如果存在一种合法的方案,就必然会有一个未购买并且价格最低的物品,并且这种方案的情况下,已经没有足够的钱来购买这个物品了。那么我们就以这个未购买的最小价格物品为基准点来统计答案,保证不重不漏。

    我们先将所有物品按照价格从小到大排序,每次枚举一个物品(i),令这个物品作为上述未购买且价格最低的物品,然后,我们强制取第(1)(i-1)个物品(如果不取,就和(i)个物品最小矛盾),并对第(i+1)(n)个物品用(0/1)背包统计方案数,那么(sum_{jin [m-a[i]+1,m]}f[j])即为本次的答案。

    我们发现每一次(dp)都是对后面连续的若干个物品进行计算,于是我们就可以放弃背包中滚动数组的做法,从(n)开始倒序(dp),并记录每一个阶段的(dp)值,然后在上述统计答案的过程中直接使用预处理的(dp)值即可。这样,时间复杂度就从(O(n^3))优化到了(O(n^2))

    (Code:)

    #include <bits/stdc++.h>
    using namespace std;
    const int N = 1020 , M = 1020;
    const long long Mod = 1000000007;
    int n,m,a[N];
    long long f[N][M],ans,sum[N];
    inline void input(void)
    {
    	scanf("%d%d",&n,&m);
    	for (int i=1;i<=n;i++)
    		scanf("%d",&a[i]);
    } 
    inline void dp(void)
    {
    	memset( f , 0 , sizeof f );
    	f[n+1][0] = 1;
    	for (int i=n;i>=1;i--)
    		for (int j=0;j<=m;j++)
    			if ( j >= a[i] ) f[i][j] = ( f[i+1][j] + f[i+1][j-a[i]] ) % Mod;
    			else f[i][j] = f[i+1][j];
    	for (int k=1;k<=n;k++)
    	{
    		for (int i=m;i>=sum[k-1];i--)
    			f[k+1][i] = f[k+1][i-sum[k-1]];
    		for (int i=0;i<min(sum[k-1],1LL*m+1);i++)
    			f[k+1][i] = 0;
    		for (int i=m;i>=m-a[k]+1;i--)
    			ans = ( ans + f[k+1][i] ) % Mod;
    	}
    }
    int main(void)
    {
    	input();
    	sort( a+1 , a+n+1 );
    	for (int i=1;i<=n;i++)
    		sum[i] = sum[i-1] + a[i];
    	dp();
    	if ( sum[n] > m ) printf("%lld
    ",ans%Mod);
    	else puts("1");
    	return 0;
    }
    

    T2 Fseq

    Description

    ​ 一个长度为 N+M 的数列,里面有 N 个+1,M 个-1 如果一个这样的数列被称作 F 序列(Fadeness) , 当且仅当它的任意前缀和均非 负。
    for example :
    1,-1,1,1,-1 is a Fadeness
    1,-1,-1,1,1 is not because S(3) <0
    求一个数列是 Fadensee 的概率。

    Input Format

    ​ 第一行, Test , 表示测试数据的组数。 每个数据 有两个数 N,M

    Output Format

    ​ 对于每组数据,输出一个实数(保留到小数点后 6 位)

    Sample Input

    3 
    1 0 
    0 1
    1 1
    

    Sample Output

    1.000000 
    0.000000 
    0.500000
    

    Hint

    30%的数据: (Test<=10),(0<=N,M<=1000).
    100%的数据: ( Test<=9008 ), ( 0<=N,M<=20000 ).

    解析

    和选举定理相似,我们直接转换成图形问题,然后应用反射原理即可。

    我们将数列取数形象地转化为坐标系中的移动问题,一开始在原点,每取一个(1)就向右上方走一步,横纵坐标各加(1),取一个(-1)就向右下方走一步,横坐标加(1),纵坐标减(1),那么最终折线会走到点((n+m,n-m))

    那么一旦折线接触到了直线(y=-1),这个序列就是不合法的。容易得出折线共有(C_{n+m}^m)种,那么问题在于求出不合法的折线数。

    应用反射原理,由点((0,-2))为起点,向右上走(n+1)次,向右下走(m-1)次,最终到点((n+m,n-m))的所有折线就是不合法的折线数,我们同样可以得出这样的折线有(C_{n+m}^{m-1})种。

    那么答案就是(frac{C_{n+m}^m-C_{n+m}^{m-1}}{C_{n+m}^m}),化简一下就是(1-frac{m}{n+1})

    (Code:)

    #include<bits/stdc++.h>
    using namespace std;
    const int T = 10020 , N = 20020;
    int t,n,m;
    int main(void)
    {
    	scanf("%d",&t);
    	while ( t-- )
    	{
    		scanf("%d%d",&n,&m);
    		if ( m > n ) puts("0.000000");
    		else printf("%.6lf
    ", 1.0 - 1.0 * m / (n+1) );
    	}
    	return 0;
    }
    

    <后记>

  • 相关阅读:
    集群任务管理系统SGE的简明教程
    三代全长转录组测序和组装
    quota
    基因组转座元件
    单倍体
    什么是寒武纪、志留纪、泥盆纪、白垩纪
    对组装之后的三代基因组进行polish
    使用 PhyML 构建进化树
    PAML软件中的mcmctree命令估算物种分歧时间
    OrthoMCL及orthofinder 两种软件进行聚类分析
  • 原文地址:https://www.cnblogs.com/Parsnip/p/11140899.html
Copyright © 2020-2023  润新知