codevs1409 拦截导弹2

Posted 2020-08-11 №〓→龙光←

【问题描述】
一场战争正在 A 国与 B 国之间如火如荼的展开。
B 国凭借其强大的经济实力开发出了无数的远程攻击导弹,B 国的领导人希望,通过这些导弹直接毁灭 A 国的指挥部,从而取得战斗的胜利!当然,A 国人民不会允许这样的事情发生,所以这个世界上还存在拦截导弹。
现在,你是一名 A 国负责导弹拦截的高级助理。B 国的导弹有效的形成了三维立体打击,我们可以将这些导弹的位置抽象三维中间的点(大小忽略),为了简单起见,我们只考虑一个瞬时的状态,即他们静止的状态。
拦截导弹设计非常精良,可以精准的引爆对方导弹而不需要自身损失,但是 A 国面临的一个技术难题是,这些导弹只懂得直线上升。精确的说,这里的直线上升指 xyz 三维坐标单调上升。给所有的 B 国导弹按照 1 至 N 标号,一枚拦截导弹可以打击的对象可以用一个 xyz 严格单调上升的序列来表示,例如:B 国导弹位置:(0, 0, 0) (1, 1, 0) (1, 1, 1), (2, 2, 2)一个合法的打击序列为:{1, 3, 4}一个不合法的打击序列为{1, 2, 4}
A 国领导人将一份导弹位置的清单交给你,并且向你提出了两个最简单不过的问题(假装它最简单吧):
1.一枚拦截导弹最多可以摧毁多少 B 国的导弹?
2.最少使用多少拦截导弹才能摧毁 B 国的所有导弹?
不管是为了个人荣誉还是国家容易,更多的是为了饭碗,你,都应该好好的把这个问题解决掉!
【输入文件】
第一行一个整数 N 给出 B 国导弹的数目。
接下来 N 行每行三个非负整数 Xi, Yi, Zi 给出一个导弹的位置,你可以假定任意两个导弹不会出现在同一位置。
【输出文件】
输出文件有且仅有两行。
第一行输出一个整数 P,表示一枚拦截导弹之多能够摧毁的导弹数。
第二行输出一个整数 Q,表示至少需要的拦截导弹数目。
【输入输出样例】
输入
4
0 0 0
1 1 0
1 1 1
2 2 2
输出
3
2
【数据范围】
所有的坐标都是[0,10 6 ]的整数
对于 30%的数据满足 N < 31
对于 50%的数据满足 N < 101
对于 100%的数据满足 N < 1001

---

 

我们按照先后顺序连出一张DAG，第一问求的就是最长链，DP，搜都没有什么问题。第二问问要多少枚导弹，这不就是问这张DAG需要多少条路才能被完全覆盖，问题转换为可以走重复路径最小路径覆盖问题，网络流或匈牙利算法解决（若果你不会最小路径覆盖：http://www.cnblogs.com/Dragon-Light/p/5604865.html）

---

 

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstdlib>
 5 #include<algorithm>
 6 #include<vector>
 7 #include<cmath>
 8 #include<ctime>
 9 #include<cstring>
10 #define llg int
11 #define maxn 3010
12 #define md 1001000
13 using namespace std;
14 int dl[md+10],n,m,a[maxn][maxn],bj[maxn],ans=0,x,used[maxn],girl[maxn],y,i,j,k,f[maxn][maxn],val[maxn],deep[maxn],head,tail;
15 vector<llg>g[maxn];
16 struct node
17 {
18     llg a1,a2,a3;
19 }c[maxn];
20     
21 bool find(llg x)
22 {
23     for (llg i=1;i<=n*2;i++)
24     {
25         if (a[x][i] && !used[i])
26         {
27             used[i]=1;
28             if (!girl[i] || find(girl[i]))
29             {
30                 girl[i]=x;
31                 return true;
32             }
33         }
34     }
35     return false;
36 }
37 
38 llg BFS(llg x)
39 {
40     for (llg i=1;i<=n;i++) bj[i]=0;
41     dl[1]=x,head=0,tail=1; llg maxl=1,w; deep[x]=1;
42     do
43     {
44         head=(head % md)+1;
45         x=dl[head]; bj[x]=0;
46         maxl=max(maxl,deep[x]);
47         w=g[x].size();
48         for (llg i=0;i<w;i++)
49             if (deep[g[x][i]]<deep[x]+1)
50             {
51                 deep[g[x][i]]=deep[x]+1;
52                 if (!bj[g[x][i]])
53                 {
54                     deep[g[x][i]]=deep[x]+1;
55                     tail=(tail % md)+1;
56                     dl[tail]=g[x][i];
57                 }
58             }
59     }while (head!=tail);
60     return maxl;
61 }
62 
63 int main()
64 {
65     //yyj("bomb");
66     cin>>n;
67     for (i=1;i<=n;i++) scanf("%d%d%d",&c[i].a1,&c[i].a2,&c[i].a3);
68     for (i=1;i<=n;i++)
69         for (j=1;j<=n;j++)
70             if (c[i].a1<c[j].a1 && c[i].a2<c[j].a2 && c[i].a3<c[j].a3)
71             {
72                 f[i][j]=1;
73                 g[i].push_back(j);
74                 val[j]++;
75             }
76     for (i=1;i<=n;i++)
77         if (val[i]==0) ans=max(BFS(i),ans);
78     cout<<ans<<endl;
79 /*    for (k=1;k<=n;k++)
80         for (i=1;i<=n;i++)
81             for (j=1;j<=n;j++)
82                 f[i][j]=f[i][j] || (f[i][k] && f[k][j]);*/
83     for (i=1;i<=n;i++)
84         for (j=1;j<=n;j++)
85             if (i!=j && f[i][j]) a[i][n+j]=1;
86     ans=0;
87     for (i=1;i<=n*2;i++)
88     {
89         memset(used,0,sizeof(used));
90         if (find(i)) ans++;
91     }
92    cout<<n-ans;
93     return 0;
94 }