http://www.lydsy.com/JudgeOnline/problem.php?id=4241 (题目链接)

  看到题目就联想到了【bzoj2809】 Apio2012—dispatching。想了想权值分块+莫队,发现不好维护块内最值,又看了看80s的时间,于是怒水一发线段树+莫队,结果先WA后TLE,不断TLE,无论怎么改常数都不行,难道nlogn*sqrt(n)就是过不了吗!!不爽,蒯个题解,再见!

题意

  求区间加权众数。

solution

  貌似是分块,离散化之后,用mx[i][j]表示第i块到第j块的答案,cnt[i][j]表示前i块数字j的个数。

代码

// bzoj4241
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf 2147483640
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;
inline LL getint() {
int f,x=0;char ch=getchar();
while (ch<='0' || ch>'9') {if (ch=='-') f=-1;else f=1;ch=getchar();}
while (ch>='0' && ch<='9') {x=x*10+ch-'0';ch=getchar();}
return x*f;
} const int maxn=100010;
struct data {int w,id;}a[maxn];
int b[maxn],pos[maxn],L[400],R[400],cnt[400][maxn],num[maxn],st[maxn];
LL mx[400][400],ans;
int n,m,block,tot,top; bool cmp(data a,data b) {
return a.w<b.w;
}
LL cal(int l,int r) {
LL ans=0;
top=0;
for (int i=l;i<=r;i++) {
if (!num[b[i]]) st[++top]=b[i];
num[b[i]]++;
ans=max(ans,(LL)num[b[i]]*a[b[i]].w);
}
while (top) num[st[top--]]=0;
return ans;
}
void build() {
block=sqrt(n);tot=(n-1)/block+1;
for (int i=1;i<=n;i++) {
pos[i]=(i-1)/block+1;
if (!L[pos[i]]) L[pos[i]]=i;
R[pos[i]]=i;
}
for (int i=1;i<=n;i++) cnt[pos[i]][b[i]]++;
for (int i=2;i<=tot;i++)
for (int j=1;j<=n;j++) cnt[i][j]+=cnt[i-1][j];
}
int main() {
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++) scanf("%d",&a[i].w),a[i].id=i;
sort(a+1,a+1+n,cmp);
b[a[1].id]=1;
for (int i=2;i<=n;i++) {
if (a[i].w==a[i-1].w) b[a[i].id]=b[a[i-1].id];
else b[a[i].id]=i;
}
build();
for (int i=1;i<=tot;i++) {
LL ans=0;
memset(num,0,sizeof(num));
for (int j=L[i];j<=n;j++) {
num[b[j]]++;
ans=max(ans,(LL)num[b[j]]*a[b[j]].w);
if (j==R[pos[j]]) mx[i][pos[j]]=ans;
}
}
memset(num,0,sizeof(num));
for (int i=1;i<=m;i++) {
int l,r;
scanf("%d%d",&l,&r);
ans=0;
if (pos[l]==pos[r]) printf("%lld\n",cal(l,r));
else {
ans=mx[pos[l]+1][pos[r]-1];
top=0;
for (int i=l;i<=R[pos[l]];i++)
if (!num[b[i]]) st[++top]=b[i],num[b[i]]=cnt[pos[r]-1][b[i]]-cnt[pos[l]][b[i]];
for (int i=L[pos[r]];i<=r;i++)
if (!num[b[i]]) st[++top]=b[i],num[b[i]]=cnt[pos[r]-1][b[i]]-cnt[pos[l]][b[i]];
for (int i=l;i<=R[pos[l]];i++) num[b[i]]++;
for (int i=L[pos[r]];i<=r;i++) num[b[i]]++;
while (top) {
ans=max(ans,(LL)a[st[top]].w*num[st[top]]);
num[st[top--]]=0;
}
printf("%lld\n",ans);
}
}
return 0;
}

  

最新文章

  1. 创建 Pool &amp; VIP - 每天5分钟玩转 OpenStack(122)
  2. JS this指向
  3. Masonry使用注意事项
  4. Webform(Linq增删改查)
  5. web.config中httpRunTime的属性
  6. ruby on rails 里使用SideKiq 做后台任务
  7. 收藏:关于UseSubmitBehavior和OnClientClick同时使用,导致无法触发后台事件的问题
  8. 集成友盟分享SDK报错
  9. .net web部署(IIS Express &amp;&amp; Nancy Self-Hosting)
  10. Modernizr.js介绍与使用
  11. final的深入理解 - final数据
  12. MySQL数据库中文变问号
  13. (十六)java中的String
  14. 原生js写轮播图效果
  15. Caused by: Java.lang.NoSuchMethodError: javax.persistence.JoinColumn.foreignKey()Ljavax/persistence/ForeignKey;
  16. go 多维度 Map 的数据存取
  17. (转载)windows下安装配置Xampp
  18. STDIN_FILENO的作用及与stdin 的区别
  19. Mybatis系列(一):Mybatis入门
  20. HDU 4054 Number String

热门文章

  1. css 内容超过容器宽度,checkbox等控件不会随着内容延伸
  2. mybatis 3.x 缓存Cache的使用
  3. 字符串匹配(hash算法)
  4. centos设置静态IP
  5. ASP.NET MVC3入门教程之第一个WEB应用程序
  6. 文本 To 音频
  7. C语言字符串与数字相互转换
  8. 【SQL语句】 - Ctrl+3 查询表属性的存储过程
  9. MVC 中的 ispostback
  10. iOS 自定义NavigationBar右侧按钮rightBarButtonItem--button