hdu 3335 Divisibility

最小路径覆盖

先把相同的数字去掉即去重，按整除或者被整除关系建立有向图，可知这个有向图一定是无环的（DAG），转化为最小路径覆盖模型

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define MAX 1010

bool vis[MAX];
int mat[MAX];
long long a[MAX];
int first[MAX];
struct edge
{
   int u,v,next;
}e[MAX*MAX];
int nume;

void add(int u , int v)
{
   e[nume].u=u; e[nume].v=v;
   e[nume].next=first[u]; first[u]=nume++;
}

bool dfs(int u)
{
   for(int k=first[u]; k!=-1; k=e[k].next)
      if(!vis[e[k].v])
      {
         int v=e[k].v;
         vis[v]=true;
         if(mat[v]==-1 || dfs(mat[v]))
         {
            mat[v]=u;
            return true;
         }
      }
   return false;
}

int main()
{
   int T,n;
   scanf("%d",&T);
   while(T--)
   {
      scanf("%d",&n);
      for(int i=1; i<=n; i++)
         scanf("%I64d",&a[i]);
      sort(a+1,a+n+1);
      int m=1;
      for(int i=2; i<=n; i++)
         if(a[i]!=a[i-1])
            a[++m]=a[i];
      n=m;

      nume=0;
      memset(first,-1,sizeof(first));
      for(int i=1; i<=n; i++)
         for(int j=1; j<=n; j++)
            if(i!=j && !(a[i]%a[j]))
             add(i,j);

      memset(mat,-1,sizeof(mat));
      int ccount=0;
      for(int i=1; i<=n; i++)
      {
         memset(vis,false,sizeof(vis));
         if(dfs(i))
            ccount++;
      }
      printf("%d\n",n-ccount);
   }
   return 0;
}