結果
| 問題 |
No.5002 stick xor
|
| コンテスト | |
| ユーザー |
 akakimidori
|
| 提出日時 | 2018-05-29 02:09:41 |
| 言語 | C90 (gcc 12.3.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 / |
| 純コード判定しない問題か言語 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| other | AC * 32 |
ソースコード
#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;
}