結果

問題 No.5002 stick xor
ユーザー akakimidoriakakimidori
提出日時 2018-05-29 02:09:41
言語 C90
(gcc 11.4.0)
結果
AC  
実行時間 185 ms / 1,000 ms
コード長 6,205 bytes
コンパイル時間 5,873 ms
実行使用メモリ 956 KB
スコア 45,699
最終ジャッジ日時 2018-05-29 02:09:49
ジャッジサーバーID
（参考情報） 		judge8 /
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 145 ms
956 KB
testcase_01 AC 135 ms
956 KB
testcase_02 AC 149 ms
952 KB
testcase_03 AC 146 ms
952 KB
testcase_04 AC 133 ms
952 KB
testcase_05 AC 133 ms
956 KB
testcase_06 AC 130 ms
952 KB
testcase_07 AC 145 ms
952 KB
testcase_08 AC 137 ms
952 KB
testcase_09 AC 138 ms
956 KB
testcase_10 AC 133 ms
952 KB
testcase_11 AC 142 ms
952 KB
testcase_12 AC 150 ms
956 KB
testcase_13 AC 135 ms
956 KB
testcase_14 AC 185 ms
956 KB
testcase_15 AC 132 ms
952 KB
testcase_16 AC 147 ms
956 KB
testcase_17 AC 141 ms
956 KB
testcase_18 AC 143 ms
952 KB
testcase_19 AC 133 ms
956 KB
testcase_20 AC 129 ms
952 KB
testcase_21 AC 140 ms
956 KB
testcase_22 AC 140 ms
956 KB
testcase_23 AC 179 ms
952 KB
testcase_24 AC 141 ms
956 KB
testcase_25 AC 125 ms
956 KB
testcase_26 AC 152 ms
956 KB
testcase_27 AC 143 ms
956 KB
testcase_28 AC 147 ms
952 KB
testcase_29 AC 141 ms
956 KB
testcase_30 AC 143 ms
956 KB
testcase_31 AC 133 ms
952 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

typedef long long int int64;

#define MAX(a,b) ((a)>(b)?(a):(b))
#define MIN(a,b) ((a)<(b)?(a):(b))
#define ABS(a) ((a)>(0)?(a):-(a))

#define POS(i,j) ((i)*(n+1)+(j))

typedef struct length{
  int l;
  int index;
} len;

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

void flip(int l,char *s,int n,int *t){
  int a=t[0];
  int b=t[1];
  if(a==t[2]){
    int i;
    for(i=0;i<l;i++){
      int index=POS(a,b+i);
      s[index]=s[index]=='0'?'1':'0';
    }
  } else {
    int i;
    for(i=0;i<l;i++){
      int index=POS(a+i,b);
      s[index]=s[index]=='0'?'1':'0';
    }
  }
  return;
}

int func_put(int l,char *s,int n,int *t){
  int a=t[0];
  int b=t[1];
  if(a==t[2]){
    int cnt=0;
    int i;
    for(i=0;i<l;i++){
      int index=POS(a,b+i);
      cnt+=(s[index]=='0'?1:-1);
    }
    return cnt;
  } else {
    int cnt=0;
    int i;
    for(i=0;i<l;i++){
      int index=POS(a+i,b);
      cnt+=(s[index]=='0'?1:-1);
    }
    return cnt;
  }
}

