有一個二維陣列a[5000][5000];
a[0]={0,0,0,1,0,0,0,0,1,1,*******0,1,1,0};
a[1]={0,0,1,1,0,1,0,0,0,0,*******1,0,1,0};
*********
a[4998]={1,0,0,1,0,1,0,1,0,0,*******0,0,0,0};
a[4999]={0,1,1,1,0,1,0,0,0,0,*******0,0,1,1};
元素的排列是隨機的,a[*]中1的個數,位置也是隨機的;
需要找到一個b[8],使得符合((((a[*][b[0]]=1)&&(a[*][b[1]]=1))&&((a[*][b[2]]=1)&&(a[*][b[3]]=1)))&&(((a[*][b[4]]=1)&&(a[*][b[5]]=1))&&((a[*][b[6]]=1)&&(a[*][b[7]]=1))))這個條件的a[*]的個數是最多的,這樣的演算法應該是怎么樣的?望賜教!
使用過兩個FOR回圈后+加5000中選8個的組合演算法,但是消耗時間太久了,有沒有什么快速的辦法?
uj5u.com熱心網友回復:
計算量不算太大啊,貼出你的代碼。Talk is cheap; show me the code.
uj5u.com熱心網友回復:
main()
{
int a[5000][5000]={{0,1,0,0,0,0,1,1,1,1,0,***},{1,0,1,0,1,0,1,0,1,1,1,0,1,0,***},{0,0,0,0,1,1,1,1,1,0,0,***},****};
int v[8],bestv[8],b,fitnum,maxfitnum;
int i;
maxfitnum;
for(v[0]=0;v[0]<4993;v[0]++)
{for(v[1]=0;v[1]<4994;v[1]++)
{for(v[2]=0;v[2]<4995;v[2]++)
{for(v[3]=0;v[3]<4996;v[3]++)
{for(v[4]=0;v[4]<4997;v[4]++)
{for(v[5]=0;v[5]<4998;v[5]++)
{for(v[6]=0;v[6]<4999;v[6]++)
{for(v[7]=0;v[7]<5000;v[7]++)
{fitnum=0;
for(i=0;i<5000;i++)
{b=a[i][v[0]]+a[i][v[1]]+a[i][v[2]]+a[i][v[3]]+a[i][v[4]]+a[i][v[5]]+a[i][v[6]]+a[i][v[7]];
if(b==8)fitnum++;
};
if(fitnum>maxfitnum)
{bestv[0]=v[0];bestv[1]=v[1];bestv[2]=v[2];bestv[3]=v[3];bestv[4]=v[4];bestv[5]=v[5];bestv[6]=v[6];bestv[7]=v[7];
};
};
};
};
};
};
};
};
};
}
uj5u.com熱心網友回復:
main()
{
int a[5000][5000]={{0,1,0,0,0,0,1,1,1,1,0,***},{1,0,1,0,1,0,1,0,1,1,1,0,1,0,***},{0,0,0,0,1,1,1,1,1,0,0,***},****};
int v[8],bestv[8],b,fitnum,maxfitnum;
int i;
maxfitnum=0;
for(v[0]=0;v[0]<4993;v[0]++)
{for(v[1]=0;v[1]<4994;v[1]++)
{for(v[2]=0;v[2]<4995;v[2]++)
{for(v[3]=0;v[3]<4996;v[3]++)
{for(v[4]=0;v[4]<4997;v[4]++)
{for(v[5]=0;v[5]<4998;v[5]++)
{for(v[6]=0;v[6]<4999;v[6]++)
{for(v[7]=0;v[7]<5000;v[7]++)
{fitnum=0;
for(i=0;i<5000;i++)
{b=a[i][v[0]]+a[i][v[1]]+a[i][v[2]]+a[i][v[3]]+a[i][v[4]]+a[i][v[5]]+a[i][v[6]]+a[i][v[7]];
if(b==8)fitnum++;
};
if(fitnum>maxfitnum)
{bestv[0]=v[0];bestv[1]=v[1];bestv[2]=v[2];bestv[3]=v[3];bestv[4]=v[4];bestv[5]=v[5];bestv[6]=v[6];bestv[7]=v[7];
};
};
};
};
};
};
};
};
};
}
uj5u.com熱心網友回復:
main()
{
int a[5000][5000]={{0,1,0,0,0,0,1,1,1,1,0,***},{1,0,1,0,1,0,1,0,1,1,1,0,1,0,***},{0,0,0,0,1,1,1,1,1,0,0,***},****};
int v[8],bestv[8],b,fitnum,maxfitnum;
int i;
maxfitnum=0;
for(v[0]=0;v[0]<4993;v[0]++)
{for(v[1]=v[0]+1;v[1]<4994;v[1]++)
{for(v[2]=v[1]+1;v[2]<4995;v[2]++)
{for(v[3]=v[2]+1;v[3]<4996;v[3]++)
{for(v[4]=v[3]+1;v[4]<4997;v[4]++)
{for(v[5]=v[4]+1;v[5]<4998;v[5]++)
{for(v[6]=v[5]+1;v[6]<4999;v[6]++)
{for(v[7]=v[6]+1;v[7]<5000;v[7]++)
{fitnum=0;
for(i=0;i<5000;i++)
{b=a[i][v[0]]+a[i][v[1]]+a[i][v[2]]+a[i][v[3]]+a[i][v[4]]+a[i][v[5]]+a[i][v[6]]+a[i][v[7]];
if(b==8)fitnum++;
};
if(fitnum>maxfitnum)
{bestv[0]=v[0];bestv[1]=v[1];bestv[2]=v[2];bestv[3]=v[3];bestv[4]=v[4];bestv[5]=v[5];bestv[6]=v[6];bestv[7]=v[7];
};
};
};
};
};
};
};
};
};
}
uj5u.com熱心網友回復:
for(v[0]=0;v[0]<4993;v[0]++)
{for(v[1]=v[0]+1;v[1]<4994;v[1]++)
{for(v[2]=v[1]+1;v[2]<4995;v[2]++)
{for(v[3]=v[2]+1;v[3]<4996;v[3]++)
{for(v[4]=v[3]+1;v[4]<4997;v[4]++)
{for(v[5]=v[4]+1;v[5]<4998;v[5]++)
{for(v[6]=v[5]+1;v[6]<4999;v[6]++)
{for(v[7]=v[6]+1;v[7]<5000;v[7]++)
uj5u.com熱心網友回復:
這個是比較耗時間,想知道有什么比較快速的辦法?
uj5u.com熱心網友回復:
int first1(int p) {
int min1=5000,j;
for (int i=0;i<5000;i++) {
for (j=p;j<5000;j++) {
if (a[i][j]==1) break;
}
if (j<min1) min1=j;
}
return min1;
}
void next1(int &p) {
int min1=5000,j;
for (int i=0;i<5000;i++) {
for (j=p+1;j<5000;j++) {
if (a[i][j]==1) break;
}
if (j<min1) min1=j;
}
p=min1;
}
for(v[0]=first1( 0);v[0]<4993;next1(v[0])) {
for(v[1]=first1(v[0]+1);v[1]<4994;next1(v[1])) {
for(v[2]=first1(v[1]+1);v[2]<4995;next1(v[2])) {
for(v[3]=first1(v[2]+1);v[3]<4996;next1(v[3])) {
for(v[4]=first1(v[3]+1);v[4]<4997;next1(v[4])) {
for(v[5]=first1(v[4]+1);v[5]<4998;next1(v[5])) {
for(v[6]=first1(v[5]+1);v[6]<4999;next1(v[6])) {
for(v[7]=first1(v[6]+1);v[7]<5000;next1(v[7])) {
uj5u.com熱心網友回復:
資料中只有0和1?uj5u.com熱心網友回復:
還有更好的辦法嗎?這個計算量是減少了很多,但還是不少,C(5000,8)=5000*4999*4998*4997*4996*4995*4994*4993/8/7/6/5/4/3/2/1=9,633,991,198,641,591,860,780,625,就散減少了10倍,還有963399119864159186078062,是否還有更好的演算法uj5u.com熱心網友回復:
資料中只有0和1uj5u.com熱心網友回復:
演算法需要改改,比如改為下面這樣:i從0到5000-1回圈
j從0到5000-8
從a[i][j]到a[i][5000-8]試圖找到8個出現1的位置保存到w[0..7]中,得到一種特征v[0..7]=w[0..7]-j
統計這種特征在整個a中出現的次數,并記錄最多次數
uj5u.com熱心網友回復:
演算法需要改改,比如改為下面這樣:i從0到5000-1回圈
j從0到5000-8
從a[i][j]到a[i][5000-8]試圖找到8個出現1的位置保存到w[0..7]中,得到一種特征
v[0]=w[0];
v[1]=w[1]-w[0];
v[2]=w[2]-w[0];
v[3]=w[3]-w[0];
v[4]=w[4]-w[0];
v[5]=w[5]-w[0];
v[6]=w[6]-w[0];
v[7]=w[7]-w[0];
統計這種特征在整個a中出現的次數,并記錄最多次數
uj5u.com熱心網友回復:
演算法需要改改,比如改為下面這樣:i從0到5000-1回圈
j從0到5000-8
從a[i][j]到a[i][5000-8]試圖找到8個出現1的位置保存到w[0..7]中,得到一種特征
v[0]=0;
v[1]=w[1]-w[0];
v[2]=w[2]-w[0];
v[3]=w[3]-w[0];
v[4]=w[4]-w[0];
v[5]=w[5]-w[0];
v[6]=w[6]-w[0];
v[7]=w[7]-w[0];
統計這種特征在整個a中出現的次數,并記錄最多次數
uj5u.com熱心網友回復:
沒明白為什么要從a[i][j]到a[i][5000-8]試圖找到8個出現1的位置保存到w[0..7]中,得到一種特征v[0]=w[0];
v[1]=w[1]-w[0];
v[2]=w[2]-w[0];
v[3]=w[3]-w[0];
v[4]=w[4]-w[0];
v[5]=w[5]-w[0];
v[6]=w[6]-w[0];
v[7]=w[7]-w[0];
uj5u.com熱心網友回復:
沒明白為什么要從a[i][j]到a[i][5000-8]試圖找到8個出現1的位置保存到w[0..7]中,得到一種特征v[0]=0;
v[1]=w[1]-w[0];
v[2]=w[2]-w[0];
v[3]=w[3]-w[0];
v[4]=w[4]-w[0];
v[5]=w[5]-w[0];
v[6]=w[6]-w[0];
v[7]=w[7]-w[0];
uj5u.com熱心網友回復:
優化第一步,把1的數量不足8個的那些行標記了,不參與測驗uj5u.com熱心網友回復:
一定都超過8個uj5u.com熱心網友回復:
僅供參考:(測驗有bug)#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define N 300 //5000
static char a[N][N];
int v[8],w[8];
unsigned long ulrand(void) {
return (
(((unsigned long)rand()<<24)&0xFF000000ul)
|(((unsigned long)rand()<<12)&0x00FFF000ul)
|(((unsigned long)rand() )&0x00000FFFul));
}
int first1(int i,int p) {
int j;
for (j=p;j<N;j++) {
if (a[i][j]==1) break;
}
return j;
}
void next1(int i,int &p) {
int j;
for (j=p+1;j<N;j++) {
if (a[i][j]==1) break;
}
p=j;
}
bool findpattern(int i,int &p) {
static int y=-1;
int k,q;
bool r=false;
if (y!=i) {
y=i;
w[0]=first1(y,0);
if (w[0]<N-8) {
q=w[0];
for (k=1;k<8;k++) {
next1(y,q);
if (q<N) {
w[k]=q;
} else break;
}
if (k>=8) {
v[0]=0;
for (k=1;k<8;k++) v[k]=w[k]-w[0];
p=w[0];
return true;
}
}
} else {
for (k=0;k<8-1;k++) w[k]=w[k+1];
q=w[7]+1;
next1(y,q);
if (q<N) {
w[7]=q;
v[0]=0;
for (k=1;k<8;k++) v[k]=w[k]-w[0];
p=w[0];
return true;
}
}
return false;
}
int main() {
int bestv[8],b,fitnum,maxfitnum;
int i,j,y,p;
srand(time(NULL));
// srand(0);
for (i=0;i<N;i++) {
for (int j=0;j<N;j++) {
if (ulrand()%(N*N/1000)==0) a[i][j]=1;
// printf("%d",a[i][j]);
}
// printf("\n");
}
maxfitnum=0;
for (y=0;y<N;y++) {
p=0;
while (1) {
if (findpattern(y,p)) {
fitnum=0;
for (i=0;i<N;i++) {
for (j=first1(i,0);j<N-v[7];next1(i,j)) {
b=a[i][v[0]+j]+a[i][v[1]+j]+a[i][v[2]+j]+a[i][v[3]+j]+a[i][v[4]+j]+a[i][v[5]+j]+a[i][v[6]+j]+a[i][v[7]+j];
if (b==8) fitnum++;
}
}
if (fitnum>maxfitnum) {
bestv[0]=v[0];bestv[1]=v[1];bestv[2]=v[2];bestv[3]=v[3];bestv[4]=v[4];bestv[5]=v[5];bestv[6]=v[6];bestv[7]=v[7];
}
p++;
} else break;
}
}
printf("\nbestv:");
for (i=0;i<8;i++) {
printf("%6d ",bestv[i]);
}
return 0;
}
uj5u.com熱心網友回復:
編譯報錯:
uj5u.com熱心網友回復:
改了幾個地方,但還是有bug:#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define N 5000
static char a[N][N];
int v[8],w[8];
unsigned long ulrand(void) {
return (
(((unsigned long)rand()<<24)&0xFF000000ul)
|(((unsigned long)rand()<<12)&0x00FFF000ul)
|(((unsigned long)rand() )&0x00000FFFul));
}
int first1(int i) {
int j;
for (j=0;j<N;j++) {
if (a[i][j]==1) break;
}
return j;
}
void next1(int i,int &p) {
int j;
if (p>=N) {p=N; return;}
for (j=p+1;j<N;j++) {
if (a[i][j]==1) break;
}
p=j;
}
bool findpattern(int i,int &p) {
static int y=-1;
int k,q;
if (p>=N) return false;
if (y!=i) {
y=i;
w[0]=first1(y);
if (w[0]<N-8) {
q=w[0];
for (k=1;k<8;k++) {
next1(y,q);
if (q<N) {
w[k]=q;
} else break;
}
if (k>=8) {
v[0]=0;
for (k=1;k<8;k++) v[k]=w[k]-w[0];
p=w[0];
return true;
}
}
} else {
for (k=0;k<8-1;k++) w[k]=w[k+1];
q=w[7]+1;
next1(y,q);
if (q<N) {
w[7]=q;
v[0]=0;
for (k=1;k<8;k++) v[k]=w[k]-w[0];
p=w[0];
return true;
}
}
return false;
}
int main() {
int bestv[8],b,fitnum,maxfitnum;
int i,j,y,p;
srand(time(NULL));
// srand(0);
for (i=0;i<N;i++) {
for (int j=0;j<N;j++) {
if (ulrand()%(N*N/1000)==0) a[i][j]=1;
// printf("%d",a[i][j]);
}
// printf("\n");
}
maxfitnum=0;
for (y=0;y<N;y++) {
p=0;
while (1) {
if (findpattern(y,p)) {
fitnum=0;
for (i=0;i<N;i++) {
for (j=first1(i);j<N-v[7];next1(i,j)) {
b=a[i][v[0]+j]+a[i][v[1]+j]+a[i][v[2]+j]+a[i][v[3]+j]+a[i][v[4]+j]+a[i][v[5]+j]+a[i][v[6]+j]+a[i][v[7]+j];
if (b==8) fitnum++;
}
}
if (fitnum>maxfitnum) {
bestv[0]=v[0];bestv[1]=v[1];bestv[2]=v[2];bestv[3]=v[3];bestv[4]=v[4];bestv[5]=v[5];bestv[6]=v[6];bestv[7]=v[7];
}
p++;
if (p>N-8) break;
} else break;
}
}
printf("\nbestv:");
for (i=0;i<8;i++) {
printf("%6d ",bestv[i]);
}
return 0;
}
uj5u.com熱心網友回復:
變數個數一定是a[5000][5000]這樣大,小了就不用發帖了!uj5u.com熱心網友回復:
希望有個沒有BUG的uj5u.com熱心網友回復:
基本沒bug了:#include <conio.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <windows.h>
#define N 5000
static char a[N][N];
int v[8],w[8];
unsigned long ulrand(void) {
return (
(((unsigned long)rand()<<24)&0xFF000000ul)
|(((unsigned long)rand()<<12)&0x00FFF000ul)
|(((unsigned long)rand() )&0x00000FFFul));
}
int first1(int i) {
int j;
for (j=0;j<N;j++) {
if (a[i][j]==1) break;
}
return j;
}
void next1(int i,int &p) {
int j;
if (p>=N) {p=N; return;}
for (j=p+1;j<N;j++) {
if (a[i][j]==1) break;
}
p=j;
}
bool findpattern(int i,int &p) {
static int y=-1;
int k,q;
if (p>=N) return false;
if (y!=i) {
y=i;
w[0]=first1(y);
if (w[0]<N-8) {
q=w[0];
for (k=1;k<8;k++) {
next1(y,q);
if (q<N) {
w[k]=q;
} else break;
}
if (k>=8) {
v[0]=0;
for (k=1;k<8;k++) {
v[k]=w[k]-w[0];
if (v[k]>=N) {
DebugBreak();
}
}
p=w[0];
return true;
}
}
} else {
for (k=0;k<8-1;k++) w[k]=w[k+1];
q=w[7]+1;
next1(y,q);
if (q<N) {
w[7]=q;
v[0]=0;
for (k=1;k<8;k++) {
v[k]=w[k]-w[0];
if (v[k]>=N) {
DebugBreak();
}
}
p=w[0];
return true;
}
}
return false;
}
int main() {
int bestv[8],b,fitnum,maxfitnum;
int i,j,y,p;
// srand(time(NULL));
srand(0);
for (i=0;i<N;i++) {
for (int j=0;j<N;j++) {
if ((ulrand()%N) < (N*1/100)) a[i][j]=1;
// printf("%d",a[i][j]);
}
// printf("\n");
}
maxfitnum=0;
for (y=0;y<N;y++) {
printf("%04d/%04d\n",y+1,N);
p=0;
while (1) {
cprintf(" %04d/%04d\b\b\b\b\b\b\b\b\b\b\b\b\b",p,N);
if (findpattern(y,p)) {
fitnum=0;
for (i=0;i<N;i++) {
for (j=first1(i);j<N-v[7];next1(i,j)) {
b=a[i][v[0]+j]+a[i][v[1]+j]+a[i][v[2]+j]+a[i][v[3]+j]+a[i][v[4]+j]+a[i][v[5]+j]+a[i][v[6]+j]+a[i][v[7]+j];
if (b==8) fitnum++;
}
}
if (fitnum>maxfitnum) {
bestv[0]=v[0];bestv[1]=v[1];bestv[2]=v[2];bestv[3]=v[3];bestv[4]=v[4];bestv[5]=v[5];bestv[6]=v[6];bestv[7]=v[7];
maxfitnum=fitnum;
}
p++;
if (p>N-8) break;
} else break;
}
cprintf("\r\n");
}
if (maxfitnum==0) {
printf("\nCan not find pattern!");
} else {
printf("\nbestv:");
for (i=0;i<8;i++) {
printf("%6d ",bestv[i]);
}
}
return 0;
}
uj5u.com熱心網友回復:
運算的次數好像也沒有理想中的減少很多!uj5u.com熱心網友回復:
我可以等后來者和我PK!
uj5u.com熱心網友回復:
回帖testuj5u.com熱心網友回復:
資料中只有0和1?uj5u.com熱心網友回復:
用動態規劃試試uj5u.com熱心網友回復:
資料中只有0和1uj5u.com熱心網友回復:
資料中只有0和1uj5u.com熱心網友回復:
有具體的例子或代碼嗎?uj5u.com熱心網友回復:
int array[8]{ 1,0,0,1,1,0,1,0 };
int tmp[8];
memset(tmp, 0, sizeof(int) * 8);
//找3個1
int sum = 0;
int i = 0;
for (; i < 8; i++)
{
sum += array[i];
if (sum == 3)
{
tmp[i] = 1;
break;
}
}
for (int j = i+1; j < 8; j++) {
if(array[j])
{
tmp[j] = tmp[j - 1] * 3+tmp[j-1];
}
else
{
tmp[j] = tmp[j - 1];
}
}
轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/280768.html
標籤:C語言
下一篇:關于宏定義與指標
