数据挖掘-Apriori算法

算法  Apriori

Input:          Transaction DataBase DMinimum support threshold minsup

Output      Frequent pattern L

(1) L1=search_frequent_1-itemsets( D );

(2) for(k=2;Lk-1φ;k++do

(3) begin

(4)    Ck=apriori-gen(Lk-1);

(5)    for all transactions t D do

(6)    begin

(7)      Ct=subset(Ckt);

(8)      for all candidates c Ct do

(9)        c.count++;

(10)    end

(11)    Lk ={c Ck|c.countminsup}

(12) end

(13) Answer L=kLk;

Procedure Search_frequent_1-itemsets( D )

(1) begin

(2)  for all transactions t D do

(3)  begin

(4)    for each item ik t do

(5)      ik.count++;

(6)  end

(7)  L1 ={ i | i.countminsup}

(8)  return L1;

(9) end

Procedure apriori_gen(Lk)

(1) begin

(2)   for each itemset l1 Ldo

(3)     for each itemset l2 Ldo

(4)     begin

(5)       if l1[1]=l2[1]) l1[2]=l2[2])  l1[k-1]=l2[k-1]) l1[k]<l2[k]) then

(6)       begin

(7)          c= l1 l2;

(8)          if Is_include_infrenquent_subset(c,Lkthen

(9)             delete c;

(10)         else add c to Ck+1 ;

(11)       end

(12)      end

(13)    return Ck+1 ;

(14) end

Procedure Is_include_infrenquent_subset(c,Lk)

(1)begin

(2)  for each k-subset s of c

(3)     if s Lk then

(4)       reture TURE;

(5)  return FALSE ;

(6)end

在主程序中,第一步首先扫描整个交易数据库D,统计每个项目(item)的支持数,计算其支持度,将支持度大于等于最小支持度minsup的项目构成的集合放入到L中;从第2步到第11步,用k-1频繁项目集构成的Lk-1生成候选集的集合Ck,以便从中生成Lk,其中apriori_gen函数(4)用来从Lk-1中生成Ck,然后对数据库进行扫描(5),对于数据库中的每一个交易,subset函数用来发现此交易包含的所有候选集(7),并为这些候选集的计数器加1(8-9)。最后满足minsup的候选集被放入到Lk中。

apriori_gen过程完成两种操作:并(join)和剪枝(prune)。在并运算步骤中,Lk-1 Lk-1 进行并运算生成潜在的候选集(2-7),条件l1[k-1]<l2[k-1]保证不会有重复的候选集生成(5)。在剪枝步骤中(8-10),利用性质2.1,删除那些存在子集不是频繁项目集的候选集,测试子集是否为频繁项目集由过程Is_include_infrenquent_subset完成。

为了清楚的阐述Apriori算法的挖掘过程,现举例如下:

【例1设事务数据库D如表2.1所示,D中包含4个事务,即|D|=4,最小支持数mincount=2,即最小支持度minsup=2/4=50%。挖掘频繁项目集的具体过程如下所述:C1={{A},{B},{C},{D},{F}},第一次循环产生L1={{A},{B},{C},{F}},由Apriori_gen(L1)生成C2,扫描数据库,计算C2中每个候选集得到L2。依此循环,得到L3。整个挖掘过程如图2.1所示。

事务数据库D

Tid

事务

100

200

300

400

B,C,F

A,C,D

B,F

A,B,C,F

   

 

1 Apriori算法的执行过程

在找到了事务数据库中的所有频繁项集后,利用这些频繁项集可以产生关联规则,产生关联规则的步骤如下:

(1) 对于每个频繁项目集l,产生l的所有非空子集。

(2) 对于l的每个非空子集m,如果support(l)/support(m)≥minconf,则输出规则m (l-m)”

例如,在上例中产生的频繁项目集l{B,C,F}l的非空子集有{B,C}{B,F}{C,F}{B}{C}{F},则运用上述产生关联规则的方法可以得到以下关联规则:

          B C    confidence=(2/4)/(4/4)=1

          B F    confidence=(2/4)/(3/4)=0.667

          C F    confidence=(2/4)/(2/4)=1

          F B    confidence=(2/4)/(3/4)= 0.667

          C B    confidence=(2/4)/(3/4)= 0.667

          B C    confidence=(2/4)/(3/4)= 0.667

 

 

 

源代码  apriori.c

//

#include<stdio.h>
typedef  struct 
{
int item[100];  //数据项
} D_Node;      //数据库D


typedef  struct 
{
int item[100];  //数据项,用item[0]保存支持度

} C_Node; //候选集

typedef  struct 
{
int item[100];  //数据项,用item[0]保存支持度
} L_Node;//频繁集

C_Node C[100][100];
L_Node L[100][100];
D_Node D[100];

int min_supp;  //最小支持度

 

void  InPut()
{
 int i,j,n,n1;
 printf("请输入最小支持度:");
 scanf("%d",&min_supp);
  printf("请输入交易集的大小");
  scanf("%d",&D[0].item[0]);
  n=D[0].item[0];

     for(i=1;i<=n;i++)  //for1
     {   
        printf("请输入交易[%d]中记录的个数(n)",i);
         scanf("%d",&n1);
         D[i].item[0]=n1;

           for(j=1;j<=n1;j++)  //for2
           {
              printf("请输入交易[%d]中记录项,直接输入数字:",i);
              scanf("%d",&D[i].item[j]);                  
           }//for2
                 
           }  //for1
     
     }//end of InPut


void  C1()
{
 //功能:扫描数据集D生成1项候选集C1
 //输入:数据集D
 //输出1项候选集C1
 //初始条件 数据集D 非空
   int i,j,k;
   int no=1,temp=0;
   C[1][0].item[0]=0;  //1 项集的个数,在本算法中,用C[n][k].item[0]来保存候选集Cn的第k项的支持度
   if(D[0].item[0]!=0)
   {  
      C[1][1].item[1]=D[1].item[1];
          
   }
      
    for(i=1;i<=D[0].item[0];i++)  //for1  
    { 
  
         for(j=1;j<=D[i].item[0];j++)  //for2   
      {
               temp=1;
               for(k=1;k<=no;k++)  //for3
      { 
                  if(C[1][k].item[1]==D[i].item[j])
            {
                     C[1][k].item[0]++;  //支持度加1
                      temp=0;  //
                           
                  }   //if

               }//end for3
            
          
              if(temp)//生成新的项集
           { 
                    C[1][++no].item[1]=D[i].item[j];
                    C[1][no].item[0]=1;     
           } 
       
         }//end for2
   
    } // end  for1

    C[1][0].item[0]=no;//数据项的个数 

 //end of  C1()


void Cn( int n)
 
 //用频繁集Ln-1为基础,通过连接得到n项候选集Cn
 
 int i,j,k,p,q,s,t,num;
 int no=0,temp=0,count;

 C[n][0].item[0]=0;  //初始化

//printf("in Cn(%d) n=%d/n",n,n);
//printf("in Cn(%d) C[%d][0].item[0]=%d/n",n,n,C[n][0].item[0]);

num=L[n-1][0].item[0];  //num是Ln-1项集的数据个数

   for(i=1;i<=num;i++)  
  
      for(j=i+1;j<=num;j++)   //for2
         {
   
      temp=1;  //测试是否满足联结条件
              if(n>2)//if 1
              {
                  for(k=1;k<n-1;k++)     //for3   
                 {
      if(L[n-1][i].item[k]!=L[n-1][j].item[k])
                   {  temp=0;
                     break; }//if 1

                    }//end for3

                 }//end if1

                      if(temp==1)//满足联结条件
                     {
                          // printf("in if 2  no=%d/n",no);
                     no++;

                    for(p=1;p<=n-1;p++)  
                    C[n][no].item[p]=L[n-1][i].item[p];
                    C[n][no].item[p]=L[n-1][j].item[p-1];
                    C[n][no].item[0]=0;
                         for(q=1;q<=D[0].item[0];q++)  //for5  测试其支持度

       {    
        
        count=0; //count用来记数,当所测试的项存在时,count加1,当count=n时,则子集存在

                               for(s=1;C[n][no].item[s]!=0;s++)  //for6
                          {
                                      for(t=1;t<=D[q].item[0];t++)  //for7
           {
                                             if(C[n][no].item[s]==D[q].item[t])

            {   count+=1;
                      break;
            } 
           }//end for7
                                     
          }//end for 6
                           if(count==n) C[n][no].item[0]+=1;//子集存在,第no项的支持度加1
                   
       }//end for5
                           
                           
                      C[n][0].item[0]+=1;
       }//end if2
           }//end for2

    
   
}//end of Cn()

void L1()
  
    int i,j,k;
 j=0;
   L[1][0].item[0]=0;
   
   //printf("C[1][0].item[0]=%d/n",C[1][0].item[0]);

   for(i=1;i<=C[1][0].item[0];i++)
   {    
       if(C[1][i].item[0]>=min_supp)
       {
       j+=1;
        for(k=1;k<=1;k++)
         L[1][j].item[k]=C[1][i].item[k];
        L[1][j].item[0]=C[1][i].item[0];
 // printf("L[1][%d].item[1]=%d   ",j,L[1][j].item[1]);  测试功能时加的
 // printf("  -------------%d/n",L[1][j].item[0]);
  
       }
   }//end for1
   L[1][0].item[0]=j;
}//end of L1()

