JLU数据结构上机实验5

7-1 图的深度优先搜索I

原题呈现：

解题思路 ： targin算法的原理部分+时间戳+保存搜索树中的边

以下是根据这一思路写出的代码：

#include<bits/stdc++.h>
using namespace std;

int dfn[100005], tim, ed[100005];
int n, m, cnt, ans;
vector<int>e[100005];

struct edge{
    int u, v;  
    edge(int a, int b) :u(a), v(b) {}
};
bool cmp(const edge& a, const edge& b) {
    return a.v < b.v;
}
vector<edge> res;

void targin(int x)
{
    dfn[x] = ++tim;
    for (int u = 0; u < e[x].size(); ++u){
        if (!dfn[e[x][u]]){
            targin(e[x][u]);
            res.push_back(edge(x, e[x][u]));
        }
    }
    ed[x] = ++tim;
}

int main()
{
    scanf("%d %d", &n, &m);
    int u, v;
    for (int i = 1; i <= m; ++i){
        scanf("%d %d", &u, &v);
        e[u].push_back(v);
        e[v].push_back(u);
    }
    for (int i = 1; i <= n; ++i){
        if (!dfn[i])
            targin(i);
    }
    for (int i = 1; i <= n; ++i){
        printf("%d %d\n", dfn[i], ed[i]);
    }
    printf("%d\n", res.size());
    sort(res.begin(), res.end(), cmp);
    for (int i = 0; i < res.size(); ++i){
        printf("%d %d\n", res[i].u, res[i].v);
    }
}

解后反思： 通过每条边上两端点时间戳的大小比较可以得到边的相关性质。

7-2 圆

原题呈现：

解题思路 ： 考虑从一个结点开始用bfs，直到队空证明一个块中的全部圆已被找到；找到下一个没有入过队的结点入队，继续bfs…重复以上过程，直到所有节点都入过队。其中数组pd表示一个数组是否入过队。

以下是根据这一思路写出的代码：

#include<bits/stdc++.h>
using namespace std;

long long x[8010][3], res=0, pd[8010];
queue<int> q;

int main() {
	int n;
	scanf("%d", &n);
	for (int i = 0; i < n; i++)
		scanf("%lld %lld %lld", &x[i][0], &x[i][1], &x[i][2]);

    q.push(0);
	while(!q.empty()){	
        int tmp = q.front();
		q.pop();
		pd[tmp] = 1;
		for (int j = 0; j < n; j++) {
;			if (!pd[j] && ((x[tmp][0] - x[j][0]) * (x[tmp][0] - x[j][0]) + (x[tmp][1] - x[j][1]) * (x[tmp][1] - x[j][1]) <= (x[tmp][2] + x[j][2]) * (x[tmp][2] + x[j][2]))) {
				pd[j] = 1;
                q.push(j);
			}
		}
		if (q.empty()) {
			res++;
            for (int j = 0; j < n; j++) {
			    if (!pd[j]) {
				    q.push(j);
				    break;
		    	}
	     	}
		}		
	}
	printf("%lld", res);
}

解后反思： 1.判断两圆是否属于一个块时用距离的平方作比较，不要开根以避免精度损失。 2.x数组开成long long，避免计算距离平方时超出int精度范围。 3.除该种做法以外，还可以用并查集的思想。

7-3 供电

原题呈现：

解题思路： 该题与平时所见的用到最小生成树的题目唯一的不同之处在于：可以在某点处付出一定代价，使之不必再与其它点相连接。因而考虑给dist数组赋初值为每点的cost。此时面对另一个问题：至少需要一个点建供电站，并与其他未建供电站的点相连接，因而选出cost最小的点，从该点出发进行prim算法即可。

以下是根据这一思路写出的代码：

#include<bits/stdc++.h>
using namespace std;
#define maxn 10005
#define inf 1000001

int vis[maxn], dist[maxn],cnt = 1;
vector<pair<int, int> >v[maxn];
long long ans = 0;

int main() {
	int n, m;
	scanf("%d %d", &n, &m);
	int mincost = inf;
	int startnode = 0;
	for (int i = 1; i <= n; i++) {
		scanf("%d", &dist[i]);
	}
	for (int i = 1; i <= n; i++) {
		if (mincost > dist[i]) {
			mincost = dist[i];
			startnode = i;
		}
	}
    vis[startnode] = 1;
	ans += mincost;
	for (int i = 1; i <= m; i++) {//邻接链表存储图
		int node1, node2, len;
		scanf("%d %d %d", &node1, &node2, &len);
		v[node1].push_back(make_pair(node2, len));
		v[node2].push_back(make_pair(node1, len));
	}

	int temp = startnode;
	//下面正式开始prime算法
	while (cnt < n) {
		for (int i = 0; i < v[temp].size(); i++) {
			if (v[temp][i].second < dist[v[temp][i].first])
				dist[v[temp][i].first] = v[temp][i].second;//松弛
		}
		int min = inf;
		int midnode = 0;
		for (int j = 1; j <= n; j++) {
			if (vis[j] == 1)
				continue;
			if (min > dist[j]){
				min = dist[j];
				midnode = j; 
			}
		}
		ans += min;
		cnt++;
		vis[midnode] = 1;
		temp = midnode;
	}
	printf("%lld\n", ans);
	return 0;
}

解后反思： 这种方法用时较长，因而可以考虑以下两种改进方式：1.堆优化 2.尝试kruskal算法

7-4 发红包

原题呈现：

解题思路： 由于要花钱最小，因而对没有要求的人只发下限888RMB，对有要求的只比他的比较对象多发1RMB，故把要钱人少的人发完后才轮到要钱人多的。由此想到拓扑排序。若形成环路，证明不能满足所有员工要求，输出-1。

以下是根据这一思路写出的代码：

#include<bits/stdc++.h>
using namespace std;

vector<long> v[20010];
long _in[20010]; // 记录index位置入边的数量
long res[20010]; // 记录的是相对888的数值大小

int main() {
	long n = 0, m = 0, from = 0, to = 0, cnt = 0;
	scanf("%d %d", &n, &m);

	for (long i = 0; i < m; i++) {
		scanf("%d %d", &from, &to);
		_in[from]++;
		v[to].push_back(from);
	}
	// 找到入边为0的即可删除相应的边，并通过差分数组更新值
	for (long i = 0; i < n; i++) {
		long k = 0;
		// 遍历查找
		for (long j = 1; j <= n; j++) {
			if (_in[j] == 0) {
				cnt++;
				_in[j] = -1;
				k = j;
				break;
			}
		}

		long len = v[k].size();
		for (long j = 1; j <= len; j++) {
			long index = v[k][j - 1];
			_in[index]--;
			res[index] = max(res[index], res[k] + 1);
		}
	}
	// 所有的边都能消去，即无闭环
	if (cnt == n) {
		long sum = 0;
		for (long i = 1; i <= n; i++) {
			sum += res[i];
		}
		printf("%d", sum + 888 * n);
	}
	// 出现了闭环
	else {		
		printf("-1");
	}
}

解后反思： 出度为0的点代表没有要求，也就是可以直接确定发888RMB。因而该题与平时按入度进行拓扑排序的题目不同，要按出度进行拓扑排序。