POJ1061

王子野心 — Mon, 08 Aug 2011 15:43:00 GMT

对于给定的x,y,m,n,L;

当m=n的时候必定无解。

可以假设m>n;

那么对于答案a

有(am-x)-(an-y)=kL

可以化成a(m-n)+y=kL+x

也可化成a(m-n)-kL=y-x

显然，如果(y-x)%gcd(m-n,L)!=0,无解

令Q=a(m-n)+y=kL+x

则Q%(m-n)=y

Q%L=x

我们可以枚举a，求出Q，判断是否有Q%L=x

这种方法会TLE

更好的方法是枚举k，判断Q%(m-n)=y

需要注意的是，m-n有可能小于y

王子野心 2011-08-08 23:43 发表评论

求网络最大流的Dinic模板，过了几个题，但还不敢保证完全正确

王子野心 — Thu, 28 Jul 2011 09:35:00 GMT

1 #include<iostream>
2 #include<cstdio>
3 #include<cstring>
4 #define min(x,y) x
5 using namespace std;
6 const int N=1100;
7 const int M=11000;
8 const int INF=100000000;
9 int g[N],dis[N],q[N];
10 int en[M],val[M],next[M],rev[M];
11 int n,m,S,T,tot,ans;
12
13 int addedge(int a,int b,int c){
14     tot++;en[tot]=b;val[tot]=c;next[tot]=g[a];g[a]=tot;
15     tot++;en[tot]=a;val[tot]=0;next[tot]=g[b];g[b]=tot;
16     rev[g[a]]=g[b];
17     rev[g[b]]=g[a];
18     return 0;
19 }
20
21 int init(){
22     int i,a,b,c;
23     if(scanf("%d%d",&n,&m)==EOF)
24         return 0;
25     S=1;T=n;
26     for(i=1;i<=n;i++)
27         g[i]=0;
28     for(i=1;i<=m;i++){
29         scanf("%d%d%d",&a,&b,&c);
30         addedge(a,b,c);
31     }
32     return 1;
33 }
34
35 int build(){
36     int hp=0,tp=0;
37     memset(dis,0,sizeof(dis));
38     q[0]=S;dis[S]=1;
39     for(hp=0;hp<=tp;hp++){
40         int k=q[hp];
41         for(int p=g[k];p;p=next[p]){
42             if(val[p]>0 && dis[en[p]]==0){
43                 q[++tp]=en[p];
44                 dis[q[tp]]=dis[k]+1;
45                 if(q[tp]==T)return 1;
46             }
47         }
48     }
49     return 0;
50 }
51
52 int  find(int s,int aug){
53     int flow=0;
54     if(s==T || aug==0) return aug;
55     for(int p=g[s];p;p=next[p]){
56         if(val[p]>0 && dis[en[p]]-1==dis[s]){
57             int k=find(en[p],min(val[p],aug));
58             flow+=k;
59             aug-=k;
60             val[p]-=k;
61             val[rev[p]]+=k;
62             if(aug==0) return flow;
63         }
64     }
65     return flow;
66 }
67
68 int work(){
69     while(build()){
70         ans+=find(S,INF);
71     }
72     printf("%d\n",ans);
73     return 0;
74 }
75
76 int main(){
77     while(init())
78         work();
79     return 0;
80 }
81

王子野心 2011-07-28 17:35 发表评论

C++博客-半壁江山の圣灵玉-随笔分类-龙神手札（一）

POJ1061

求网络最大流的Dinic模板，过了几个题，但还不敢保证完全正确