void Ln(int n)
  
    int i,j,k;
 Cn(n);
    j=0;
    L[n][0].item[0]=0;

   // printf("in Ln(%d) C[%d][0].item[0]=%d/n",n,n,C[n][0].item[0]);

   for(i=1;i<=C[n][0].item[0];i++)  //for 1
   {
       if(C[n][i].item[0]>=min_supp)
       {
         j+=1;
        for(k=1;k<=n;k++)
         L[n][j].item[k]=C[n][i].item[k];
        L[n][j].item[0]=C[n][i].item[0]; 
       }  //end if


   }//end for1

  

L[n][0].item[0]=j; //保存数据的个数

}//end of Ln(int n)


  void  OutPut(int n)
  {
     int i,j,k;
     printf("频繁项目集L%d如下:/n",n);
     k=L[n][0].item[0];
         if(k!=0)
         { 
             for(i=1;i<=k;i++)
          {
                printf("{");
               for(j=1;j<=n;j++)      
        printf("  I%d ",L[n][i].item[j]);
        printf("}        支持度:%d/n",L[n][i].item[0]);
           
      
             }//for
            
          }
         else                printf("项目集为空/n");
           
      
  }


  void main()
  {
 int i;
 int n=1;
    InPut();
 C1();//初始化,生成1项候选集C1
 L1();//得到1项频繁集L1
    while(L[n][0].item[0]!=0)
 {
  n+=1;
        Ln(n);
 }
    for(i=1;i<=n;i++)
        OutPut(i);

 char ch;
   scanf("%d",&i);
  
  }

 

 

