博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
证明 poj 1014 模优化修剪,部分递归 有错误
阅读量:7101 次
发布时间:2019-06-28

本文共 3294 字,大约阅读时间需要 10 分钟。

这个问题是存在做。我发现即使是可行的一个问题,但不一定正确。

大部分数据疲软,因为主题。

题目大意:有6堆石头,权重分别为1 2 3 4 5 6,要求输入 每堆个数 ,求能否够平分石头使得两堆价值同样。

网上对这道题的做法就两种,当中有错误的版本号。却也能够AC。

起初这让我等菜鸟感慨代码的简洁,但无法得出正确性的证明

接下来就对两种方法的错误性进行证明。

1.多重背包

#include #include
#include
#include
#include
#include
using namespace std;#define MAXN 100+60000int v[MAXN];int a[MAXN/3];int b[7] = {1, 60, 30, 20, 15, 12, 10};int N,T,n,sum;/*int direct[4][2]={-1,0,1,0,0,-1,0,1}; int dp[210][210][210];*/int max(int a,int b){ return a>b?a:b;}int main(){ int i,j,k,flag,casenum; casenum=0; while(1) { casenum++; memset(v,0,sizeof(v)); flag=0; sum=0; n=0; for(i=0;i<6;i++) { scanf("%d",&k); if(k) flag=1; k=k%b[i+1]; sum+=k*(i+1); for(j=1;j<=k;j++) a[n++]=i+1; } if(flag==0) break; if(sum&1) { printf("Collection #%d:\nCan't be divided.\n\n",casenum); continue; } flag=0; sum/=2; v[0]=1; for(i=0;i
=a[i];j--) { v[j]+= v[j-a[i]]; if(v[sum]) { flag=1;break; } if(flag) break; } } if(flag) printf("Collection #%d:\nCan be divided.\n\n",casenum); else printf("Collection #%d:\nCan't be divided.\n\n",casenum); } return 0;}

状态定义的是有几种方法能够转到这里来

k=k%b[i+1];

这句是一种优化,起初看到,认为非常奇妙,但并不理解为什么能够这样做。

后来证明是错误的,证明例如以下:

取模优化是错误的,以下证明优化一堆的情况

1.1a+2b+3c+4d+5e+6f
2.60*m*t+     1a+2b+3c+4d+5e+6f(t是某堆石子的个数,m是某堆石子的权重)
证明优化正确即证明1 是 2 式是充分必要条件
当1成立时候,自然得到2成立(60能够分到两堆)
当2成立有两种情况,
第一种情况,2可分,1的部分本身可分,那么60*m*t 这部分本来分掉就好
另外一种情况。2可分。1的部分本身不可分,须要将60*m*t这部分拆解分到两人才可行
由此得证将某个拆分掉是不可行的,可是不排除每堆都优化会遇到碰巧可行的情况
最后举个样例给大家
1. 0 0 0 0 66 5 -> 0 0 0 0 6 5   ture
2. 60 0 0 0 0 1 -> 0 0 0 0 0 1   fault

优化还是用2进制的方法优化吧(1,2,4,...,2^(k-1),n[i]-2^k+1,且k是满足n[i]-2^k+1>0的最大整数。

比如。假设n[i]为13。就将这样的物品分成系数分别为1,2,4,6的四件物品)

  1. 为何网上有些转移方程为v[i][j]=max(v[i-1][j],v[j-a[i]]+a[i])? 
  2. 答:能够看到j-a[i]表明与a[i]互补的状态,事实上为j,从全部的J角度来看。并未改变,这是v[0]=0 

2. dfs版本号(转载于大牛Blog)

//Memory Time //452K 0MS /*DFS*/#include
using namespace std;int n[7]; //价值为i的物品的个数int SumValue; //物品总价值int HalfValue; //物品平分价值bool flag; //标记能否平分SumValuevoid DFS(int value,int pre){ if(flag) return; if(value==HalfValue) { flag=true; return; } for(int i=pre;i>=1;i--) { if(n[i]) { if(value+i<=HalfValue) { n[i]--; DFS(value+i,i); if(flag) break; } } } return;}int main(int i){ int test=1; while(cin>>n[1]>>n[2]>>n[3]>>n[4]>>n[5]>>n[6]) { SumValue=0; //物品总价值 for(i=1;i<=6;i++) SumValue+=i*n[i]; if(SumValue==0) break; if(SumValue%2) //sum为奇数,无法平分 { cout<<"Collection #"<
<<':'<
这个版本号dfs写的非常好,当中
这个深度优先有两个长处值得思量

1.为什么没有回溯。而是直接减去了数量n[i]--;

答:两个人选择,必定是将这部分分为两份,假设不选择到最接近的数字,那剩余的则是更接近的

2.从大到小选择?

答:可能有多个小的能够用一个大的数字直接替换掉

----------------------------------------------------------------------------------------------------------------------------------

可是存在问题。

其本质是使用了贪心的策略,但无法满足有些“跳跃”的要求

eg: 0 0 3 0 3 1  须要选取的数字是不连续的,事实上还要有回溯的。

避免这个问题能够用这个版本号

void divide(int cur_value, int cur_index)  {          // set break point          if (flag)                  return;          if (cur_value == half_value)          {                  flag = true;                  return;          }          if (cur_value > half_value || cur_index >= max_index)                  return;          divide(cur_value+array[cur_index], cur_index+1);          divide(cur_value, cur_index+1);  }

看来学习或谨慎小心,信的过程

版权声明:本文博客原创文章,博客,未经同意,不得转载。

你可能感兴趣的文章
java简单统计.java文件中的有效代码行,空行,注释行
查看>>
Android面试题整理
查看>>
mysql中的主从复制slave-skip-errors参数使用方法
查看>>
Linux安装JIRA6.3.6以及安装破解汉化插件
查看>>
一个HTTP需要经过哪些步骤
查看>>
Finite State Transducers
查看>>
慧聪电子网战略升级 玩转电子产业供应链服务之道
查看>>
Javascript定时器(三)——setTimeout(func, 0)
查看>>
Git基础入门(七)Git撤销操作和远程仓库管理
查看>>
以毒攻毒?牛津大学研究人员用VR治愈被迫害妄想症
查看>>
巧用Powercfg命令 - 玩转Windows 7中的电源管理
查看>>
Java工具创建密钥库,用于Unity 3D打包、签名、发布
查看>>
《你不知道的JavaScript》整理(二)——this
查看>>
提升windows 2000的启动速度
查看>>
iftop工具
查看>>
java-第二章-华氏温度转摄氏温度
查看>>
html查看器android
查看>>
从零打造B/S 自动化运维平台 (一、自动化运维平台的应用及业务流程)
查看>>
SQL Server 2012 AlwaysOn高可用配置之三:安装“故障转移群集”功能
查看>>
shell中使用FTP
查看>>