BZOJ3550 ONTAK2010 Vacation

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/hanyuweining/article/details/84852426

传送门（没有是因为我没有权限号）

单纯形暴力构造吼啊

对每个元素限制一下不能选多于一个 再对每N个元素限制一下不能多于K个

一共是3N+2N-1=5N-1个约束

跑就行了2333

（单纯形在不需要ID的时候尽量不要写 因为需要开N+M的空间 RE惨惨 非常容易忘）

附代码。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 20021225
#define ll long long
#define mxn 210
#define db double
#define eps 1e-8
using namespace std;

db M[mxn*6][mxn*3];
int n,m;//id[mxn*6];
void privot(int x,int y)
{
	//swap(id[n+x],id[y]);
	db tmp=1.0/M[x][y];M[x][y]=1.0;
	for(int i=0;i<=m;i++)	M[x][i]*=tmp;
	for(int i=0;i<=n;i++)
	{
		if(i==x||abs(M[i][y])<eps)	continue;
		db t=M[i][y];M[i][y]=0;
		for(int j=0;j<=m;j++)
			M[i][j]-=M[x][j]*t;
	}
}

bool simplex()
{
	while(1)
	{
		int x=0,y=0;db mn=1e15;
		for(int i=1;i<=m;i++)	if(M[0][i]>eps){y=i;break;}
		if(!y)	return true;
		for(int i=1;i<=n;i++)
			if(M[i][y]>eps&&M[i][0]/M[i][y]<mn)
				mn=M[i][0]/M[i][y],x=i;
		if(!x){printf("Unbounded\n");return false;}
		privot(x,y);
	}
	return true;
}

int main()
{
	int nn,nnn,k;
	//freopen("9.in","r",stdin);
	scanf("%d%d",&nn,&k);nnn=3*nn;n=5*nn+1;m=nnn;
	for(int i=1;i<=nnn;i++)	M[i][i]=M[i][0]=1.0;//3n
	for(int i=1;i<=nnn;i++)
	{
		scanf("%lf",&M[0][i]);
		if(i+nn-1<=nnn)
		{
			for(int j=0;j<nn;j++)
				M[i+nnn][i+j]=1.0;
			M[i+nnn][0]=k;
		}
	}
	if(simplex())
		printf("%.0lf\n",-M[0][0]);
	return 0;
}