2018.7.11DP优化考试

考题：

NO1：给你一个n个数的序列，数字两两不同，要你把这些数升序排列，有两种操作：1.交换相邻的两个数。2，交换三个数。问你操作的最少次数……

1<n<100000

 1，暴力……

2，考虑操作2只会在奇数位的数之间、或偶数位的数之间进行交换。
那么将序列排序后，记录排序后的位次。
如果一个数排序后的位次和原来的奇偶性不同，那么它就需要调换。
统计出需要调换的数的数量num。

答案就是num/2。（根本不是DP，还DP优化测试呢！）

代码就算了，太简单了。

NO2:所有城市的上方有一个出发点，所有城市的下方有一个终点，旅行者们可以从出发点出 发穿越所有层城市到达终点。 这些城市有 m 层，从高往低为 1 到 m 层，每层都有 n 个城市。从出发点到第一层的第 i 个城市的花费为

 对于 30% 的数据，1 ≤

1,暴力……

2，f[i][j+b[k]]%p+=f[i][j]

第一层要特殊处理

最后一层也要

3，矩阵快速幂……

#include<bits/stdc++.h>

using namespace std;

int n,l,p;

int a[1000001],b[1000001],c[1000001],f[1000001];

struct cow

{

int a[101][101],n,m;

void init(int c,int b)

{

n=c;m=b;

memset(a,0,sizeof(a));

}

friend cow operator *(cow a,cow b)

{

cow c;c.init(a.n,b.m);

for(int i=0;i<c.n;i++)

  for(int j=0;j<c.m;j++)

    for(int k=0;k<a.m;k++)

    {

    c.a[i][j]+=1ll*a.a[i][k]*b.a[k][j]%1000000007;

    if(c.a[i][j]>=1000000007)c.a[i][j]-=1000000007;

    }

return c;

}

}qq,pp;

int sub(int x, int y) {

x = x - y;

if(x < 0) x += p;

return x;

}

int main()

{

cin>>n>>l>>p;

for(int i=1;i<=n;i++){cin>>a[i];if(a[i]>=p)a[i]=a[i]-p;}

for(int i=1;i<=n;i++){cin>>b[i];if(b[i]>=p)b[i]=b[i]-p;}

for(int i=1;i<=n;i++){cin>>c[i];if(c[i]>=p)c[i]=c[i]-p;}

    pp.init(p,p);

for(int i=0;i<p;i++)

  for(int j=1;j<=n;j++)

  pp.a[i][sub(i,b[j])]++;

qq.init(p,1);

for(int i=1;i<=n;i++)qq.a[a[i]][0]++;

int u=l-2;

cow res;res.init(pp.n,pp.m);

for(int i=0;i<pp.n;i++)res.a[i][i]=1;

while(u)

    {

    if(u&1)res=res*pp;

    pp=pp*pp;

    u>>=1;

    }

    pp=res;

    qq=pp*qq;

    for(int i=0;i<p;i++)f[i]=qq.a[i][0];

int ans=0;

    for(int i=1;i<=n;i++)

    {

    int t=c[i]+b[i];

    if(t>=p)t=t-p;

    t=-t;

    if(t<0)t+=p;

    ans+=f[t];

if(ans>=100000007)ans-=1000000007;

    }

    cout<<ans;

return 0;

}

NO3：定义一个 A 和 B 的公共子序列是“完美的”，当且仅当这个公共子序列包含了正整数序列 C。换句话说，C 是这个公共子序列的一个子串（连续）。 求 A 和 B 最长的“完美的”公共子序列的长度。 

对于 30% 的数据，1 ≤

1，暴力……

2，暴力（似乎有些欠打）

f[i][j]表示A[1~i],b[1~j]的LCS

f[i][j]表示A[n~i],b[n~j]的LCS

处理出na[i]表示A到i要进过多少到C串

处理出nb[i]表示B到i要进过多少到C串

答案……f[i][j]+g[na[i]][nb[i]]

#include<bits/stdc++.h>

using namespace std;

int n,m,k;

int a[3001],b[3001],c[3001],f[3001][3001]={},g[3001][3001]={},p[3001],q[3001];

int main()

{

cin>>n;for(int i=1;i<=n;i++)cin>>a[i];

cin>>m;for(int i=1;i<=m;i++)cin>>b[i];

cin>>k;for(int i=1;i<=k;i++)cin>>c[i];

for(int i=1;i<=m;i++)

  for(int j=1;j<=n;j++)

    if(a[i]==b[j])f[i][j]=f[i-1][j-1]+1;

          else f[i][j]=max(f[i][j-1],f[i-1][j]);

if(!k){cout<<f[m][n];return 0;}

for(int i=n;i>=1;i--)

  for(int j=m;j>=1;j--)

    if(a[i]==b[j])g[i][j]=g[i+1][j+1]+1;

          else g[i][j]=max(g[i][j+1],g[i+1][j]);

for(int i=1;i<=n;i++)p[i]=n+1;

for(int i=1;i<=m;i++)q[i]=m+1;

for(int i=1;i<=n;i++)

{

int o=1;

for(int j=i;j<=n;j++)

{

  if(a[j]==c[o])o++;

      if(o==k+1){p[i]=j;break;}

}

}

for(int i=1;i<=m;i++)

{

int o=1;

for(int j=i;j<=m;j++)

{

  if(b[j]==c[o])o++;

      if(o==k+1){q[i]=j;break;}

}

}

int ans=-1;

for(int i=1;i<=n;i++)

  for(int j=1;j<=m;j++)

if(p[i]!=n+1&&q[i]!=m+1)

  ans=max(ans,f[i-1][j-1]+k+g[p[i]+1][q[j]+1]);

cout<<ans;

return 0;

}