习题日常第五练

1矩阵取数游戏（题目链接）

这个题的数据确实大，一开始准备用高精度，看过别人的题解以后发现int128也可以过，只是它的读入和输出需要另外的函数来完成。首先把2的幂次方都算出来备用，然后进行动态规划，因为只能从头或者尾取数，所以进行动态规划的时候，现在的长度代表的和等于头乘以相应的2的幂次方加上剩余和与尾乘以相应2的幂次方加上剩余和之中大的那一个。具体代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
using namespace std;
__int128 a[90][90],b[90][90],dp[100][100],zp[1000],lin[1000];
inline int read()
{
    int s=0,f=1;
    char c=getchar();
    while(c<'0' || c>'9')
    {
        if(c=='-') f=-1;
        c=getchar();
    }
    while(c>='0' && c<='9')
    {
        s=s*10+c-48;
        c=getchar();
    }
    return s*f;
}
void Print(__int128 x)
{
    if(x==0) return;
    else if(x) Print(x/10);
    putchar(x%10+'0');
}
int main()
{
    ios::sync_with_stdio(false);
    int n,m,i,j,k;
    __int128 ans=0;
    n=read();
    m=read();
    zp[0]=1;
    for(i=1;i<=m;i++)
        zp[i]=zp[i-1]*2;
    for(i=0;i<n;i++)
    {
        memset(dp,0,sizeof(dp));
        for(j=0;j<m;j++)
        {
            a[i][j]=read();
            dp[j][j]=a[i][j]*zp[m];
        }
        for(k=1;k<=m-1;k++)
        {
            for(j=0;j+k<m;j++)
            {
                dp[j][j+k]=max(dp[j+1][j+k]+zp[m-k]*a[i][j],dp[j][j+k-1]+zp[m-k]*a[i][j+k]);
            }
        }
        ans+=dp[0][m-1];
    }
    if(ans==0)
        cout<<0;
    else
    Print(ans);
  return 0;
}

2邮票面值设计（题目链接）

这个题一开始确实无从下手，知道这是个搜索题，但是不会如何确定1到max 之间每一个数字都要能满足条件。看过别人的题解后发现可以用动态规划来维护每一个值所需要的最小票数。每一次搜索的值应当是上次选过的值加一到前面能拼成的最大值加一，搜索的目标就呼之欲出了。具体代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#define ll long long int
const int mod=1e9+7;
using namespace std;
int n,k,dp[50100],a[100],b[100],ans[100],anp=0;
int dpe(int sum,int de)
{
    memset(dp,0x3f3f,sizeof(dp));
    dp[0]=0;
    for(int i=1;i<=de;i++)
    {
        for(int j=b[i];j<=n*sum;j++)
        {
            dp[j]=min(dp[j],dp[j-b[i]]+1);
        }
    }
    for(int i=1;i<=n*sum;i++)
    {
        if(dp[i]>n)
            return i-1;
    }
    return n*sum;
}
void dfs(int lis1,int lis2,int de)
{
    if(de>k)
    {
        if(lis2>anp)
        {
            anp=lis2;
            for(int i=1;i<=k;i++)
                ans[i]=b[i];
        }
        return;
    }
    int i,j;
    for(i=lis1+1;i<=lis2+1;i++)
    {
        b[de]=i;
        int p=dpe(i,de);
        dfs(i,p,de+1);
    }
    return;
}
int main()
{
    ios::sync_with_stdio(false);
    int i,j;
    cin>>n>>k;
    dfs(0,0,1);
    for(i=1;i<=k;i++)
        cout<<ans[i]<<" ";
    cout<<endl<<"MAX="<<anp;
  return 0;
}

3产生数（题目链接）

这个题一开始我想用并查集求每个数字能替换几次，后来发现是有方向的，所以就改用floyd算法求出每个数字能被替换几次。数据比较大，所以用高精度来存最终的值。具体代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#define ll long long int
const int mod=1e9+7;
using namespace std;
int n,k,dp[50100],a[100],b[100],ans[100],anp=0,vis[20][20];
int sum[200];
void cheng(int x)
{
    int i;
    for(i=0;i<100;i++)
        sum[i]*=x;
    for(i=0;i<100;i++)
    {
        sum[i+1]+=sum[i]/10;
        sum[i]%=10;
    }
    return;
}
string s;
int main()
{
    ios::sync_with_stdio(false);
    int i,j,p,q;
    cin>>s>>k;
    memset(vis,false,sizeof(vis));
    memset(a,0,sizeof(a));
    memset(sum,0,sizeof(sum));
    for(i=0;i<10;i++)
        vis[i][i]=true;
    for(i=0;i<k;i++)
    {
        int ap,bp;
        cin>>ap>>bp;
        vis[ap][bp]=1;
    }
     for(p=0;p<10;p++)
     {
         for(i=0;i<=9;i++)
         {
             for(j=0;j<=9;j++)
             {
                 if(vis[i][p]&&vis[p][j])
                    vis[i][j]=1;
             }
         }
     }
     for(i=0;i<10;i++)
     {
         for(j=0;j<10;j++)
         {
            if(vis[i][j])
                a[i]++;
         }
     }
    sum[0]=1;
     for(i=0;i<s.length();i++)
     {
         cheng(a[s[i]-'0']);
     }
     for(i=100;i>=0;i--)
        if(sum[i]!=0)
        break;
     for(;i>=0;i--)
        cout<<sum[i];
  return 0;
}