POJ 2443 Set Operation -- 位运算 + bitset

/*
	http://poj.org/problem?id=2443 
	POJ2443 Set Operation -- 位运算 + bitset
	询问两个元素在不在同一个集合内
	(在 ? yes : no)
*/

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
#define CLR(c,v) memset(c,v,sizeof(c));
int const N = 1e5 +10;

double loge[N];
double nCr(int n,int r){
	if(n==0)return 0.;
	return loge[n] - loge[r] - loge[n-r];
}
int a[N],sum[N];

int main(){
	//freopen("in.txt","r",stdin);
	loge[0] = 0;
	for(int i = 1 ; i < N ; i++)
		loge[i] += loge[i-1] + log(i*1.0);
	
	int n,k;
	while(cin >> n >> k,n||k){
		CLR(sum,0);
		for(int i = 0 ; i < n ; i++){
			scanf("%d",&a[i]);
			sum[i+1] += sum[i] + a[i];
		}
		double ans = 0;
		for(int i = k ; i <= n ; i++){ // seq length
			double p = exp( nCr(i-2,k-2) - nCr(n,k) );
			ans += p * ( (sum[n]-sum[i-1])-(sum[n-i+1]-sum[0]) );
		}
		printf("%.3f\n" , ans);
	}
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
POJ - 3616是一个题目,题目描述如下: 给定一组区间,每个区间有一个权重,要求选择一些区间,使得这些区间的右端点都小于等于k,并且权重之和最大。请问最大的权重和是多少? 解决这个问题的思路是使用动态规划。首先,将区间按照左端点从小到大进行排序。然后,定义一个dp数组,dp[i]表示右端点小于等于i的所有区间所能得到的最大权重。 接下来,遍历每一个区间,对于每个区间i,将dp[i]初始化为区间i的权重。然后,再遍历i之前的每个区间j,如果区间j的右端点小于等于k,并且区间j的权重加上区间i的权重大于dp[i],则更新dp[i]为dp[j]加上区间i的权重。 最后,遍历整个dp数组,找到最大的权重和,即为所求的答案。 下面是具体的代码实现: ```cpp #include <cstdio> #include <cstring> #include <algorithm> using namespace std; struct interval{ int start, end, weight; }; interval intervals[10005]; int dp[10005]; int n, m, k; bool compare(interval a, interval b) { if (a.start == b.start) { return a.end < b.end; } else { return a.start < b.start; } } int main() { while(~scanf("%d %d %d", &n, &m, &k)) { memset(dp, 0, sizeof dp); for (int i = 0; i < m; i++) { scanf("%d %d %d", &intervals[i].start, &intervals[i].end, &intervals[i].weight); } sort(intervals, intervals + m, compare); for (int i = 0; i < m; i++) { dp[i] = intervals[i].weight; for (int j = 0; j < i; j++) { if (intervals[j].end <= k && dp[j] + intervals[i].weight > dp[i]) { dp[i] = dp[j] + intervals[i].weight; } } } int maxWeight = 0; for (int i = 0; i < m; i++) { maxWeight = max(maxWeight, dp[i]); } printf("%d\n", maxWeight); } } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值