【BZOJ】【3275】Numbers

网络流/最小割

---

　　Orz了Jiry_2神犇，蒟蒻网络流建模什么的完全不会啊T_T

　　按奇偶性来分组实在太巧妙了……然后相关的点之间连边表示只能选其一，来求最小割……

 /**************************************************************
     Problem: 3275
     User: Tunix
     Language: C++
     Result: Accepted
     Time:1252 ms
     Memory:1708 kb
 ****************************************************************/

 //BZOJ 3275
 #include<cmath>
 #include<vector>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 #define pb push_back
 using namespace std;
 inline int getint(){
     int v=,sign=; char ch=getchar();
     while(ch<''||ch>''){ if (ch=='-') sign=-; ch=getchar();}
     while(ch>=''&&ch<=''){ v=v*+ch-''; ch=getchar();}
     return v*sign;
 }
 const int N=,M=,INF=~0u>>;
 typedef long long LL;
 /******************tamplate*********************/
 int n,m,a[],tot,ans;
 struct edge{
     int to,v;
 };
 int gcd(int a,int b){if (!b) return a;else return gcd(b,a%b);}
 bool judge(int a,int b){
     int s=a*a+b*b;
     int q=int(sqrt(s));
     if (q*q!=s) return ;
     if (gcd(a,b)!=) return ;
     return ;
 }
 struct Net{
     edge E[M];
     int head[N],next[M],cnt;
     bool vis[][];
     void ins(int x,int y,int v){
         E[++cnt]=(edge){y,v};
         next[cnt]=head[x]; head[x]=cnt;
     }
     void add(int x,int y,int v){
         ins(x,y,v); ins(y,x,);
     }
     int s,t,cur[N],d[N],Q[N];
     void init(){
         n=getint();cnt=;
         tot=ans=;
         s=; t=n+;
         F(i,,n){
             a[i]=getint();tot+=a[i];
             if (a[i]&) add(s,i,a[i]);
             else add(i,t,a[i]);
         }
         F(i,,n) if (a[i]&)
             F(j,,n) if ((a[j]&)==)
                 if (judge(a[i],a[j])) add(i,j,INF);
     }
     bool mklevel(){
         memset(d,-,sizeof d);
         d[s]=;
         int l=,r=-;
         Q[++r]=s;
         while(l<=r){
             int x=Q[l++];
             for(int i=head[x];i;i=next[i])
                 if (d[E[i].to]==- && E[i].v){
                     d[E[i].to]=d[x]+;
                     Q[++r]=E[i].to;
                 }
         }
         return d[t]!=-;
     }
     int dfs(int x,int a){
         if (x==t) return a;
         int flow=;
         for(int &i=cur[x];i && flow<a;i=next[i])
             if (E[i].v && d[E[i].to]==d[x]+){
                 int f=dfs(E[i].to,min(a-flow,E[i].v));
                 E[i].v-=f;
                 E[i^].v+=f;
                 flow+=f;
             }
         if (!flow) d[x]=-;
         return flow;
     }
     void Dinic(){
         while(mklevel()){
             F(i,s,t) cur[i]=head[i];
             ans+=dfs(s,INF);
         }
     }
 }G1;
 int main(){
 #ifndef ONLINE_JUDGE
     freopen("3275.in","r",stdin);
     freopen("3275.out","w",stdout);
 #endif
     G1.init(); G1.Dinic();
     printf("%d\n",tot-ans);
     return ;
 }

3275: Number

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 495  Solved: 217
[Submit][Status][Discuss]

Description

有N个正整数，需要从中选出一些数，使这些数的和最大。
若两个数a,b同时满足以下条件，则a,b不能同时被选
1:存在正整数C，使a*a+b*b=c*c
2:gcd(a,b)=1

Input

第一行一个正整数n，表示数的个数。

第二行n个正整数a1,a2,?an。

 

 

Output

最大的和。

 

Sample Input

5
3 4 5 6 7

Sample Output

22

HINT

n<=3000。

Source

网络流

[Submit][Status][Discuss]