bzoj4514 [Sdoi2016]数字配对（网络流）

Description

有 n 种数字，第 i 种数字是 ai、有 bi 个，权值是 ci。

若两个数字 ai、aj 满足，ai 是 aj 的倍数，且 ai/aj 是一个质数，

那么这两个数字可以配对，并获得 ci×cj 的价值。

一个数字只能参与一次配对，可以不参与配对。

在获得的价值总和不小于 0 的前提下，求最多进行多少次配对。

 

Input

第一行一个整数 n。

第二行 n 个整数 a1、a2、……、an。

第三行 n 个整数 b1、b2、……、bn。

第四行 n 个整数 c1、c2、……、cn。

 

 

Output

 一行一个数，最多进行多少次配对

 

Sample Input

3
2 4 8
2 200 7
-1 -2 1

Sample Output

4

 

HINT

 n≤200，ai≤10^9，bi≤10^5，∣ci∣≤10^5

 

Source

鸣谢Menci上传

 

 

主要是要注意到质因数有几个（2^2算两个），这样个数为奇数只能与个数为偶数配对，个数为偶数只能和个数为奇数配对，然后就是一个二分图，随便建一建图跑最大费用流就好。

program rrr(input,output);
const
  inf=123456789012345;
type
  etype=record
     t,c,next,rev:longint;
     w:int64;
  end;
var
  e:array[0..100000]of etype;
  num,a,d,q,fre,frv:array[0..220]of longint;
  s:array[0..100010]of boolean;
  p:array[0..100010]of longint;
  c,dis:array[0..220]of int64;
  inq:array[0..220]of boolean;
  n,m,i,j,x,b,cnt,h,t,ans:longint;
  w,f:int64;
function min(a,b:int64):int64;
begin
   if a<b then exit(a) else exit(b);
end;
procedure ins(x,y,c:longint;w:int64);
begin
   inc(cnt);e[cnt].t:=y;e[cnt].c:=c;e[cnt].w:=w;e[cnt].next:=a[x];a[x]:=cnt;
end;
procedure add(x,y,c:longint;w:int64);
begin
   ins(x,y,c,w);e[cnt].rev:=cnt+1;ins(y,x,0,-w);e[cnt].rev:=cnt-1;
end;
procedure spfa;
begin
   for i:=1 to n do dis[i]:=-inf;dis[0]:=0;
   h:=0;t:=1;q[1]:=0;inq[0]:=true;
   while h<>t do
      begin
         inc(h);if h>210 then h:=1;
         i:=a[q[h]];
         while i<>0 do
            begin
               if (e[i].c>0) and (dis[q[h]]+e[i].w>dis[e[i].t]) then
                  begin
                     dis[e[i].t]:=dis[q[h]]+e[i].w;
                     fre[e[i].t]:=i;frv[e[i].t]:=q[h];
                     if not inq[e[i].t] then
                        begin
                           inc(t);if t>210 then t:=1;
                           q[t]:=e[i].t;inq[e[i].t]:=true;
                        end;
                  end;
               i:=e[i].next;
            end;
         inq[q[h]]:=false;
      end;
end;
begin
   assign(input,'r.in');assign(output,'r.out');reset(input);rewrite(output);
   fillchar(s,sizeof(s),true);s[1]:=false;
   for i:=2 to 320 do if s[i] then
      begin
         j:=i+i;while j<=100000 do begin s[j]:=false;j:=j+i; end;
      end;
   m:=0;for i:=2 to 100000 do if s[i] then begin inc(m);p[m]:=i; end;
   readln(n);
   for i:=1 to n do read(num[i]);
   for i:=1 to n do
      begin
         x:=num[i];j:=1;d[i]:=0;
         while x>1 do begin while x mod p[j]=0 do begin inc(d[i]);x:=x div p[j]; end;inc(j);if j>m then break; end;
         if x>1 then inc(d[i]);
      end;
   fillchar(a,sizeof(a),0);cnt:=0;
   for i:=1 to n do begin read(b);if d[i] mod 2=0 then add(i,n+1,b,0) else add(0,i,b,0); end;
   for i:=1 to n do read(c[i]);
   for i:=1 to n do for j:=i+1 to n do
      if (abs(d[i]-d[j])=1) and ((num[i] mod num[j]=0) or (num[j] mod num[i]=0)) then
         begin
            if d[i] mod 2=0 then add(j,i,100000000,c[i]*c[j]) else add(i,j,100000000,c[i]*c[j]);
         end;
   ans:=0;inc(n);w:=0;
   while true do
      begin
         spfa;
         if dis[n]=-inf then break;
         i:=n;f:=100000000;
         while i<>0 do begin f:=min(f,e[fre[i]].c);i:=frv[i]; end;
         if w+f*dis[n]<0 then begin ans:=ans+w div (-dis[n]);break; end
         else begin ans:=ans+f;w:=w+f*dis[n]; end;
         i:=n;while i<>0 do begin dec(e[fre[i]].c,f);inc(e[e[fre[i]].rev].c,f);i:=frv[i]; end;
      end;
   write(ans);
   close(input);close(output);
end.