int func(int l,char *s,const int n,int *a){
  int max=-l-1;
  int i,j;
  for(j=0;j<n;j++){
    int c=0;
    for(i=0;i<l;i++){
      c+=(s[POS(i,j)]=='0'?-1:1);
    }
    i=0;
    if(c>max){
      max=c;
      a[0]=i;
      a[1]=j;
      a[2]=i+l-1;
      a[3]=j;
    }
    while(i+l<n){
      c-=(s[POS(i,j)]=='0'?-1:1);
      c+=(s[POS(i+l,j)]=='0'?-1:1);
      i++;
      if(c>max){
	max=c;
	a[0]=i;
	a[1]=j;
	a[2]=i+l-1;
	a[3]=j;
      }
    }
  }
  for(i=0;i<n;i++){
    int c=0;
    for(j=0;j<l;j++){
      c+=(s[POS(i,j)]=='0'?-1:1);
    }
    j=0;
    if(c>max){
      max=c;
      a[0]=i;
      a[1]=j;
      a[2]=i;
      a[3]=j+l-1;
    }
    while(j+l<n){
      c-=(s[POS(i,j  )]=='0'?-1:1);
      c+=(s[POS(i,j+l)]=='0'?-1:1);
      j++;
      if(c>max){
	max=c;
	a[0]=i;
	a[1]=j;
	a[2]=i;
	a[3]=j+l-1;
      }
    }
  }
  flip(l,s,n,a);
  return max;
}

int score(const char *s,const int n){
  int res=0;
  int i,j;
  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      res+=(s[POS(i,j)]=='0'?1:0);
    }
  }
  return res;
}

void print(const char *s,const int n){
  printf("score : %d\n",score(s,n));
  int i,j;
  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      putchar(s[POS(i,j)]);
    }
    putchar('\n');
  }
  putchar('\n');
  return;
}

void swap(len *array,int x,int y){
  len t=array[x];
  array[x]=array[y];
  array[y]=t;
  return;
}

void shuffle(len *array,int len){
  int i;
  for(i=0;i<len;i++) swap(array,i,rand()%len);
  return;
}

int func_algo6(int l1,int l2,char *s,int n,int *a,int *b){
  int max=-l1-1;
  int i,j,k;
  for(i=0;i<n;i++){
    for(j=0;j+l1-1<n;j++){
      for(k=0;k+l2-1<l1;k++){
	int c=0;
	int x;
	for(x=j;x<j+k;x++) c+=(s[POS(i,x)]=='0'?-1:1);
	for(x=j+k+l2;x<j+l1;x++) c+=(s[POS(i,x)]=='0'?-1:1);
	if(c>max){
	  max=c;
	  a[0]=i;
	  a[1]=j;
	  a[2]=i;
	  a[3]=j+l1-1;
	  b[0]=i;
	  b[1]=j+k;
	  b[2]=i;
	  b[3]=j+k+l2-1;
	}
      }
      for(;k<l1 && j+k+l2-1<n;k++){
	int c=0;
	int x;
	for(x=j;x<j+k;x++) c+=(s[POS(i,x)]=='0'?-1:1);
	for(x=l1;x<k+l2;x++) c+=(s[POS(i,x)]=='0'?-1:1);
	if(c>max){
	  max=c;
	  a[0]=i;
	  a[1]=j;
	  a[2]=i;
	  a[3]=j+l1-1;
	  b[0]=i;
	  b[1]=j+k;
	  b[2]=i;
	  b[3]=j+k+l2-1;
	}
      }
    }
  }
  for(i=0;i+l1-1<n;i++){
    for(j=0;j<n;j++){
      for(k=0;k+l2-1<l1;k++){
	int c=0;
	int x;
	for(x=i;x<i+k;x++) c+=(s[POS(x,j)]=='0'?-1:1);
	for(x=i+k+l2;x<i+l1;x++) c+=(s[POS(x,j)]=='0'?-1:1);
	if(c>max){
	  max=c;
	  a[0]=i;
	  a[1]=j;
	  a[2]=i+l1-1;
	  a[3]=j;
	  b[0]=i+k;
	  b[1]=j;
	  b[2]=i+k+l2-1;
	  b[3]=j;
	}
      }
      for(;k<l1 && i+k+l2-1<n;k++){
	int c=0;
	int x;
	for(x=i;x<i+k;x++) c+=(s[POS(x,j)]=='0'?-1:1);
	for(x=l1;x<k+l2;x++) c+=(s[POS(x,j)]=='0'?-1:1);
	if(c>max){
	  max=c;
	  a[0]=i;
	  a[1]=j;
	  a[2]=i+l1-1;
	  a[3]=j;
	  b[0]=i+k;
	  b[1]=j;
	  b[2]=i+k+l2-1;
	  b[3]=j;
	}
      }
    }
  }
  flip(l1,s,n,a);
  flip(l2,s,n,b);
  return max;
}

