BZOJ

洛谷

\(Solution\)

很显然的建二分图后跑最大费用流，但有个问题是一个数是只能用一次的，这样二分图两部分都有这个数。

那么就用两倍的。如果\(i\)可以向\(j'\)连边，\(j\)也向\(i'\)连边，如果上一次走了\(i->j'\)，那么这一次一定走\(j->i'\)。

每次跑最大费用流，直至有一次费用变成负，然后加上当前正权值能抵消它的流量，最后总流量除以2就可以了。

\(Another Solution\)

两个数能匹配首先要能整除，其次它们所有质因子的次数和一定只相差1.

于是可以按这个次数的奇偶性建二分图，分别连源点汇点，一遍最大流就可以了。

ps：sb的我硬是把\(i\to j\)的费用设成\(0\)，然后源点汇点连\(\log b[i]\)的边，然后再取幂。。最后写完发现不对，不会处理上面的问题。。贪心也调不对。

数字匹配\(O(sqrt(n))\)直接枚举就可以，闲的写什么MillerRabin+set。。

考试的时候一个题写了三百多行 还不算颓

时隔一年后的考试写出来了2333（虽然写过还是道水题）

第一种方法：

  1. //3212kb	112ms
  2. #include <set>
  3. #include <queue>
  4. #include <cmath>
  5. #include <cstdio>
  6. #include <cctype>
  7. #include <algorithm>
  8. #define gc() getchar()
  9. typedef long long LL;
  10. const int N=410,M=100005,INF=1e9;
  11. const int P[9]={2,3,5,7,13,17,31};
  13. int n,src,des,A[N],B[N],C[N],Enum,H[N],nxt[M],fr[M],to[M],cap[M],pre[N];
  14. bool inq[N];
  15. LL dis[N],cost[M];
  16. std::queue<int> q;
  17. std::set<int> is_p[2];
  19. inline int read()
  20. {
  21. 	int now=0,f=1;register char c=gc();
  22. 	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
  23. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  24. 	return now*f;
  25. }
  26. inline void AddEdge(int u,int v,int w,LL c)
  27. {
  28. 	to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w, cost[Enum]=c;
  29. 	to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0, cost[Enum]=-c;
  30. }
  31. int FP(LL x,int k,int p)
  32. {
  33. 	LL t=1;
  34. 	for(; k; k>>=1,x=x*x%p)
  35. 		if(k&1) t=t*x%p;
  36. 	return (int)t;
  37. }
  38. bool Miller_Rabin(int x)
  39. {
  40. 	if(x==2) return 1;
  41. 	if(!(x&1)||x==1) return 0;
  42. 	for(int i=0; i<7; ++i)
  43. 		if(x==P[i]) return 1;
  44. 		else if(!(x%P[i])) return 0;
  45. 	int u=x-1,t=0;
  46. 	while(!(u&1)) u>>=1,++t;
  47. 	for(int i=0; i<7; ++i)
  48. 	{
  49. 		LL now=FP(P[i],u,x),las;
  50. 		for(int j=1; j<=t; ++j)
  51. 		{
  52. 			las=now, (now*=now)%=x;
  53. 			if(now==1&&las!=1&&las!=x-1) return 0;
  54. 		}
  55. 		if(now!=1) return 0;
  56. 	}
  57. 	return 1;
  58. }
  59. bool SPFA()
  60. {
  61. 	for(int i=1; i<=des; ++i) dis[i]=-1e12;
  62. 	dis[src]=0, q.push(src);
  63. 	while(!q.empty())
  64. 	{
  65. 		int x=q.front(); q.pop();
  66. 		inq[x]=0;
  67. 		for(int v,i=H[x]; i; i=nxt[i])
  68. 			if(cap[i]&&dis[to[i]]<dis[x]+cost[i])
  69. 			{
  70. 				dis[v=to[i]]=dis[x]+cost[i], pre[v]=i;
  71. 				if(!inq[v]) q.push(v),inq[v]=1;
  72. 			}
  73. 	}
  74. 	return dis[des]>-1e12;
  75. }
  76. LL Solve()
  77. {
  78. 	int flow=0,mn;
  79. 	LL sum_cost=0;
  80. 	while(SPFA())
  81. 	{//注意这只是一个二分图。。
  82. 		mn=1e9;
  83. 		for(int i=des; i!=src; i=fr[pre[i]]) mn=std::min(cap[pre[i]],mn);
  84. 		if(sum_cost+1ll*dis[des]*mn>=0){//dis[des]就是费用。
  85. 			for(int i=des; i!=src; i=fr[pre[i]]) cap[pre[i]]-=mn,cap[pre[i]^1]+=mn;
  86. 			sum_cost+=1ll*dis[des]*mn, flow+=mn;
  87. 		}
  88. 		else{
  89. 			flow+=(int)(sum_cost/std::abs(dis[des]));//加上剩余正价值能补多少负值。
  90. 			break;
  91. 		}
  92. 	}
  93. 	return flow>>1;
  94. }
  96. int main()
  97. {
  98. 	n=read(),Enum=1,src=0,des=n<<1|1;
  99. 	for(int i=1; i<=n; ++i) A[i]=read();
  100. 	for(int i=1; i<=n; ++i) B[i]=read();
  101. 	for(int i=1; i<=n; ++i) C[i]=read();
  102. 	for(int i=1; i<=n; ++i) AddEdge(src,i,B[i],0),AddEdge(i+n,des,B[i],0);
  103. 	for(int i=1; i<=n; ++i)
  104. 		for(int j=1; j<=n; ++j)
  105. 			if(A[i]>A[j]&&!(A[i]%A[j]))
  106. 			{
  107. 				int t=A[i]/A[j];
  108. 				if(is_p[0].count(t)) continue;
  109. 				if(is_p[1].count(t)) AddEdge(i,j+n,INF,1ll*C[i]*C[j]),AddEdge(j,i+n,INF,1ll*C[i]*C[j]);
  110. 				else{
  111. 					bool f=Miller_Rabin(t);
  112. 					if(f) is_p[1].insert(t),AddEdge(i,j+n,INF,1ll*C[i]*C[j]),AddEdge(j,i+n,INF,1ll*C[i]*C[j]);
  113. 					else is_p[0].insert(t);
  114. 				}
  115. 			}
  116. 	printf("%lld",Solve());
  118. 	return 0;
  119. }

