[暴力DP][暴力STL]SNOI2017 && LOJ#2256. 英雄联盟

列出DP方程
fi,j=min{fi1,jk+k×Cik|jkKi }
fi,j 表示前 i 个物品产生 j 中方案时的最小价格

j 很大,有 1018 ,但是 j 只有2,3,5,7三种质数,所以 j 能取到的值上限为 (logm)4 ,而且远小于这个上限。
因为懒,所以不继续优化下去了……DP数组直接用map存一下,转移时遍历一遍。

然而本地跑了10多秒,开了O2跑进3秒,交上去T了最后一个点…

map不行就换unordered_map,硬生生跑进0.5s……

#include 
#include 
#include 
#include 
#define fi first
#define se second

using namespace std;
using namespace tr1;

typedef long long ll;

unordered_map f,g;

int n,k[210],c[210];
ll m,ans;

inline void fix(ll &x,ll y){
  if(x>y) x=y;
}

int main(){
  scanf("%d%lld\n",&n,&m);
  for(int i=1;i<=n;i++) scanf("%d",&k[i]);
  for(int i=1;i<=n;i++) scanf("%d",&c[i]);
  f[0]=0; ans=1LL<<60;
  for(int i=1;i<=n;i++){
    for(auto v : f){
      if(v.fi==0){
    for(int j=1;j<=k[i];j++)
      if(g.count(j)) fix(g[j],j*c[i]);
      else g[j]=j*c[i];
    g[0]=0;
    continue;
      }
      if(g.count(v.fi)) fix(g[v.fi],v.se);
      else g[v.fi]=v.se;
      for(int j=2;j<=k[i];j++){
    if(v.fi*j>=m){
      fix(ans,v.se+j*c[i]);
      break;
    }
    if(v.se+j*c[i]>=ans) break;
    if(g.count(v.fi*j)) fix(g[v.fi*j],v.se+c[i]*j);
    else g[v.fi*j]=v.se+c[i]*j;
      }
    }
    swap(f,g); g.clear();
  }
  printf("%lld\n",ans);
  return 0;
}

你可能感兴趣的:(乱搞,DP,STL)