void algo9(char *s,int n,len *array,int k,int *ans,int *tmp){
  const int parameter=4;
  const int use_parameter=5;
  int i;
  for(i=0;i<k && array[i].l>=parameter;i++);
  const int m=i;
  len *now=(len *)malloc(sizeof(len)*m);
  for(i=0;i<m;i++) now[i]=array[i];
  len *used=(len *)malloc(sizeof(len)*m);
  len *next=(len *)malloc(sizeof(len)*m);
  int nowLen=m;
  int flag=1;
  while(flag && nowLen>0){
    //print(s,n);
    flag=0;
    shuffle(now,nowLen);
    int usedLen=0;
    int nextLen=0;
    for(i=0;i<nowLen;i++){
      int t=func(now[i].l,s,n,ans+4*now[i].index);
      if(t<=now[i].l/use_parameter){
	flip(now[i].l,s,n,ans+4*now[i].index);
	next[nextLen++]=now[i];
      } else {
	flag=1;
	used[usedLen++]=now[i];
      }
    }
    for(i=0;i<usedLen;i++){
      int l=used[i].l;
      int index=used[i].index;
      if(func_put(l,s,n,ans+4*index)<=used[i].l/use_parameter){
	//printf("back\n");
	flip(l,s,n,ans+4*index);
	next[nextLen++]=used[i];
      }
    }
    len *swapL=now;
    now=next;
    next=swapL;
    nowLen=nextLen;
  }
  qsort(now,nowLen,sizeof(len),cmp);
  /*
  {
    int cnt[26]={0};
    for(i=0;i<nowLen;i++) cnt[now[i].l]++;
    for(i=25;i>=parameter;i--) printf("%d : %d\n",i,cnt[i]);

  }
  */
  const int f=4;
  if(nowLen%f!=0){
    for(i=nowLen/f*f;i<nowLen;i++){
      func(now[i].l,s,n,ans+4*now[i].index);
    }
  }
  nowLen=nowLen/f*f;
  int diff=nowLen/f;
  int j;
  for(j=0;j<f;j+=2){
    //print(s,n);
    for(i=j*diff;i<(j+1)*diff;i++){
      func_algo6(now[i].l,now[i+diff].l,s,n,ans+4*now[i].index,ans+4*now[i+diff].index);
    }
  }
  //print(s,n);
  for(i=m;i<k;i++){
    func(array[i].l,s,n,ans+4*array[i].index);
  }
  //print(s,n);
  return;
}

void run(void){
  int n,k;
  scanf("%d%d",&n,&k);
  char *s=(char *)calloc((n+1)*n,sizeof(char));
  len *array=(len *)malloc(sizeof(len)*k);
  int *ans=(int *)malloc(sizeof(int)*4*k);
  int *tmp=(int *)malloc(sizeof(int)*4*k);
  int i;
  for(i=0;i<k;i++){
    int t;
    scanf("%d",&t);
    array[i].l=t;
    array[i].index=i;
  }
  qsort(array,k,sizeof(len),cmp);
  for(i=0;i<n;i++) scanf("%s",s+POS(i,0));
  algo9(s,n,array,k,ans,tmp);
  for(i=0;i<k;i++){
    printf("%d %d %d %d\n",ans[i*4+0]+1,ans[i*4+1]+1,ans[i*4+2]+1,ans[i*4+3]+1);
  }
  return;
}

int main(int argc,char **argv){
  run();
  return 0;
}
0