第二种方法（19.2.17）：

  1. /*
  2. 2840kb	28ms
  3. 常规做建二分图跑费用流很容易，但因为拆了点每对数之间可能选两次。但是既然选了一次那就一定会选第二次。最后费用即将变成负数时判一下还能选多少个，然后流量除以二即可。
  4. 因为匹配条件是$i$可整除$j$且结果是质数，那么$i,j$之间质因子次数一定相差一，可以直接判，同时根据质因子次数奇偶性可以分别连源点汇点，就没有之前的问题了。
  5. 另外这张二分图里费用就是$dis[T]$。
  6. */
  7. #include <queue>
  8. #include <cstdio>
  9. #include <cctype>
  10. #include <algorithm>
  11. #define gc() getchar()
  12. typedef long long LL;
  13. const int N=205,M=(N*N+N)*2,INF=1<<30;
  14. const LL INFll=1ll<<60;
  16. int T,Enum,H[N],nxt[M],fr[M],to[M],cap[M],pre[N];
  17. LL dis[N],cost[M];
  19. inline int read()
  20. {
  21. 	int now=0,f=1;register char c=gc();
  22. 	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
  23. 	for(;isdigit(c);now=now*10+c-48,c=gc());
  24. 	return now*f;
  25. }
  26. inline void AE(int u,int v,int w,LL c)
  27. {
  28. 	to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w, cost[Enum]=c;
  29. 	to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0, cost[Enum]=-c;
  30. }
  31. int Calc(int x)
  32. {
  33. 	int t=0;
  34. 	for(int i=2; 1ll*i*i<=x; ++i)
  35. 		if(!(x%i))
  36. 		{
  37. 			x/=i, ++t;
  38. 			while(!(x%i)) x/=i, ++t;
  39. 		}
  40. 	if(x!=1) ++t;
  41. 	return t;
  42. }
  43. bool SPFA()
  44. {
  45. 	static bool inq[N];
  46. 	static std::queue<int> q;
  47. 	for(int i=1; i<=T; ++i) dis[i]=-INFll;
  48. 	q.push(0), dis[0]=0;
  49. 	while(!q.empty())
  50. 	{
  51. 		int x=q.front(); q.pop();
  52. 		inq[x]=0;
  53. 		for(int i=H[x],v; i; i=nxt[i])
  54. 			if(dis[v=to[i]]<dis[x]+cost[i] && cap[i])
  55. 				pre[v]=i, dis[v]=dis[x]+cost[i], !inq[v]&&(q.push(v),inq[v]=1);
  56. 	}
  57. 	return dis[T]>-INFll;
  58. }
  59. int Solve()
  60. {
  61. 	int res=0; LL now=0;
  62. 	while(SPFA())
  63. 	{
  64. 		int mn=INF;
  65. 		for(int i=T; i; i=fr[pre[i]]) mn=std::min(mn,cap[pre[i]]);
  66. 		if(now+dis[T]*mn>=0)
  67. 		{
  68. 			now+=dis[T]*mn, res+=mn;
  69. 			for(int i=T; i; i=fr[pre[i]]) cap[pre[i]]-=mn, cap[pre[i]^1]+=mn;
  70. 		}
  71. 		else
  72. 		{
  73. 			res+=now/std::abs(dis[T]);//abs!!
  74. 			break;
  75. 		}
  76. 	}
  77. 	return res;
  78. }
  80. int main()
  81. {
  82. 	static int A[N],B[N],C[N],tm[N];
  84. 	freopen("pair.in","r",stdin);
  85. 	freopen("pair.out","w",stdout);
  87. 	int n=read(); Enum=1, T=n+1;
  88. 	for(int i=1; i<=n; ++i) tm[i]=Calc(A[i]=read());
  89. 	for(int i=1; i<=n; ++i) B[i]=read();
  90. 	for(int i=1; i<=n; ++i) C[i]=read();
  91. 	for(int i=1; i<=n; ++i)
  92. 	{
  93. 		if(tm[i]&1)
  94. 		{
  95. 			AE(0,i,B[i],0);
  96. 			for(int j=1; j<=n; ++j)
  97. 				if(!(tm[j]&1) && tm[i]-tm[j]==1 && !(A[i]%A[j])) AE(i,j,INF,1ll*C[i]*C[j]);
  98. 		}
  99. 		else
  100. 		{
  101. 			AE(i,T,B[i],0);
  102. 			for(int j=1; j<=n; ++j)
  103. 				if(tm[j]&1 && tm[i]-tm[j]==1 && !(A[i]%A[j])) AE(j,i,INF,1ll*C[i]*C[j]);
  104. 		}
  105. 	}
  106. 	printf("%d\n",Solve());
  108. 	return 0;
  109. }