--------------------------------------------------------------------------------------

FAST apriori.cpp

//

#include <stdio.h>
#include <string.h>

typedef  struct 
{
char item[10];  //数据项
int min_supp_count;//最小支持度数
} C_Node;      //候选集


typedef  struct 
{
char item[10];  //数据项
int min_supp_count;//最小支持度数
} L_Node;     //频繁集

char D[10][10];
L_Node L[100];
C_Node C[100];
int min_supp_count=2;
int num=100;
void  InPut()
{

 strcpy(D[1],"abe");
 strcpy(D[2],"bd");
 strcpy(D[3],"bc");
 strcpy(D[4],"abd");
 strcpy(D[5],"ac");
 strcpy(D[6],"bc");
 strcpy(D[7],"ac");
 strcpy(D[8],"abce");
 strcpy(D[9],"abc");

 
}//end of InPut
int  * DB=new int[num];

 

void suppDB()
{
  int m='e';
  int n;
  int k;
  for (int i=1;i<=9;i++)
  {
   n=strlen(D[i]);
   DB [i]=0;
   for (int j=0;j<n;j++)
   {
    k=1; 
    DB [i]+=k<<(int)(m-D[i][j]);  
   }  
  }

 

}

void check_supp(int num,int no)
{
int i,j,k,m;
    int check;
   m='e';
   
     for(i=1;i<=num;i++)
  {    check=0;
       C[i].min_supp_count=0;
       for (j=0;j<no;j++)
    {
   k=1;
   check+=(int)(k<<(m-C[i].item[j]));
      }
   for (j=1;j<=9;j++)
   {
    if (check==(check&DB[j]))
    {
      C[i].min_supp_count+=1;//子集存在,支持度数加1
    }
   }
  }

}

