problem

luogu-P2065

solution

這個(gè)拿走一組共兩張卡片的操作其實(shí)就是一個(gè)匹配。

直接兩個(gè)數(shù)的最大公約數(shù)大于 $1$ 就建一條邊，跑二分圖匹配最大流即可。

然而如果直接枚舉兩個(gè)數(shù)然后算他們的 $\text{gcd}$ ，時(shí)間復(fù)雜度 $O(T{n}^2\log V)$ 會(huì) $\text{TLE}$。

又不能預(yù)處理任意兩個(gè)數(shù)的 $\text{gcd}$，這個(gè)數(shù)字 $\in (1,1e7)$ 實(shí)在是太大了。

這里很巧妙，任何數(shù)都可以拆成若干個(gè)質(zhì)數(shù)的冪，兩個(gè)數(shù)的 $\text{gcd}>1$ 無非就是兩個(gè)數(shù)有相同的某些個(gè)質(zhì)因子。

所以我們可以預(yù)處理出 $1e7$ 以內(nèi)的所有質(zhì)數(shù)，這并不多，然后將一個(gè)數(shù)和其質(zhì)因子連邊。

只顯然每張卡片只能拿走一次，流量與源匯點(diǎn)之間為 $1$ 即可。

這樣時(shí)間復(fù)雜度就是 $O(Tn\sqrt{n})$ ，網(wǎng)絡(luò)流就沒考慮反正是 $O($ 能過 $)$。

code

#include <bits/stdc++.h> using namespace std; #define maxn 100005 #define maxm 800005 int T, m, n, s, t, cnt = -1, tot; struct node { int to, nxt, flow; }E[maxm]; int b[maxn], r[maxn], head[maxn], cur[maxn], dep[maxn]; queue < int > q;int gcd( int x, int y ) {if( ! y ) return x;else return gcd( y, x % y ); }void addedge( int u, int v, int w ) {E[++ cnt] = { v, head[u], w }, head[u] = cnt;E[++ cnt] = { u, head[v], 0 }, head[v] = cnt; }bool bfs() {memset( dep, 0, sizeof( dep ) );memcpy( cur, head, sizeof( head ) );dep[s] = 1; q.push( s );while( ! q.empty() ) {int u = q.front(); q.pop();for( int i = head[u];~ i;i = E[i].nxt ) {int v = E[i].to;if( ! dep[v] and E[i].flow ) {dep[v] = dep[u] + 1;q.push( v );}}}return dep[t]; }int dfs( int u, int cap ) {if( u == t or ! cap ) return cap;int flow = 0;for( int i = cur[u];~ i;i = E[i].nxt ) {cur[u] = i; int v = E[i].to;if( dep[v] == dep[u] + 1 ) {int w = dfs( v, min( cap, E[i].flow ) );if( ! w ) continue;E[i ^ 1].flow += w;E[i].flow -= w;flow += w;cap -= w;if( ! cap ) break;}}return flow; }int dinic() {int ans = 0;while( bfs() ) ans += dfs( s, 1e9 );return ans; }#define maxp 10000005 int prime[maxp], vis[maxp], num[maxp]; int cntp; void sieve() {for( int i = 2;i < 1e7;i ++ ) {if( ! vis[i] ) prime[++ cntp] = i;for( int j = 1;j <= cntp and i * prime[j] < 1e7;j ++ ) {vis[i * prime[j]] = 1;if( i % prime[j] == 0 ) break;}} } void divide( int id, int val ) {for( int i = 1;i <= cntp and prime[i] <= val;i ++ ) if( val % prime[i] == 0 ) {if( ! num[prime[i]] ) num[prime[i]] = ++ tot;if( id <= m ) addedge( id, num[prime[i]], 1e9 );else addedge( num[prime[i]], id, 1e9 );while( val % prime[i] == 0 ) val /= prime[i];} }int main() {sieve(); tot = 1001;scanf( "%d", &T );while( T -- ) {memset( head, -1, sizeof( head ) ); cnt = -1;scanf( "%d %d", &m, &n );s = 0, t = n + m + 1;for( int i = 1;i <= m;i ++ ) scanf( "%d", &b[i] );for( int i = 1;i <= n;i ++ ) scanf( "%d", &r[i] );for( int i = 1;i <= m;i ++ ) divide( i, b[i] );for( int i = 1;i <= n;i ++ ) divide( i + m, r[i] );for( int i = 1;i <= m;i ++ ) addedge( s, i, 1 );for( int i = 1;i <= n;i ++ ) addedge( i + m, t, 1 );printf( "%d\n", dinic() );}return 0; }

總結(jié)

以上是生活随笔為你收集整理的[TJOI2011] 卡片（网络流 + 质因子优化建图）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。