弃疗的暴力+最大流（18.4.1）：

  1. #include <set>
  2. #include <queue>
  3. #include <cstdio>
  4. #include <cctype>
  5. #include <algorithm>
  6. #define gc() getchar()
  7. typedef long long LL;
  8. const int N=410,M=100005,INF=1e9;
  9. const int P[9]={2,3,5,7,13,17,31};
  11. int n,range,src,des,A[N],B[N],C[N],Enum,cur[N],H[N],nxt[M],fr[M],to[M],cap[M],pre[N],lev[N],num[N],que[N];
  12. std::queue<int> q;
  13. std::set<int> is_p[2];
  15. inline int read()
  16. {
  17. 	int now=0,f=1;register char c=gc();
  18. 	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
  19. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  20. 	return now*f;
  21. }
  22. inline void AddEdge(int u,int v,int w)
  23. {
  24. 	printf("%d->%d %d\n",u,v,w);
  25. 	to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w;
  26. 	to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0;
  27. }
  28. int FP(LL x,int k,int p)
  29. {
  30. 	LL t=1;
  31. 	for(; k; k>>=1,x=x*x%p)
  32. 		if(k&1) t=t*x%p;
  33. 	return (int)t;
  34. }
  35. bool Miller_Rabin(int x)
  36. {
  37. 	if(x==2) return 1;
  38. 	if(!(x&1)||x==1) return 0;
  39. 	for(int i=0; i<7; ++i)
  40. 		if(x==P[i]) return 1;
  41. 		else if(!(x%P[i])) return 0;
  42. 	int u=x-1,t=0;
  43. 	while(!(u&1)) u>>=1,++t;
  44. 	for(int i=0; i<7; ++i)
  45. 	{
  46. 		LL now=FP(P[i],u,x),las;
  47. 		for(int j=1; j<=t; ++j)
  48. 		{
  49. 			las=now, (now*=now)%=x;
  50. 			if(now==1&&las!=1&&las!=x-1) return 0;
  51. 		}
  52. 		if(now!=1) return 0;
  53. 	}
  54. 	return 1;
  55. }
  56. namespace Subtask1
  57. {
  58. 	int Ans;
  59. 	bool can[13][13],vis[13];
  61. 	void DFS(int x,int tm,LL sum)
  62. 	{//O(n!)
  63. 		if(x>n)
  64. 			if(sum>=0 && tm>Ans) Ans=tm;
  65. 			else ;
  66. 		else{
  67. 			if(vis[x]) DFS(x+1,tm,sum);
  68. 			else{
  69. 				vis[x]=1;
  70. 				for(int i=1; i<=n; ++i)
  71. 					if(can[x][i]&&!vis[i])
  72. 						vis[i]=1, DFS(x+1,tm+1,sum+1ll*C[x]*C[i]), vis[i]=0;
  73. 				vis[x]=0;
  74. 				DFS(x+1,tm,sum);
  75. 			}
  76. 		}
  77. 	}
  78. 	void Solve()
  79. 	{
  80. 		for(int i=1; i<=n; ++i) A[i]=read();
  81. 		for(int i=1; i<=n; ++i) B[i]=read();
  82. 		for(int i=1; i<=n; ++i) C[i]=read();
  83. 		for(int i=1; i<=n; ++i)
  84. 			for(int j=1; j<=n; ++j)
  85. 				if(!(A[i]%A[j])&&i!=j)
  86. 				{
  87. 					int t=A[i]/A[j];
  88. 					if(is_p[0].count(t)) continue;
  89. 					if(is_p[1].count(t)) can[i][j]=1;
  90. 					else{
  91. 						bool f=Miller_Rabin(t);
  92. 						if(f) is_p[1].insert(t),can[i][j]=1;
  93. 						else is_p[0].insert(t);
  94. 					}
  95. 				}
  96. 		DFS(1,0,0);
  97. 		printf("%d",Ans);
  98. 	}
  99. }
  100. bool BFS()
  101. {
  102. 	for(int i=src; i<des; ++i) lev[i]=des+1;
  103. 	lev[des]=0, que[0]=des; int h=0,t=1;
  104. 	while(h<t)
  105. 	{
  106. 		int x=que[h++];
  107. 		for(int i=H[x]; i; i=nxt[i])
  108. 			if(lev[to[i]]==des+1 && cap[i^1])
  109. 				lev[to[i]]=lev[x]+1, que[t++]=to[i];
  110. 	}
  111. 	return lev[src]<=des;
  112. }
  113. void Augment(){
  114. 	for(int i=des; i!=src; i=fr[pre[i]])
  115. 		--cap[pre[i]], ++cap[pre[i]^1];
  116. }
  117. int ISAP()
  118. {
  119. 	if(!BFS()) return 0;
  120. 	for(int i=src; i<=des; ++i) ++num[lev[i]],cur[i]=H[i];
  121. 	int x=src,res=0; bool f=0;
  122. 	while(lev[src]<=des)
  123. 	{
  124. 		if(x==des) x=src,f=0,++res,Augment();
  125. 		bool can=0;
  126. 		for(int i=cur[x]; i; i=nxt[i])
  127. 			if((to[i]==des&&f)||(to[i]!=des &&  cap[i]))
  128. 			{
  129. 				printf("%d->%d lev:%d f:%d ",x,to[i],lev[x],f);
  130. 				if(x>n&&to[i]<=n) f=1;
  131. 				can=1, cur[x]=i, pre[x=to[i]]=i;
  132. 				printf("now:%d\n",f);
  133. 				break;
  134. 			}
  135. 		if(!can)
  136. 		{
  137. 			int mn=des;
  138. 			for(int i=H[x]; i; i=nxt[i])
  139. 				if((to[i]==des&&f)||(to[i]!=des&&cap[i])) mn=std::min(mn,lev[to[i]]);
  140. 			if(!--num[lev[x]]) break;
  141. 			printf("Change:%d %d->",x,lev[x]);
  142. 			++num[lev[x]=mn+1];
  143. 			printf("%d\n",lev[x]);
  144. 			cur[x]=H[x];
  145. 			if(x!=src)
  146. 			{
  147. 				if(x<=n&&fr[pre[x]]>n) f=0;
  148. 				x=fr[pre[x]];
  149. 			}
  150. 		}
  151. 	}
  152. 	return res;
  153. }
  155. int main()
  156. {
  157. 	freopen("pair.in","r",stdin);
  158. 	freopen("pair.out","w",stdout);
  160. 	n=read();
  161. 	if(n<=10) {Subtask1::Solve(); return 0;}
  162. 	Enum=1,src=0,des=n<<1|1;
  163. 	for(int i=1; i<=n; ++i) A[i]=read();
  164. 	for(int i=1; i<=n; ++i) B[i]=read();
  165. 	for(int i=1; i<=n; ++i) C[i]=read();
  166. 	for(int i=1; i<=n; ++i) AddEdge(src,i,INF),AddEdge(i+n,des,INF),AddEdge(i,i+n,B[i]);
  167. 	for(int i=1; i<=n; ++i)
  168. 		for(int j=1; j<=n; ++j)
  169. 			if(!(A[i]%A[j])&&i!=j)
  170. 			{
  171. 				int t=A[i]/A[j];
  172. 				if(is_p[0].count(t)) continue;
  173. 				if(is_p[1].count(t)) AddEdge(i+n,j,INF);
  174. 				else{
  175. 					bool f=Miller_Rabin(t);
  176. 					if(f) is_p[1].insert(t),AddEdge(i+n,j,INF);
  177. 					else is_p[0].insert(t);
  178. 				}
  179. 			}
  180. 	printf("%d",ISAP());
  182. 	return 0;
  183. }

