题目描述

Given a list of N integers with absolute values no larger than 10 15, find a non empty subset of these numbers which minimizes the absolute value of the sum of its elements. In case there are multiple subsets, choose the one with fewer elements.

输出

The input contains multiple data sets, the first line of each data set contains N <= 35, the number of elements, the next line contains N numbers no larger than 10 15 in absolute value and separated by a single space. The input is terminated with N = 0

输入

For each data set in the input print two integers, the minimum absolute sum and the number of elements in the optimal subset.

样例

样例输入

1
10
3
20 100 -100
0

样例输出

10 1
0 2

分析

一句话题意:多组数据,以 0 为结尾,给你n个数,求出这n个数的一个非空子集,使子集中的数加和的绝对值最小,在此基础上子集中元素的个数应最小。

N的最大值为35,暴力枚举肯定是要超时的,因为你要把235种情况都枚举出来,即使30秒的时限也是不够用的

暴力枚举

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn=50;
ll a[maxn],ans;
int n,cnt;
ll mabs(ll x){
if(x>0) return x;
return -x;
}
void dfs(int now,ll tot,int xuan){
if(xuan && mabs(tot)<ans){
cnt=xuan;
ans=mabs(tot);
}
else if(xuan && ans==mabs(tot)) cnt=min(cnt,xuan);
for(int i=now;i<=n;i++){
dfs(i+1,tot+a[i],xuan+1);
}
}
int main(){
while(scanf("%d",&n)!=EOF && n!=0){
ans=0x3f3f3f3f3f3f3f3f;
cnt=40;
for(int i=1;i<=n;i++){
scanf("%lld",&a[i]);
}
dfs(1,0,0);
printf("%lld %d\n",ans,cnt);
}
return 0;
}

其实我们可以考虑把n个数分成两部分,每个部分最多有218种结果

我们先把前一半预处理处理出来,然后再枚举后一半

对于每一个值,我们都要从前一半中找到与它的相反数最接近的数,然后两数相加取最小值

需要注意的是,最优解也有可能出现在前一半或后一半,不一定是前一半中的数与后一半中的数相加的结果

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<map>
using namespace std;
typedef long long ll;
const int maxn=50;
int n,cnt=0;
map<ll,ll> ma;
//map记录当前值所选的最少元素个数
ll a[maxn],jl[1<<20],ans=1e18;
//jl数组存储前一半的结果,ans是最终结果
ll abss(ll x){
if(x>0) return x;
return -x;
}//绝对值函数
ll solve(ll num){
if(num<=jl[1]) return jl[1];
else if(num>=jl[cnt]) return jl[cnt];
else {
ll l=1,r=cnt,mids;
while(l<=r){
mids=(l+r)/2;
if(jl[mids]==num) break;
else if(jl[mids]<num) l=mids+1;
else r=mids-1;
}
if(jl[mids]==num){
return num;
}
else if(abss(jl[l]-num)==abss(jl[r]-num)){
if(ma[jl[l]]<ma[jl[r]]) return jl[l];
return jl[r];
}
else if(abss(jl[l]-num)>=abss(jl[r]-num)){
return jl[r];
} else {
return jl[l];
}
}
}//模板,查找与x最接近的元素,要注意的是如果有多个元素与x最接近,取数的个数最小的那一个
int main(){
while(scanf("%d",&n)!=EOF && n!=0){
cnt=0,ans=1e18;
ma.clear();
memset(a,0,sizeof(a));
memset(jl,0,sizeof(jl));
for(int i=1;i<=n;i++){
scanf("%lld",&a[i]);
}//初始化+输入
if(n==1){
printf("%lld 1\n",abss(a[1]));
continue;
}//如果只有一个数,特判一下,因为前一半没有数,无法二分
int le=n/2,re=n-le;
int mmax=(1<<le)-1;
ll id=0x3f3f3f3f3f3f3f3f;
for(int i=1;i<=mmax;i++){
ll now=0,js=0;
for(int j=1;j<=le;j++){
if(i&(1<<(j-1))) now+=a[j],js++;
}
if(ma[now]) ma[now]=min(ma[now],js);
//如果当前值已经出现过,元素个数取较小的
else ma[now]=js,jl[++cnt]=now;
//没有出现过,建立映射关系
if(abss(now)<ans){
ans=abss(now);
id=js;
}
//如果答案更优,更新答案
else if(abss(now)==ans){
id=min(id,js);
}
//如果答案相同,元素个数取较小的
}
//枚举前一半数,同时记录答案
sort(jl+1,jl+cnt+1);
//排序
for(int i=1;i<(1<<re);i++){
ll now=0,js=0;
for(int j=1;j<=re;j++){
if(i&(1<<(j-1))){
now+=a[j+n/2],js++;
}
}
if(abss(now)<ans){
ans=abss(now);
id=js;
}
else if(abss(now)==ans){
id=min(id,js);
}
//同上
ll aft=solve(-now);
//二分找到与相反数最接近的数
if(ans>abss(aft+now)){
ans=abss(aft+now);
id=js+ma[aft];
}//如果答案更优,更新答案
else if(ans==abss(aft+now)) id=min(id,js+ma[aft]);
//如果和之前的答案相同,个数取较小的
}
printf("%lld %d\n",ans,id);
}
return 0;
}

