首页 > 代码库 > 动态规划——邮局、区间、坐标
动态规划——邮局、区间、坐标
施工中ing...
问题 A: P1045
时间限制: 1 Sec 内存限制: 128 MB题目描述
题目很简单,给出N个数字,不改变它们的相对位置,在中间加入K个乘号和N-K-1个加号,(括号随便加)使最终结果尽量大。因为乘号和加号一共就是N-1个了,所以恰好每两个相邻数字之间都有一个符号。例如: N=5, K=2,5个数字分别为1、2、3、4、5,可以加成: 1*2*(3+4+5)=24 1*(2+3)*(4+5)=45 (1*2+3)*(4+5)=45 ……
输入
输入文件共有二行,第一行为两个有空格隔开的整数,表示N和K,其中(2< =N< =15, 0< =K< =N-1)。第二行为 N个用空格隔开的数字(每个数字在0到9之间)。
输出
输出文件仅一行包含一个整数,表示要求的最大的结果 最后的结果< =maxlongint
样例输入
样例输出
提示
对于30%的数据,N< = 10; 对于全部的数据,N < = 100。
邮局动规。
我们可以构造出他的结构特征f(i,j)表示前i个数用了j个乘号的最大结果。
然后是递归式:f(i,j)=f(k,j-1)*Sum(k+1,i);
1 #include<iostream> 2 #include<algorithm> 3 #include<cstdio> 4 #include<cstring> 5 #include<cmath> 6 #include<cstdlib> 7 #include<vector> 8 using namespace std; 9 typedef long long ll; 10 typedef long double ld; 11 typedef pair<int,int> pr; 12 const double pi=acos(-1); 13 #define rep(i,a,n) for(int i=a;i<=n;i++) 14 #define per(i,n,a) for(int i=n;i>=a;i--) 15 #define Rep(i,u) for(int i=head[u];i;i=Next[i]) 16 #define clr(a) memset(a,0,sizeof(a)) 17 #define pb push_back 18 #define mp make_pair 19 #define fi first 20 #define sc second 21 #define pq priority_queue 22 #define pqb priority_queue <int, vector<int>, less<int> > 23 #define pqs priority_queue <int, vector<int>, greater<int> > 24 #define vec vector 25 ld eps=1e-9; 26 ll pp=1000000007; 27 ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} 28 ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} 29 void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); } 30 //void add(int x,int y,int z){ v[++e]=y; next[e]=head[x]; head[x]=e; cost[e]=z; } 31 int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1}; 32 ll read(){ ll ans=0; char last=‘ ‘,ch=getchar(); 33 while(ch<‘0‘ || ch>‘9‘)last=ch,ch=getchar(); 34 while(ch>=‘0‘ && ch<=‘9‘)ans=ans*10+ch-‘0‘,ch=getchar(); 35 if(last==‘-‘)ans=-ans; return ans; 36 } 37 #define N 1000 38 ll a[N],f[N][N],s[N]; 39 int main() 40 { 41 int n=read(),K=read(); 42 for (int i=1;i<=n;i++) 43 a[i]=read(); 44 for (int i=1;i<=n;i++) 45 s[i]=s[i-1]+a[i]; 46 for (int i=1;i<=n;i++) f[i][0]=s[i]; 47 for (int i=2;i<=n;i++){ 48 for (int j=1;j<=min(i-1,K);j++) 49 for (int k=1;k<=i-1;k++) 50 f[i][j]=max(f[i][j],f[k][j-1]*(s[i]-s[k])); 51 } 52 printf("%lld",f[n][K]); 53 return 0; 54 }
问题 B: P1047
时间限制: 1 Sec 内存限制: 128 MB题目描述
今年是国际数学联盟确定的“2000——世界数学年”,又恰逢我国著名数学家华罗庚先生诞辰90周年。在华罗庚先生的家乡江苏金坛,组织了一场别开生面的数学智力竞赛的活动,你的一个好朋友XZ也有幸得以参加。活动中,主持人给所有参加活动的选手出了这样一道题目: 设有一个长度N的数字串,要求选手使用K个乘号将它分成K+1个部分,找出一种分法,使得这K+1个部分的乘积能够为最大。 同时,为了帮助选手能够正确理解题意,主持人还举了如下的一个例子: 有一个数字串: 312,当N=3,K=1时会有以下两种分法: 1)3*12=36 2)31*2=62 这时,符合题目要求的结果是: 31*2=62 现在,请你帮助你的好朋友XZ设计一个程序,求得正确的答案。
输入
程序的输入共有两行: 第一行共有2个自然数N,K (6< =N< =40,0< =K< =5) 第二行是一个K度为N的数字串。
输出
结果输出到文件,相对于输入,应输出所求得的最大乘积(一个自然数)。
样例输入
样例输出
我们可以构造出他的结构特征f(i,j)表示前i个数用了j个乘号的最大结果。
然后是递归式:f(i,j)=f(k,j-1)*S(k+1,i);
1 #include<iostream> 2 #include<algorithm> 3 #include<cstdio> 4 #include<cstring> 5 #include<cmath> 6 #include<cstdlib> 7 #include<vector> 8 using namespace std; 9 typedef long long ll; 10 typedef long double ld; 11 typedef pair<int,int> pr; 12 const double pi=acos(-1); 13 #define rep(i,a,n) for(int i=a;i<=n;i++) 14 #define per(i,n,a) for(int i=n;i>=a;i--) 15 #define Rep(i,u) for(int i=head[u];i;i=Next[i]) 16 #define clr(a) memset(a,0,sizeof(a)) 17 #define pb push_back 18 #define mp make_pair 19 #define fi first 20 #define sc second 21 #define pq priority_queue 22 #define pqb priority_queue <int, vector<int>, less<int> > 23 #define pqs priority_queue <int, vector<int>, greater<int> > 24 #define vec vector 25 ld eps=1e-9; 26 ll pp=1000000007; 27 ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} 28 ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} 29 //void add(int x,int y,int z){ v[++e]=y; next[e]=head[x]; head[x]=e; cost[e]=z; } 30 int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1}; 31 ll read(){ ll ans=0; char last=‘ ‘,ch=getchar(); 32 while(ch<‘0‘ || ch>‘9‘)last=ch,ch=getchar(); 33 while(ch>=‘0‘ && ch<=‘9‘)ans=ans*10+ch-‘0‘,ch=getchar(); 34 if(last==‘-‘)ans=-ans; return ans; 35 } 36 #define N 100 37 struct Hp{ 38 int len,nu[1000]; 39 Hp(){ 40 len=0; memset(nu,0,sizeof(nu)); 41 } 42 friend Hp operator +(Hp a,Hp b){ 43 Hp c; 44 if (a.len<b.len) c.len=b.len; else c.len=a.len; 45 for (int i=1;i<=c.len;i++){ 46 c.nu[i]+=a.nu[i]+b.nu[i]; 47 if (c.nu[i]>9) c.nu[i+1]++,c.nu[i]-=10; 48 } 49 if (c.nu[c.len]>0) c.len++; 50 return c; 51 } 52 friend Hp operator *(Hp a,Hp b){ 53 Hp c; 54 for (int i=1;i<=a.len;i++) 55 for (int j=1;j<=b.len;j++){ 56 c.nu[i+j-1]+=a.nu[i]*b.nu[j]; 57 } 58 c.len=a.len+b.len; 59 for (int i=1;i<c.len;i++){ 60 if (c.nu[i]>9) c.nu[i+1]+=c.nu[i]/10,c.nu[i]%=10; 61 } 62 while (c.nu[c.len]==0&&c.len>1) c.len--; 63 return c; 64 } 65 void out(){ 66 for (int i=len;i>0;i--) 67 printf("%d",nu[i]); 68 puts(""); 69 } 70 }; 71 Hp max_(Hp a,Hp b){ 72 if (a.len<b.len) return b; 73 if (a.len>b.len) return a; 74 for (int i=a.len;i>0;i--){ 75 if (a.nu[i]>b.nu[i]) return a; 76 if (a.nu[i]<b.nu[i]) return b; 77 } 78 return a; 79 } 80 Hp s[N][N],a[N],f[N][N]; 81 char str[N]; 82 int main() 83 { 84 int n=read(),K=read(); 85 scanf("%s",str+1); 86 for (int i=1;i<=n;i++) a[i].len=1,a[i].nu[1]=str[i]-‘0‘; 87 Hp t; t.len=2; t.nu[1]=0; t.nu[2]=1; 88 for (int i=1;i<=n;i++) 89 for (int j=i;j<=n;j++) s[i][j]=s[i][j-1]*t+a[j]; 90 for (int i=1;i<=n;i++) f[i][0]=s[1][i]; 91 for (int i=1;i<=n;i++) 92 for (int j=1;j<=min(i-1,K);j++){ 93 for (int k=1;k<=i-1;k++) 94 f[i][j]=max_(f[i][j],f[k][j-1]*s[k+1][i]); 95 } 96 f[n][K].out(); 97 return 0; 98 }
问题 C: 括号序列
时间限制: 1 Sec 内存限制: 128 MB题目描述
定义如下规则序列(字符串):
1.空序列是规则序列;
2.如果S是规则序列,那(S)和[S]也是规则序列;
3.如果A和B都是规则序列,那么AB也是规则序列。
例如,下面的字符串都是规则序列:
(), [], (()), ([]), ()[], ()[()]
这几个不是规则序列:
(, [, ], )(, ([()
现在,给出一些有‘(‘ , ‘)‘ , ‘[‘ , ‘]‘组成的序列,请添加尽量少的括号,得到一个规则序列,并输出该序列的长度。
输入
输入一个有‘(‘ , ‘)‘ , ‘[‘ , ‘]‘组成的序列,长度不超过200
输出
输出规则后的字串长度
样例输入
样例输出
提示
([()])//6
区间动规。
我们可以构造出他的结构特征f(i,j)表示i~j最少需要添多少括号。
然后是递归式:f(i,j)=min(f(i,k)+f(k+1,j));如果两边匹配的话就=min(f[i+1][j-1]);
1 #include<iostream> 2 #include<algorithm> 3 #include<cstdio> 4 #include<cstring> 5 #include<cmath> 6 #include<cstdlib> 7 #include<vector> 8 using namespace std; 9 typedef long long ll; 10 typedef long double ld; 11 typedef pair<int,int> pr; 12 const double pi=acos(-1); 13 #define rep(i,a,n) for(int i=a;i<=n;i++) 14 #define per(i,n,a) for(int i=n;i>=a;i--) 15 #define Rep(i,u) for(int i=head[u];i;i=Next[i]) 16 #define clr(a) memset(a,0,sizeof(a)) 17 #define pb push_back 18 #define mp make_pair 19 #define fi first 20 #define sc second 21 #define pq priority_queue 22 #define pqb priority_queue <int, vector<int>, less<int> > 23 #define pqs priority_queue <int, vector<int>, greater<int> > 24 #define vec vector 25 ld eps=1e-9; 26 ll pp=1000000007; 27 ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} 28 ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} 29 void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); } 30 //void add(int x,int y,int z){ v[++e]=y; next[e]=head[x]; head[x]=e; cost[e]=z; } 31 int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1}; 32 ll read(){ ll ans=0; char last=‘ ‘,ch=getchar(); 33 while(ch<‘0‘ || ch>‘9‘)last=ch,ch=getchar(); 34 while(ch>=‘0‘ && ch<=‘9‘)ans=ans*10+ch-‘0‘,ch=getchar(); 35 if(last==‘-‘)ans=-ans; return ans; 36 } 37 int f[300][300]; 38 char str[300]; 39 bool match(int i,int j){ 40 if (str[i]==‘(‘ && str[j]==‘)‘) return 1; 41 if (str[i]==‘[‘ && str[j]==‘]‘) return 1; 42 return 0; 43 } 44 int main() 45 { 46 scanf("%s",str+1); int n=strlen(str+1); 47 for (int i=1;i<=n;i++) 48 for (int j=1;j<=n;j++) f[i][j]=1e+08; 49 for (int i=1;i<=n;i++) f[i][i]=1,f[i][i-1]=0; 50 for (int p=2;p<=n;p++) 51 for (int i=1;i+p-1<=n;i++){ 52 int j=i+p-1; 53 if (match(i,j)) f[i][j]=min(f[i][j],f[i+1][j-1]); 54 for (int k=i;k<j;k++) 55 f[i][j]=min(f[i][j],f[i][k]+f[k+1][j]); 56 } 57 printf("%d",n+f[1][n]); 58 return 0; 59 }
问题 D: 加分二叉树
时间限制: 1 Sec 内存限制: 128 MB题目描述
设一个n个节点的二叉树tree的中序遍历为(l,2,3,…,n),其中数字1,2,3,…,n为节点编号。每个节点都有一个分数(均为正整数),记第j个节点的分数为di,tree及它的每个子树都有一个加分,任一棵子树subtree(也包含tree本身)的加分计算方法如下:
subtree的左子树的加分× subtree的右子树的加分+subtree的根的分数
若某个子树为主,规定其加分为1,叶子的加分就是叶节点本身的分数。不考虑它的空
子树。
试求一棵符合中序遍历为(1,2,3,…,n)且加分最高的二叉树tree。要求输出;
(1)tree的最高加分
(2)tree的前序遍历
输入
第1行:一个整数n(n<30),为节点个数。
第2行:n个用空格隔开的整数,为每个节点的分数(分数<100)。
输出
第1行:一个整数,为最高加分(结果不会超过4,000,000,000)。
第2行:n个用空格隔开的整数,为该树的前序遍历。
样例输入
样例输出
1 #include<iostream> 2 #include<algorithm> 3 #include<cstdio> 4 #include<cstring> 5 #include<cmath> 6 #include<cstdlib> 7 #include<vector> 8 using namespace std; 9 typedef long long ll; 10 typedef long double ld; 11 typedef pair<int,int> pr; 12 const double pi=acos(-1); 13 #define rep(i,a,n) for(int i=a;i<=n;i++) 14 #define per(i,n,a) for(int i=n;i>=a;i--) 15 #define Rep(i,u) for(int i=head[u];i;i=Next[i]) 16 #define clr(a) memset(a,0,sizeof(a)) 17 #define pb push_back 18 #define mp make_pair 19 #define fi first 20 #define sc second 21 #define pq priority_queue 22 #define pqb priority_queue <int, vector<int>, less<int> > 23 #define pqs priority_queue <int, vector<int>, greater<int> > 24 #define vec vector 25 ld eps=1e-9; 26 ll pp=1000000007; 27 ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} 28 ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} 29 void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); } 30 //void add(int x,int y,int z){ v[++e]=y; next[e]=head[x]; head[x]=e; cost[e]=z; } 31 int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1}; 32 ll read(){ ll ans=0; char last=‘ ‘,ch=getchar(); 33 while(ch<‘0‘ || ch>‘9‘)last=ch,ch=getchar(); 34 while(ch>=‘0‘ && ch<=‘9‘)ans=ans*10+ch-‘0‘,ch=getchar(); 35 if(last==‘-‘)ans=-ans; return ans; 36 } 37 #define N 100 38 ll f[N][N],a[N]; 39 int s[N][N],n,num; 40 void dfs(int i,int j){ 41 if (i>j) return; 42 int x=s[i][j]; num++; 43 if (num<n) printf("%d ",x); 44 else printf("%d\n",x); 45 dfs(i,x-1); 46 dfs(x+1,j); 47 } 48 int main() 49 { 50 n=read(); 51 for (int i=1;i<=n;i++) a[i]=read(); 52 for (int i=1;i<=n;i++) f[i][i]=a[i],s[i][i]=i; 53 for (int i=0;i<=n;i++) f[i+1][i]=1; 54 for (int p=2;p<=n;p++) 55 for (int i=1;i+p-1<=n;i++){ 56 int j=i+p-1; 57 for (int k=i;k<=j;k++) 58 { 59 ll sum=f[i][k-1]*f[k+1][j]+a[k]; 60 if (sum>f[i][j]){ 61 f[i][j]=sum; s[i][j]=k; 62 } 63 } 64 } 65 printf("%lld\n",f[1][n]); 66 dfs(1,n); 67 return 0; 68 }
问题 E: 小明的喷漆计划
时间限制: 5 Sec 内存限制: 128 MB题目描述
小明极其喜欢涂鸦,总是在墙上涂上各种颜色的漆。现在小明得到一个任务,需要喷涂一段空白围墙,且单位长度内的颜色都是相同的。小明有一种喷涂工具,它可以给任意长度的一段墙面涂上任意颜色的漆,这样的操作计为一次操作。小明要完成这个任务,又想使得操作次数尽量少,就请你帮他解决这个问题吧。
输入
有多组输入数据。
每组包含一个长度不超过100的字符串(均由小写字母组成),代表需要涂鸦的墙面目标状态
输出
至少需要几次操作,可达到目标
样例输入
样例输出
(str[i]==str[j]) {
f[i][j]=min(f[i][j],f[i][j-1]);f[i][j]=min(f[i][j],f[i-1][j]);}1 #include<iostream> 2 #include<algorithm> 3 #include<cstdio> 4 #include<cstring> 5 #include<cmath> 6 #include<cstdlib> 7 #include<vector> 8 using namespace std; 9 typedef long long ll; 10 typedef long double ld; 11 typedef pair<int,int> pr; 12 const double pi=acos(-1); 13 #define rep(i,a,n) for(int i=a;i<=n;i++) 14 #define per(i,n,a) for(int i=n;i>=a;i--) 15 #define Rep(i,u) for(int i=head[u];i;i=Next[i]) 16 #define clr(a) memset(a,0,sizeof(a)) 17 #define pb push_back 18 #define mp make_pair 19 #define fi first 20 #define sc second 21 #define pq priority_queue 22 #define pqb priority_queue <int, vector<int>, less<int> > 23 #define pqs priority_queue <int, vector<int>, greater<int> > 24 #define vec vector 25 ld eps=1e-9; 26 ll pp=1000000007; 27 ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} 28 ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} 29 void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); } 30 //void add(int x,int y,int z){ v[++e]=y; next[e]=head[x]; head[x]=e; cost[e]=z; } 31 int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1}; 32 ll read(){ ll ans=0; char last=‘ ‘,ch=getchar(); 33 while(ch<‘0‘ || ch>‘9‘)last=ch,ch=getchar(); 34 while(ch>=‘0‘ && ch<=‘9‘)ans=ans*10+ch-‘0‘,ch=getchar(); 35 if(last==‘-‘)ans=-ans; return ans; 36 } 37 int f[200][200]; 38 char str[200]; 39 void work() 40 { 41 int n=strlen(str+1); 42 for (int i=0;i<=n;i++) 43 for (int j=0;j<=n;j++) f[i][j]=pp; 44 for (int p=1;p<=n;p++) 45 for (int i=1;i+p-1<=n;i++){ 46 int j=i+p-1; 47 if (p==1) f[i][j]=1; 48 else{ 49 if (str[i]==str[j]){ 50 f[i][j]=min(f[i][j],f[i][j-1]); 51 f[i][j]=min(f[i][j],f[i-1][j]); 52 } 53 else { 54 for (int k=i;k<j;k++) 55 f[i][j]=min(f[i][j],f[i][k]+f[k+1][j]); 56 } 57 } 58 } 59 cout<<f[1][n]<<endl; 60 } 61 int main(){ 62 while (~scanf("%s",str+1)){ 63 work(); 64 } 65 }
动态规划——邮局、区间、坐标