原始状态
#include <stdio.h>
#include <vector>
using namespace std;
const int INF=1000;
const float lrate=1;
int W[INF],X[INF][INF];
int XS[INF];
int XN,WN,B=0;
int err(int p){
int sum=0;
for(int i=0;i<WN;i++){
sum+=W[i]*X[p][i];
}
sum=(sum+B)*XS[p];
if(sum<=0)
return 1;
return 0;
}
int main(){
//freopen("in.txt","r",stdin);
scanf("%d",&WN);
scanf("%d",&XN);
for(int i=0;i<XN;i++){
for(int j=0;j<WN;j++){
scanf("%d",&X[i][j]);
}
scanf("%d",&XS[i]);
}
for(int i=0;i<WN;i++){
W[i]=0;
}
for(int i=0;i<XN;i++){
if(err(i)){
for(int j=0;j<WN;j++){
W[j]+=(lrate*X[i][j]*XS[i]);
}
B+=(lrate*XS[i]);
printf("%d\t%d\t%d\t%d\t\n",i,W[0],W[0],B);
i=-1;
}
}
return 0;
}
对偶形式
#include <stdio.h>
#include <string.h>
const int INF=1000;
int X[INF][INF],Y[INF],gram[INF][INF];
int A[INF],B=0,XN,WN,lrate=1;
int gramini(){
for(int i=0;i<XN;i++){
for(int j=0;j<XN;j++){
for(int k=0;k<WN;k++){
gram[i][j]+=X[i][k]*X[j][k];
}
}
}
return 0;
}
int err(int p){
int sum=0;
for(int i=0;i<XN;i++){
sum+=(A[i]*Y[i]*gram[i][p]);
}
sum=(sum+B)*Y[p];
if(sum<=0)
return 1;
return 0;
}
int main(){
memset(gram,0,sizeof(gram));
memset(A,0,sizeof(A));
memset(X,0,sizeof(X));
freopen("in.txt","r",stdin);
scanf("%d",&WN);
scanf("%d",&XN);
for(int i=0;i<XN;i++){
for(int j=0;j<WN;j++){
scanf("%d",&X[i][j]);
}
scanf("%d",&Y[i]);
}
gramini();
for(int i=0;i<XN;i++){
if(err(i)){
A[i]++;
B+=Y[i];
printf("%d\t%d\t%d\t%d\t%d\n",i,A[0],A[1],A[2],B);
i=-1;
}
}
return 0;
}