#301288 (C) No.738 平らな農地

提出ソース

結果

問題	No.738 平らな農地
ユーザー	akakimidori
提出日時	2018-12-04 03:57:16
言語	C (gcc 12.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

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 5,248 KB
testcase_01	AC	1 ms 5,376 KB
testcase_02	AC	1 ms 5,376 KB
testcase_03	AC	1 ms 5,376 KB
testcase_04	AC	1 ms 5,376 KB
testcase_05	AC	2 ms 5,376 KB
testcase_06	AC	2 ms 5,376 KB
testcase_07	AC	1 ms 5,376 KB
testcase_08	AC	1 ms 5,376 KB
testcase_09	AC	1 ms 5,376 KB
testcase_10	AC	1 ms 5,376 KB
testcase_11	AC	1 ms 5,376 KB
testcase_12	AC	2 ms 5,376 KB
testcase_13	AC	1 ms 5,376 KB
testcase_14	AC	1 ms 5,376 KB
testcase_15	AC	42 ms 5,376 KB
testcase_16	AC	43 ms 5,376 KB
testcase_17	AC	46 ms 5,376 KB
testcase_18	AC	44 ms 5,376 KB
testcase_19	AC	52 ms 5,376 KB
testcase_20	AC	43 ms 5,376 KB
testcase_21	AC	52 ms 5,376 KB
testcase_22	AC	44 ms 5,376 KB
testcase_23	AC	48 ms 5,376 KB
testcase_24	AC	49 ms 5,376 KB
testcase_25	AC	1 ms 5,376 KB
testcase_26	AC	1 ms 5,376 KB
testcase_27	AC	1 ms 5,376 KB
testcase_28	AC	1 ms 5,376 KB
testcase_29	AC	1 ms 5,376 KB
testcase_30	AC	1 ms 5,376 KB
testcase_31	AC	1 ms 5,376 KB
testcase_32	AC	1 ms 5,376 KB
testcase_33	AC	1 ms 5,376 KB
testcase_34	AC	2 ms 5,376 KB
testcase_35	AC	1 ms 5,376 KB
testcase_36	AC	1 ms 5,376 KB
testcase_37	AC	1 ms 5,376 KB
testcase_38	AC	1 ms 5,376 KB
testcase_39	AC	1 ms 5,376 KB
testcase_40	AC	1 ms 5,376 KB
testcase_41	AC	1 ms 5,376 KB
testcase_42	AC	1 ms 5,376 KB
testcase_43	AC	1 ms 5,376 KB
testcase_44	AC	1 ms 5,376 KB
testcase_45	AC	46 ms 5,376 KB
testcase_46	AC	42 ms 5,376 KB
testcase_47	AC	44 ms 5,376 KB
testcase_48	AC	40 ms 5,376 KB
testcase_49	AC	39 ms 5,376 KB
testcase_50	AC	41 ms 5,376 KB
testcase_51	AC	47 ms 5,376 KB
testcase_52	AC	42 ms 5,376 KB
testcase_53	AC	45 ms 5,376 KB
testcase_54	AC	45 ms 5,376 KB
testcase_55	AC	46 ms 5,376 KB
testcase_56	AC	46 ms 5,376 KB
testcase_57	AC	41 ms 5,376 KB
testcase_58	AC	44 ms 5,376 KB
testcase_59	AC	45 ms 5,376 KB
testcase_60	AC	44 ms 5,376 KB
testcase_61	AC	44 ms 5,376 KB
testcase_62	AC	40 ms 5,376 KB
testcase_63	AC	47 ms 5,376 KB
testcase_64	AC	45 ms 5,376 KB
testcase_65	AC	54 ms 5,376 KB
testcase_66	AC	57 ms 5,376 KB
testcase_67	AC	34 ms 5,376 KB
testcase_68	AC	33 ms 5,376 KB
testcase_69	AC	41 ms 5,376 KB
testcase_70	AC	38 ms 5,376 KB
testcase_71	AC	12 ms 5,376 KB
testcase_72	AC	15 ms 5,376 KB
testcase_73	AC	15 ms 5,376 KB
testcase_74	AC	15 ms 5,376 KB
testcase_75	AC	40 ms 5,376 KB
testcase_76	AC	35 ms 5,376 KB
testcase_77	AC	41 ms 5,376 KB
testcase_78	AC	45 ms 5,376 KB
testcase_79	AC	45 ms 5,376 KB
testcase_80	AC	38 ms 5,376 KB
testcase_81	AC	41 ms 5,376 KB
testcase_82	AC	43 ms 5,376 KB
testcase_83	AC	15 ms 5,376 KB
testcase_84	AC	18 ms 5,376 KB
testcase_85	AC	64 ms 5,376 KB
testcase_86	AC	58 ms 5,376 KB
testcase_87	AC	37 ms 5,376 KB
testcase_88	AC	34 ms 5,376 KB
testcase_89	AC	0 ms 5,376 KB
testcase_90	AC	1 ms 5,376 KB
testcase_91	AC	1 ms 5,376 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#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;
}

yukicoder

結果

テストケース

ソースコード