void  C1()
{
 //功能:扫描数据集D生成1项候选集C1
 //输入:数据集D
 //输出1项候选集C1
 //初始条件 数据集D 非空
 strcpy(C[1].item,"a");
 strcpy(C[2].item,"b");
 strcpy(C[3].item,"c");
 strcpy(C[4].item,"d");
 strcpy(C[5].item,"e");
   
   C[0].min_supp_count=5;  //1 项候选集的个数,在本算法中,用C[0].min_supp_count来保存候选集Cn的个数  

  check_supp(5,1);
  
   
  //end of  C1()


void Cn( int n)
 
 //用频繁集Ln-1为基础,通过连接得到n项候选集Cn
 
 int i,j,k,p,num;
 int no=0,temp=0;

 C[0].min_supp_count=0;  //初始化


num=L[0].min_supp_count;  //num是Ln-1项集的数据个数

   for(i=1;i<=num;i++)  
  
      for(j=i+1;j<=num;j++)   //for2
         {
   
      temp=1;  //测试是否满足联结条件
              if(n>2)//if 1
              {
                  for(k=0;k<n-2;k++)     //for3   
                 {
      if(L[i].item[k]!=L[j].item[k])
                   {  temp=0;
                     break; }//if 1

                    }//end for3

                 }//end if1

                      if(temp==1)//满足联结条件
                     {
                          // printf("in if 2  no=%d/n",no);
                     no++;

                    for(p=0;p<=n-2;p++)  
                    C[no].item[p]=L[i].item[p];
                    C[no].item[p]=L[j].item[p-1];
                    C[no].min_supp_count=0;
                    C[0].min_supp_count+=1;
       }//end if2
           }//end for2
num=C[0].min_supp_count;
   check_supp(num,n);//测试支持度 
}//end of Cn() 


void L1()
  int n=1;
    int i,j,k;
 j=0;
   L[0].min_supp_count=0;//频繁集的个数,初始为0
    
   for(i=1;i<=C[0].min_supp_count;i++)
   {    
       if(C[i].min_supp_count>=min_supp_count)
       {
       j+=1;
        strcpy(L[j].item,C[i].item);  
  L[j].min_supp_count=C[i].min_supp_count;
       }
   }//end for1
   L[0].min_supp_count=j;///频繁集的个数,最后为j个

    printf("频繁项目集L%d如下:/n",n);
     k=L[0].min_supp_count;
         if(k!=0)
         { 
             for(i=1;i<=k;i++)
          {
                printf("{");
               for(j=0;j<n;j++)      
        printf("  %c ",L[i].item[j]);
        printf("}        支持度:%d/n",L[i].min_supp_count);
           
      
             }//for
            
          }
         else                printf("项目集为空/n");

}//end of L1()

void Ln(int n)
  
    int i,j,k;
 Cn(n);
    j=0;
    L[0].min_supp_count=0;

   
   for(i=1;i<=C[0].min_supp_count;i++)  //for 1
   {
       if(C[i].min_supp_count >=min_supp_count)
       {
         j+=1;
        strcpy(L[j].item,C[i].item);  
  L[j].min_supp_count=C[i].min_supp_count;
         
       }  //end if


   }//end for1

  

L[0].min_supp_count=j; //保存数据的个数
 printf("频繁项目集L%d如下:/n",n);
     k=L[0].min_supp_count;
         if(k!=0)
         { 
             for(i=1;i<=k;i++)
          {
                printf("{");
               for(j=0;j<n;j++)      
        printf("  %c ",L[i].item[j]);
        printf("}        支持度:%d/n",L[i].min_supp_count);
           
      
             }//for
            
          }
         else                printf("项目集为空/n");

}//end of Ln(int n)


 void main()
{

   
      int n=1;
    InPut();
 suppDB();
 C1();//初始化,生成1项候选集C1
 L1();//得到1项频繁集L1
  while(L[0].min_supp_count!=0)
 {
  n+=1;
        Ln(n);
  }
  char ch;
  printf("press any key to eixe/n");
  scanf("%c",&ch);
   
   }

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值