CF889 E Mod Mod Mod——DP

题目：http://codeforces.com/contest/889/problem/E

这题真好玩。

官方题解说得很好。

想到相邻 a[ i ] 之间的段可能可以一起维护，但是不太会。

原来是表示成 i*x+k 的形式。其中 x 是具体的值，放在 DP 数组里只要记录 “ x<= ... 的 x 都满足这个式子” 就行，因为每次经过一个 i ，范围都是保留最底部的一些之类的。

然后写了一个不行的代码。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#define ll long long
using namespace std;
ll rdn()
{
  ll ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return fx?ret:-ret;
}
ll Mx(ll a,ll b){return a>b?a:b;}
const int N=2e5+,M=N*;
int n,tot,l[N],r[N];ll a[N],dp[M],dy[M],ans;
map<ll,int> mp[N];
int main()
{
  n=rdn();
  for(int i=;i<=n;i++)a[i]=rdn();
  r[]=++tot; dp[tot]=; dy[tot]=a[]-;
  for(int i=;i<n;i++)
    {
      l[i+]=r[i]=tot;
      for(int j=l[i]+,d;j<=r[i];j++)
    {
      ll k=dp[j],r=dy[j],tp=r%a[i+];
      if(r>=a[i+])
        {
          if(!mp[i+].count(a[i+]-))
        {
          mp[i+][a[i+]-]=++tot;
          dy[tot]=a[i+]-;
        }
          d=mp[i+][a[i+]-];
          dp[d]=Mx(dp[d],k+i*(r-tp-a[i+]));
        }
      if(!mp[i+].count(tp))
        {
          mp[i+][tp]=++tot; dy[tot]=tp;
        }
      d=mp[i+][tp];
      dp[d]=Mx(dp[d],k+i*(r-tp));
    }
    }
  for(int i=l[n]+;i<=tot;i++)
    {
      ll k=dp[i],r=dy[i];
      ans=Mx(ans,n*r+k);
    }
  printf("%lld\n",ans);
  return ;
}

不能遍历 r 没有变的位置。也不能让它们占用更多空间，比如在 i 处用一个位置、又在 i+1 处用一个位置。

然后抄了一番题解。原来 map 可以这样遍历。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#define ll long long
using namespace std;
ll rdn()
{
  ll ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return fx?ret:-ret;
}
ll Mx(ll a,ll b){return a>b?a:b;}
const int N=2e5+;
int n;ll a[N],ans;
map<ll,ll> mp;
map<ll,ll>::iterator it;
int main()
{
  n=rdn();for(int i=;i<=n;i++)a[i]=rdn();
  mp[a[]-]=;
  for(int i=;i<n;i++)
    {
      while()
    {
      it=mp.end(); it--;
      ll r=(*it).first, k=(*it).second;
      if(r<a[i+])break;
      mp.erase(it);
      mp[a[i+]-]=Mx(mp[a[i+]-],k+i*(r-r%a[i+]-a[i+]));
      mp[r%a[i+]]=Mx(mp[r%a[i+]],k+i*(r-r%a[i+]));
    }
    }
  for(it=mp.begin();it!=mp.end();it++)
    ans=Mx(ans,n*(*it).first+(*it).second);
  printf("%lld\n",ans);
  return ;
}