結果

問題 No.738 平らな農地
ユーザー akakimidori
提出日時 2018-12-04 03:57:16
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 64 ms / 2,000 ms
コード長 2,068 bytes
コンパイル時間 314 ms
コンパイル使用メモリ 31,872 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-05 21:13:47
合計ジャッジ時間 5,730 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 87
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<stdio.h>
#include<stdlib.h>

typedef long long int int64;

int cmp(const void *a,const void *b){
  return *(int *)a-*(int *)b;
}

int compress(int *a,const int n){
  qsort(a,n,sizeof(int),cmp);
  int i=0;
  int j=0;
  while(i<n){
    a[j++]=a[i++];
    while(i<n && cmp(a+j-1,a+i)==0) i++;
  }
  return j;
}

int toIndex(const int *za,const int len,const int v){
  int l=0;
  int r=len;
  while(r-l>1){
    int m=(l+r)/2;
    if(cmp(za+m,&v)<=0){
      l=m;
    } else {
      r=m;
    }
  }
  return l;
}

void add(int64 *bit,int index,int64 v){
  int n=bit[0];
  for(int i=index;i<=n;i+=i&-i) bit[i]+=v;
}

int64 sum(int64 *bit,int index){
  int64 res=0;
  for(int i=index;i>0;i-=i&-i) res+=bit[i];
  return res;
}

int search(int64 *bit,int64 x){
  const int n=bit[0];
  int k=1;
  while(2*k<n) k*=2;
  int y=0;
  for(;k>0;k>>=1){
    if(y+k<=n && bit[y+k]<x){
      x-=bit[y+k];
      y+=k;
    }
  }
  return y+1;
}

int64 calc(int64 *cntBit,int64 *sumBit,int *za,int k){
  int med=(k+1)/2;
  int index=search(cntBit,med);
  int64 v=za[index-1];
  int64 s=sum(sumBit,sumBit[0]);
  int64 low=sum(sumBit,index)-v*(sum(cntBit,index)-med);
  int64 res=(med*v-low)+(s-low-v*(k-med));
  return res;
}

void run(void){
  int n,k;
  scanf("%d%d",&n,&k);
  int *a=(int *)calloc(n,sizeof(int));
  int i;
  for(i=0;i<n;i++) scanf("%d",a+i);
  int *za=(int *)calloc(n,sizeof(int));
  for(i=0;i<n;i++) za[i]=a[i];
  int len=compress(za,n);
  int64 *cntBit=(int64 *)calloc(len+1,sizeof(int64));
  int64 *sumBit=(int64 *)calloc(len+1,sizeof(int64));
  cntBit[0]=sumBit[0]=len;
  for(i=0;i<k;i++){
    int index=toIndex(za,len,a[i]);
    add(cntBit,index+1,1);
    add(sumBit,index+1,a[i]);
  }
  int64 min=calc(cntBit,sumBit,za,k);
  for(i=k;i<n;i++){
    int index=toIndex(za,len,a[i]);
    add(cntBit,index+1,1);
    add(sumBit,index+1,a[i]);
    index=toIndex(za,len,a[i-k]);
    add(cntBit,index+1,-1);
    add(sumBit,index+1,-a[i-k]);
    int64 t=calc(cntBit,sumBit,za,k);
    if(t<min) min=t;
  }
  printf("%lld\n",min);
}

int main(void){
  run();
  return 0;
}
0