草稿

最新推荐文章于 2024-07-15 17:05:45 发布

Deep_Kevin

最新推荐文章于 2024-07-15 17:05:45 发布

阅读量249

点赞数

本文链接：https://blog.csdn.net/Deep_Kevin/article/details/85201287

版权

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

int n,T,R,m;
struct node{
	int p,s,t,id;
	bool operator<(const node e)const{
		return s<e.s;
	}
}q[1010];
double a[4010][1010],eps=1e-5;
double b[4010][1010];
int ans,an[1010],last[1010];
int id[2010];

void pivot(int x,int y){
	swap(id[n+x],id[y]);
	double temp=a[x][y];
	for(int i=0;i<=n;i++) a[x][i]/=temp;
	a[x][y]=1/temp;
	for(int i=0;i<=m;i++) if(x!=i){
		temp=a[i][y];a[i][y]=0;
		for(int j=0;j<=n;j++) a[i][j]-=temp*a[x][j];
	}
}

void simplex(){
	int x,y;
	double mmin;
	while(1){
		x=y=0;
		for(int i=1;i<=n;i++) if(a[0][i]>eps) {y=i;break;}
		if(!y) break;
		mmin=(double)1e18;
		for(int i=1;i<=m;i++) if(a[i][y]>eps && a[i][0]/a[i][y]<mmin) {x=i;mmin=a[i][0]/a[i][y];}
		pivot(x,y);
	}
	return ;
}

void dfs(int x){
	for(int i=0;i<=m;i++) 
		for(int j=0;j<=n;j++) 
			a[i][j]=b[i][j];
	simplex();
	if(-a[0][0]<=ans) return ;
	if(x==n+1){
		ans=(int)-a[0][0];
		for(int i=1;i<=m;i++) an[id[n+i]]=(int)a[i][0];
		return ;
	}
	m++;b[m][x]=1;b[m][0]=0;
	dfs(x+1);
	b[m][x]=-1;b[m][0]=-1;
	dfs(x+1);m--;
}

int main(){
	scanf("%d %d %d\n",&n,&T,&R);ans=R;
	for(int i=1;i<=n;i++) scanf("%d %d %d",&q[i].s,&q[i].p,&q[i].t),q[i].id=i;
	sort(q+1,q+1+n);
	a[0][0]=-R;
	for(int i=1;i<=n;i++) a[0][i]=q[i].p;
	for(int i=1;i<=n;i++){
		a[i][i]=q[i].s;
		for(int j=1;j<i;j++) a[i][j]=-q[j].p;
		a[i][0]=R;
	}
	for(int i=1;i<=n;i++) id[i]=i;
	for(int i=1;i<=n;i++) a[n+1][i]=q[i].t;a[n+1][0]=T;
	for(int i=1;i<=n;i++) a[n+1+i][i]=1,a[n+1+i][0]=1;
	m=n+n+1;
	for(int i=0;i<=m;i++)
		for(int j=0;j<=n;j++) b[i][j]=a[i][j];
	dfs(1);
	printf("%d\n",ans);
	for(int i=1;i<=n;i++) if(an[i]) last[++last[0]]=q[i].id;
	for(int i=1;i<=last[0];i++) printf("%d ",last[i]);
}

$\sum_{i=1}^{A}P^i[m],(A<=m) \\=\frac{P^{A+1}}{P-1}[m]$

$Ans(k)=\sum_{i=1}^{n}a(i)^k$

$F(x)=\sum_{i=1}^n\sum_{j=0}a_i^j*x^j \\=\sum_{i=1}^n\frac{1}{1-a_ix} \\=n-x\sum_{i=1}^n \frac{-a_i}{1-ai_x} \\=n-x\sum_{i=1}^n [\ln(1-a_ix)]' \\=n-x(\ln(\prod_{i=1}^n(1-a_ix))) \\=n-x(\frac{1}{\prod_{i=1}^n(1-a_ix)} \times (\prod_{i=1}^n(1-a_ix)')$

$\sum_{i=0}^n [\frac{a*i}{c}]=\frac{n*(n+1)}{2}[\frac{a}{c}]$

$\sum_{a=1}^n\sum_{b=1}^n\sum_{c=1}^n\sum_{d=1}^nA_aB_bC_cD_dE_{\gcd(a,b)}F_{\gcd(b,c)}G_{\gcd(c,d)}H_{gcd(a,d)} \\=\sum_{e=1}^n\sum_{f=1}^n\sum_{g=1}^n\sum_{h=1}^nE_eF_fG_gH_h\sum_{a=1}^{\frac{n}{[e,h]}}\sum_{b=1}^{\frac{n}{[e,f]}}\sum_{c=1}^{\frac{n}{[f,g]}}\sum_{d=1}^{\frac{n}{[g,h]}}A_{a*[e,h]}B_{b*[e,f]}C_{c*[f,g]}D_{d*[g,h]} \\=$