大意: 给定n*m初始字符矩阵, 有k个新矩阵, 第$i$个矩阵是由初始矩阵区间赋值得到的, 求选择一个新矩阵, 使得其余新矩阵到它距离和最小.

字符集比较小, 可以考虑每次区间覆盖对每个字符的贡献. 区间覆盖转化为差分, 然后前缀和优化.

刚开始辅助数组开多了, 卡内存卡了好久. 看别人代码似乎可以再优化掉两个辅助数组.

#include <iostream>
#include <sstream>
#include <algorithm>
#include <cstdio>
#include <math.h>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <string.h>
#include <bitset>
#define REP(i,a,n) for(int i=a;i<=n;++i)
#define PER(i,a,n) for(int i=n;i>=a;--i)
#define hr putchar(10)
#define pb push_back
#define lc (o<<1)
#define rc (lc|1)
#define mid ((l+r)>>1)
#define ls lc,l,mid
#define rs rc,mid+1,r
#define x first
#define y second
#define io std::ios::sync_with_stdio(false)
#define endl '\n'
#define DB(a) ({REP(__i,1,n) cout<<a[__i]<<' ';hr;})
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int P = 1e9+7, INF = 0x3f3f3f3f;
ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}
ll qpow(ll a,ll n) {ll r=1%P;for (a%=P;n;a=a*a%P,n>>=1)if(n&1)r=r*a%P;return r;}
ll inv(ll x){return x<=1?1:inv(P%x)*(P-P/x)%P;}
inline int rd() {int x=0;char p=getchar();while(p<'0'||p>'9')p=getchar();while(p>='0'&&p<='9')x=x*10+p-'0',p=getchar();return x;}
//head const int M = 3e5+10, N = 1e3+10;
int n, m, k;
char s[N][N];
int l[M], l2[M], r[M], r2[M], e[M];
int c[26][N][N], d[26][N][N];
ll sum[N][N], val[26][N][N];
void upd(int c[N][N], int l, int r, int l2, int r2, int v) {
c[l][r]+=v,c[l][r2+1]-=v,c[l2+1][r]-=v,c[l2+1][r2+1]+=v;
}
ll get(ll c[N][N], int id) {
return c[l2[id]][r2[id]]-c[l[id]-1][r2[id]]-c[l2[id]][r[id]-1]+c[l[id]-1][r[id]-1];
} int main() {
scanf("%d%d%d", &n, &m, &k);
REP(i,1,n) scanf("%s", s[i]+1);
REP(i,1,n) REP(j,1,m) s[i][j]-='a';
REP(i,1,k) {
char cc;
scanf("%d%d%d%d %c", l+i,r+i,l2+i,r2+i,&cc);
e[i] = cc-'a';
REP(t,0,25) upd(d[t],l[i],r[i],l2[i],r2[i],-1);
REP(t,0,25) upd(c[t],l[i],r[i],l2[i],r2[i],abs(e[i]-t));
}
REP(i,1,n) REP(j,1,m) {
REP(t,0,25) {
d[t][i][j] += d[t][i-1][j]+d[t][i][j-1]-d[t][i-1][j-1];
c[t][i][j] += c[t][i-1][j]+c[t][i][j-1]-c[t][i-1][j-1];
val[t][i][j] = (ll)(k+d[t][i][j])*abs(t-s[i][j])+c[t][i][j];
}
sum[i][j] = val[s[i][j]][i][j]+sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
REP(t,0,25) val[t][i][j]+=val[t][i-1][j]+val[t][i][j-1]-val[t][i-1][j-1];
}
ll ans = 1e18;
REP(i,1,k) ans=min(ans,sum[n][m]-get(sum,i)+get(val[e[i]],i));
printf("%lld\n", ans);
}

最新文章

  1. ASP.NET MVC 请求路径相关参数的获取
  2. JS心得——判断一个对象是否为空
  3. redis 集群热备自动切换sentinel配置实战
  4. 自定义RadioButton样式
  5. Nginx开发从入门到精通 学习目录分享学习 (阿里著作)
  6. C# Use Pop3Client to read gmail
  7. ThreadLocal类的实现用法
  8. hadoop多次格式化后,导致datanode启动不了,怎么办?(伪分布式)
  9. [Linked List]Remove Duplicates from Sorted List II
  10. Exec sql/c
  11. 【百度地图API】如何制作班级地理通讯录?LBS通讯录
  12. Listview性能优化
  13. jenkins gradle 编译遇到tomcat异常
  14. bzoj4198 荷马史诗 哈夫曼编码
  15. Windows环境下使用python 3.x自带的CGI服务器测试cgi脚本--Python
  16. mysql GTID
  17. 多轴APM飞控调参
  18. EntityFramework Code-First 简易教程(五)-------领域类配置
  19. Hive记录-配置支持事务管理hive-site.xml
  20. [troubleshoot][automake] automake编译的时候发生死循环

热门文章

  1. RTS打卡计划第四周
  2. java判断一个字符串是否为数字型
  3. CI集成Smarty的实现方式
  4. linu逻辑分区动态调整大小
  5. ubuntu关于ssh协议登录问题
  6. Groovy脚本基础全攻略
  7. pycloudtag_wordcloud 中英文标签云
  8. 关于Win7 64位下:Navicat无法连接64位Oracle 11gR2:Cannot load OCI DLL解决方法
  9. 基于 Open vSwitch 的 OpenFlow 实践
  10. oracle增、删、改、查