【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;
}