Subset POJ - 3977(折半枚举+二分查找)的更多相关文章

  1. CSU OJ PID=1514: Packs 超大背包问题,折半枚举+二分查找。

    1514: Packs Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 61  Solved: 4[Submit][Status][Web Board] ...

  2. poj3977(折半枚举+二分查找)

    题目链接:https://vjudge.net/problem/POJ-3977 题意:给一个大小<=35的集合,找一个非空子集合,使得子集合元素和的绝对值最小,如果有多个这样的集合,找元素个数 ...

  3. Subset---poj3977(折半枚举+二分查找)

    题目链接:http://poj.org/problem?id=3977 给你n个数,找到一个子集,使得这个子集的和的绝对值是最小的,如果有多种情况,输出子集个数最少的: n<=35,|a[i]| ...

  4. POJ 3977 折半枚举

    链接: http://poj.org/problem?id=3977 题意: 给你n个数,n最大35,让你从中选几个数,不能选0个,使它们和的绝对值最小,如果有一样的,取个数最小的 思路: 子集个数共 ...

  5. POJ 3273 Monthly Expense二分查找[最小化最大值问题]

    POJ 3273 Monthly Expense二分查找(最大值最小化问题) 题目:Monthly Expense Description Farmer John is an astounding a ...

  6. 【折半枚举+二分】POJ 3977 Subset

    题目内容 Vjudge链接 给你\(n\)个数,求出这\(n\)个数的一个非空子集,使子集中的数加和的绝对值最小,在此基础上子集中元素的个数应最小. 输入格式 输入含多组数据,每组数据有两行,第一行是 ...

  7. POJ 3977 Subset(折半枚举+二分)

    SubsetTime Limit: 30000MS        Memory Limit: 65536KTotal Submissions: 6754        Accepted: 1277 D ...

  8. POJ 2549 Sumsets(折半枚举+二分)

    Sumsets Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 11946   Accepted: 3299 Descript ...

  9. Divide and conquer:Subset(POJ 3977)

    子序列 题目大意:给定一串数字序列,要你从中挑一定个数的数字使这些数字和绝对值最小,求出最小组合数 题目的数字最多35个,一看就是要数字枚举了,但是如果直接枚举,复杂度就是O(2^35)了,显然行不通 ...

随机推荐

  1. ASP.NET Core 3.1 WebApi部署到腾讯云CentOS 7+Docker

    一.准备 首先需要有一台CentOS服务器,安装最新版Docker,配置镜像加速等,安装方法网上很多,下面是一些相关指令: yum install -y yum-utils device-mapper ...

  2. MySQL查询优化利刃-EXPLAIN

    有一个 ? 遇到这样一个疑问:当where中In一个索引字段,那么在查询中还会使用到索引吗? SELECT * FROM table_name WHERE column_index in (expr) ...

  3. Remap 后的 USART1 不能发送数据

    最近在使用 STM32 的 USART1 时,遇到一点问题,记录一下. 如果 USART1 不是使用 PA9 和 PA10 作接收和发送引脚,而是 PB6 和 PB7,那么有一点就一定得注意了:此时 ...

  4. Largest Allowed Area【模拟+二分】

    Largest Allowed Area 题目链接(点击) 题目描述 A company is looking for land to build its headquarters. It has a ...

  5. 商城04——门户网站介绍&商城首页搭建&内容系统创建&CMS实现

    1.   课程计划 1.门户系统的搭建 2.显示商城首页 3.内容管理系统的实现 a)  内容分类管理 b) 内容管理 2.   门户系统的搭建 2.1. 什么是门户系统 从广义上来说,它将各种应用系 ...

  6. 大众点评cat实时监控简介及部署

    简介 背景 CAT(Central Application Tracking)是由吴其敏(前大众点评首席架构师,现携程架构负责人)主导设计基于Java开发打造的实时应用监控平台,为大众点评网提供了全面 ...

  7. ca74a_c++__文件流对象的使用-用来读写文件ifstream

    /*ca74a_c++__文件流对象的使用-用来读写文件将文件流对象绑定到文件上检查文件是否打开成功将文件流与新文件重新绑定清楚文件流的状态infile.close();//关闭流 infile.cl ...

  8. 08.DRF-反序列化

    三.反序列化使用 3.1 验证 使用序列化器进行反序列化时,需要对数据进行验证后,才能获取验证成功的数据或保存成模型类对象. 在获取反序列化的数据前,必须调用is_valid()方法进行验证,验证成功 ...

  9. opencv C++ Mat构造函数

    cv::Scalar scalar1(v); cv::Mat mat3(size,CV_8UC1,scalar1); std::cout<<mat3<<std::endl; s ...

  10. ​Shiro授权

    Shiro三种授权方式 编程式:通过写 if/else 授权代码块完成: Subject subject = SecurityUtils.getSubject(); if(subject.hasRol ...