【Plan 第五场】

超级大滚粗，不想说了。。。

A

第一眼：nima这不裸的最大全闭合子图么

然后就进坑了。。。。。。。。。。再也没跳出来

经vfk认证，此方法转化为了一个NPC问题，因此，如果“最大获利”问题的代价和收益反过来的话就不可做，除非是二分图。。。

二分图：http://blog.csdn.net/willinglive/article/details/42234539

要是想到用DP就很水了

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define rep(i,l,r) for(int i=(l),_=(r);i<=_;i++)
#define per(i,r,l) for(int i=(r),_=(l);i>=_;i--)
#define MS(arr,x) memset(arr,x,sizeof(arr))
#define INE(i,u) for(int i=head[u];~i;i=e[i].next)
#define LL long long
inline const int read()
{int r=0,k=1;char c=getchar();for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;
for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;}

int n,a[710],g;
int dp[10010];

int gcd(int a,int b){return b?gcd(b,a%b):a;}
void MIN(int &a,int b){if(a>b)a=b;}

void input()
{
	rep(i,1,n=read()) a[i]=read();
	g=a[1]; rep(i,2,n) g=gcd(g,a[i]);
}
void solve()
{
	MS(dp,1);
	dp[0]=0;
	rep(i,1,n)
	{
		per(j,10000,0)
		{
			MIN(dp[gcd(j,a[i])],dp[j]+1);
			//printf("dp[%d][%d]=%d\n",i+1,gcd(j,a[i+1]),dp[gcd(j,a[i])]);
		}
	}
	printf("%d\n",n-dp[g]);
}

int main()
{
    freopen("gcd.in","r",stdin); freopen("gcd.out","w",stdout);
    input(),solve();
    return 0;
}

B

迭代，把dp[n][1]的公式求出来

被卡精度丢了20分、、、、

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define rep(i,l,r) for(int i=(l),_=(r);i<=_;i++)
#define per(i,r,l) for(int i=(r),_=(l);i>=_;i--)
#define MS(arr,x) memset(arr,x,sizeof(arr))
#define INE(i,u) for(int i=head[u];~i;i=e[i].next)
#define LL long long
inline const int read()
{int r=0,k=1;char c=getchar();for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;
for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;}

int n,m;
double dp[1010][1010];



void input()
{
    n=read(); m=read();
}
void solve()
{
	dp[1][1]=1;
	rep(i,2,n)
	{
		double a=0.5,b=0;
		rep(j,2,i)
		{
			a/=2;
			b=b/2+dp[i-1][j-1]/3;
		}
		dp[i][1]=(1/6.0+b/2)/(1-a);
		rep(j,2,i) dp[i][j]=dp[i][j-1]/2+dp[i-1][j-1]/3;
	}
	printf("%.9f\n",dp[n][m]);
}

int main()
{
    freopen("game.in","r",stdin); freopen("game.out","w",stdout);
    input(),solve();
    return 0;
}

C

一道很神奇的DP，题目很简单

但是蕴含了各种坑。。。

要解决后效性的问题，然后滚动数组一些奇怪的问题，还有如何利用树状数组求最值，一堆转化过后就变得容易了。。。

估计现在让我做还是做不出来QAQ

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <set>
using namespace std;
#define rep(i,l,r) for(int i=(l),_=(r);i<=_;i++)
#define per(i,r,l) for(int i=(r),_=(l);i>=_;i--)
#define MS(arr,x) memset(arr,x,sizeof(arr))
#define INE(i,u) for(int i=head[u];~i;i=e[i].next)
#define LL long long
inline const int read()
{int r=0,k=1;char c=getchar();for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;
for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;}

const int N=1010;
int n,h[N];
multiset<pair<int,int> >S;
int t[N];
int bit[N][N];

void MAX(int &a,int b){if(a<b)a=b;}
void insert(int *c,int o,int x){for(;o<=n;o+=o&-o)MAX(c[o],x);}
int getmax(int *c,int o){int s=0;for(;o;o-=o&-o)MAX(s,c[o]);return s;}

void input()
{
    rep(i,1,n=read()) S.insert(make_pair(read(),i));
    int lst=-1,cnt=0;
    for(multiset<pair<int,int> >::iterator it=S.begin();it!=S.end();it++)
    {
    	if(lst==it->first) h[it->second]=cnt;
    	else h[it->second]=++cnt;
    	lst=it->first;
    }
}
void solve()
{
	rep(i,1,n)
	{
		rep(j,1,n) t[j]=getmax(bit[j],h[i])+1;
		//rep(k,1,n) t2[k]=getmax(bit[k],h[i])+1;
		
		rep(j,1,n) insert(bit[j],h[i],t[j]);
		rep(k,1,n) insert(bit[h[i]],k,t[k]);
	}
	int ans=0;
	rep(j,1,n) MAX(ans,getmax(bit[j],n));
	printf("%d\n",ans);
}

int main()
{
    freopen("music.in","r",stdin); freopen("music.out","w",stdout);
    input(),solve();
    return 0;
}