丧心病狂的费用流(其实就是sb了)（18.4.1）

  1. #include <set>
  2. #include <queue>
  3. #include <cmath>
  4. #include <cstdio>
  5. #include <cctype>
  6. #include <algorithm>
  7. #define gc() getchar()
  8. typedef long long LL;
  9. const int N=410,M=100005,INF=1e9;
  10. const int P[9]={2,3,5,7,13,17,31};
  12. int n,range,src,des,A[N],B[N],C[N],Enum,H[N],nxt[M],fr[M],to[M],cap[M],Cap[M],pre[N];
  13. double dis[N],cost[M],neg[M];
  14. bool inq[N],ignore[M];
  15. std::queue<int> q;
  16. std::set<int> is_p[2];
  18. inline int read()
  19. {
  20. 	int now=0,f=1;register char c=gc();
  21. 	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
  22. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  23. 	return now*f;
  24. }
  25. inline void AddEdge(int u,int v,int w,int c)
  26. {
  27. 	if(!c){
  28. 		to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, Cap[Enum]=cap[Enum]=w, cost[Enum]=0, ignore[Enum]=1;
  29. 		to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, Cap[Enum]=cap[Enum]=0, cost[Enum]=0, ignore[Enum]=1;
  30. 	}
  31. 	else if(c<0){
  32. 		to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, Cap[Enum]=cap[Enum]=w, cost[Enum]=-log2((double)(-c)), neg[Enum]=1;
  33. 		to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, Cap[Enum]=cap[Enum]=0, cost[Enum]=log2((double)(-c)), neg[Enum]=1;
  34. 	}
  35. 	else{
  36. 		to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, Cap[Enum]=cap[Enum]=w, cost[Enum]=log2((double)c);
  37. 		to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, Cap[Enum]=cap[Enum]=0, cost[Enum]=-log2((double)c);
  38. 	}
  39. 	printf("%d->%d %d %d %.3lf\n",u,v,w,c,cost[Enum-1]);
  40. }
  41. int FP(LL x,int k,int p)
  42. {
  43. 	LL t=1;
  44. 	for(; k; k>>=1,x=x*x%p)
  45. 		if(k&1) t=t*x%p;
  46. 	return (int)t;
  47. }
  48. bool Miller_Rabin(int x)
  49. {
  50. 	if(x==2) return 1;
  51. 	if(!(x&1)||x==1) return 0;
  52. 	for(int i=0; i<7; ++i)
  53. 		if(x==P[i]) return 1;
  54. 		else if(!(x%P[i])) return 0;
  55. 	int u=x-1,t=0;
  56. 	while(!(u&1)) u>>=1,++t;
  57. 	for(int i=0; i<7; ++i)
  58. 	{
  59. 		LL now=FP(P[i],u,x),las;
  60. 		for(int j=1; j<=t; ++j)
  61. 		{
  62. 			las=now, (now*=now)%=x;
  63. 			if(now==1&&las!=1&&las!=x-1) return 0;
  64. 		}
  65. 		if(now!=1) return 0;
  66. 	}
  67. 	return 1;
  68. }
  69. bool SPFA()
  70. {
  71. 	for(int i=1; i<=des; ++i) dis[i]=1e14;
  72. 	dis[src]=0, q.push(src);
  73. 	while(!q.empty())
  74. 	{
  75. 		int x=q.front(); q.pop();
  76. 		inq[x]=0;
  77. 		for(int v,i=H[x]; i; i=nxt[i])
  78. 			if(cap[i]&&(neg[i]||dis[to[i]]>dis[x]+cost[i]))
  79. 			{
  80. 				dis[v=to[i]]=dis[x]+cost[i], pre[v]=i;
  81. 				if(!inq[v]) q.push(v);
  82. 			}
  83. 	}
  84. 	return dis[des]<1e14;
  85. }
  86. bool debug=1;
  87. double MCMF(LL &flow)
  88. {
  89. 	int mn=1e9; double res=0,tmp1=0;
  90. 	bool f=0;
  91. 	for(int i=des; i!=src; i=fr[pre[i]])
  92. 	{
  93. 		mn=std::min(cap[pre[i]],mn);
  94. 		if(ignore[i]&&i<=range) f=1;
  95. 	}
  96. //	printf("flow:%I64d ",flow);
  97. 	if(mn>=flow) mn=flow,flow=0;
  98. 	else flow-=mn;
  99. 	if(f){
  100. 		for(int i=des; i!=src; i=fr[pre[i]])
  101. 		{
  102. 			if(debug) printf("%d<-",i),
  103. 			cap[pre[i]]-=mn, cap[pre[i]^1]+=mn;
  104. 			if(neg[i]) ;
  105. 		}
  106. 		if(debug) printf("%d\n",src);
  107. 		if(debug) printf("then:%I64d %d\n",flow,mn);
  108. 		return 0;
  109. 	}
  110. 	for(int i=des; i!=src; i=fr[pre[i]])
  111. 	{
  112. 		if(debug) printf("%d<-",i),
  113. 		cap[pre[i]]-=mn, cap[pre[i]^1]+=mn, res+=cost[pre[i]];
  114. 		if(neg[pre[i]]&&!cost[pre[i]]) tmp1+=1.0;
  115. 	}
  116. 	if(debug) printf("%d\n",src);
  117. 	if(debug) printf("then:%I64d res:%.3lf 2:%.3lf\n",flow,res,res<0 ? -1.0*mn*pow(2.0,-res) : 1.0*mn*pow(2.0,res));
  118. 	return res<0 ? -1.0*mn*pow(2.0,-res)-tmp1 : 1.0*mn*pow(2.0,res)-tmp1;
  119. }
  120. bool Check(LL flow)
  121. {
  122. 	double sigma=0;
  123. 	LL tmp=flow;
  124. 	for(int i=2; i<=Enum; ++i) cap[i]=Cap[i];
  125. 	while(SPFA())
  126. 	{
  127. 		sigma+=MCMF(flow);
  128. 		if(!flow||sigma>0) break;
  129. 	}
  130. 	if(debug) printf("End:flow:%I64d bef:%I64d sigma:%.3lf\n\n",flow,tmp,sigma);
  131. 	return !flow&&sigma<=0;
  132. }
  134. int main()
  135. {
  136. 	freopen("pair.in","r",stdin);
  137. //	freopen("pair.out","w",stdout);
  139. 	n=read(),Enum=1,src=0,des=n<<1|1;
  140. 	for(int i=1; i<=n; ++i) A[i]=read();
  141. 	for(int i=1; i<=n; ++i) B[i]=read();
  142. 	for(int i=1; i<=n; ++i) C[i]=read();
  143. 	for(int i=1; i<=n; ++i) AddEdge(src,i,INF,0),AddEdge(i+n,des,INF,0),AddEdge(i,i+n,B[i],-C[i]);
  144. //	range=Enum;
  145. 	for(int i=1; i<=n; ++i)
  146. 		for(int j=1; j<=n; ++j)
  147. 			if(!(A[i]%A[j])&&i!=j)
  148. 			{
  149. 				int t=A[i]/A[j];
  150. 				if(is_p[0].count(t)) continue;
  151. 				if(is_p[1].count(t)) AddEdge(i+n,j,INF,0);
  152. 				else{
  153. 					bool f=Miller_Rabin(t);
  154. 					if(f) is_p[1].insert(t),AddEdge(i+n,j,INF,0);
  155. 					else is_p[0].insert(t);
  156. 				}
  157. 			}
  159. //	long long l=0,r=1e10,mid,ans=0;
  160. 	long long l=1,r=15,mid,ans=0;
  161. 	while(l<=r)
  162. 		if(Check(mid=l+r>>1)) ans=mid,l=mid+1;
  163. 		else r=mid-1;
  164. 	printf("%I64d",ans);
  166. 	return 0;
  167. }

