问个关于dp的问题
查看原帖
问个关于dp的问题
110694
ChengZe楼主2020/7/26 13:35

本题中我使用以下两个dp,却代码一95,代码二100(其中ok[i]为所有的符合条件的情况),初始f[0]=0,其他0x3f3f3f。

//代码一
for(int i=1;i<=tot;i++)
for(int j=0;j<(1<<n);j++)if((j&ok[i])==ok[i])
f[j]=min(f[j],f[j-ok[i]]+1);
//代码二
for(int i=1;i<=tot;i++)
for(int j=0;j<(1<<n)j++)
f[j|ok[i]]=min(f[j|ok[i]],f[j]+1);

本人认为一与二完全相同,却并非如此  以下附上完整代码

#include <bits/stdc++.h>
using namespace std;
const double eps=1e-10;
int f[(1<<18)+10],n,m,ok[20*20],tot;
double a,b,x[20],y[20];

void getab(double x1,double y1,double x2,double y2){
	a=(y1*x2-y2*x1)/(x1*x2*(x1-x2));
	b=((y1*x2*x2-y2*x1*x1)/(x1*x2*x2-x1*x1*x2));
}

double cal(double x){
	return a*x*x+b*x;
}

void doit()
{
	tot=0;memset(f,0x3f,sizeof(f));memset(ok,0,sizeof(ok));
	f[0]=0;
	scanf("%d%d",&n,&m);
	for(int i=0;i<n;i++)scanf("%lf%lf",&x[i],&y[i]);
	for(int i=0;i<n;i++)
	for(int j=i;j<n;j++){
		if(i==j){
			ok[++tot]=(1<<i); 
			continue;
		} 
		if(x[i]==x[j])continue;
		getab(x[i],y[i],x[j],y[j]);
		if(a>=-eps)continue;
		tot++;
		for(int k=0;k<n;k++){
			if(fabs(cal(x[k])-y[k])<=eps)ok[tot]|=(1<<k);
		}
	}
	for(int i=1;i<=tot;i++)
//	for(int j=0;j<(1<<n);j++)if((j&ok[i])==ok[i])f[j]=min(f[j],f[j-ok[i]]+1);
//就是这句话不同
	for(int j=0;j<(1<<n);j++)f[j|ok[i]]=min(f[j|ok[i]],f[j]+1);
	printf("%d\n",f[(1<<n)-1]);
}

int main()
{

	int T;scanf("%d",&T);
	while(T--)doit();
	return 0;
 } 
2020/7/26 13:35
加载中...