BZOJ.4514.[SDOI2016]数字配对(费用流SPFA 二分图)的更多相关文章

  1. BZOJ 4514: [Sdoi2016]数字配对 [费用流 数论]

    4514: [Sdoi2016]数字配对 题意: 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ai/aj 是一个质数 ...

  2. 4514: [Sdoi2016]数字配对 费用流

    链接 https://www.lydsy.com/JudgeOnline/problem.php?id=4514 思路 EK直接贪心做 <0的时候加上剩余返回 二分图a->b的时候 把b- ...

  3. 图论（费用流）：BZOJ 4514 [Sdoi2016]数字配对

    4514: [Sdoi2016]数字配对 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 820  Solved: 345[Submit][Status ...

  4. BZOJ 4514: [Sdoi2016]数字配对

    4514: [Sdoi2016]数字配对 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1606  Solved: 608[Submit][Statu ...

  5. 【BZOJ4514】[Sdoi2016]数字配对 费用流

    [BZOJ4514][Sdoi2016]数字配对 Description 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ...

  6. bzoj4514: [Sdoi2016]数字配对--费用流

    看了一眼题目&数据范围,觉得应该是带下界的费用流 原来想拆点变成二分图,能配对的连边,跑二分图,可行性未知 后来看到另外一种解法.. 符合匹配要求的数要满足:质因子的个数相差为1,且两者可整除 ...

  7. 【BZOJ 4514】[Sdoi2016]数字配对 费用流

    利用spfa流的性质,我直接拆两半,正解分奇偶(妙),而且判断是否整除且质数我用的是暴力根号,整洁判断质数个数差一(其他非spfa流怎么做?) #include <cstdio> #inc ...

  8. 【BZOJ4514】【SDOI2016】数字配对 [费用流]

    数字配对 Time Limit: 10 Sec  Memory Limit: 128 MB[Submit][Status][Discuss] Description 有 n 种数字,第 i 种数字是 ...

  9. [bzoj4514]数字配对[费用流]

    今年SDOI的题,看到他们在做,看到过了一百多个人,然后就被虐惨啦... 果然考试的时候还是打不了高端算法,调了...几天 默默地yy了一个费用流构图: 源连所有点,配对的点连啊,所有点连汇... 后 ...

随机推荐

  1. Here’s just a fraction of what you can do with linear algebra

    Here’s just a fraction of what you can do with linear algebra The next time someone wonders what the ...

  2. javascript命名空间

    命名空间:命名空间有助于减少程序中所需要的全局变量的数量,并且同时还有助于避免命名冲突或过长的名字前缀. 板栗: var MYAPP = MYAPP || {}; MYAPP.namespace = ...

  3. [百度地图] 用于类似 DWZ UI 框架的 百度地图 功能封装类 [MultiZMap.js] 实例源码

    MultiZMap 功能说明 MultiZMap.js 本类方法功能大多使用 prototype 原型 实现,它是 ZMap 的多加载版本,主要用于类似 DWZ 这个 多标签的 UI 的框架: 包含的 ...

  4. es6笔记（4） Set数据结构

    概要 介绍: 集合是由一组无序且唯一的项组成的,这个数据结构使用了与有限集合相同的数学概念,应用在计算机的数据结构中. ES6提供了数据结构Set.它类似于数组,但是没有重复的值. 特点: key与v ...

  5. C# 解决VS2008在win7找不到输入序列号的地方

    1.VS2008在Windows7 打开维护界面看不到可以输序列号的地方. 因为微软把他隐藏了. 2.我们可以借用工具把他显示出来 下载地址:http://www.zlsoft.com/techbbs ...

  6. python实现两个经纬度点之间的距离和方位角

    from:http://blog.csdn.net/zhuqiuhui/article/details/53180395 1.  求两个经纬点的方位角,P0(latA, lonA), P1(latB, ...

  7. MySQL多源复制【转】

    什么是多源复制? 首先,我们需要清楚 multi-master 与multi-source 复制不是一样的. Multi-Master 复制通常是环形复制, 你可以在任意主机上将数据复制给其他主机. ...

  8. how-to-pass-a-class-variable-to-a-decorator-inside-class-definition

    https://stackoverflow.com/questions/17522706/how-to-pass-a-class-variable-to-a-decorator-inside-clas ...

  9. Java工具库：

    1.  重试框架: https://docs.spring.io/spring-batch/trunk/reference/html/retry.html <dependency> < ...

  10. 重温CSS之文档结构

    我们来看看几个基本的HTML页面 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "h ...