SRM 565

最新推荐文章于 2019-01-08 08:41:16 发布
最新推荐文章于 2019-01-08 08:41:16 发布
阅读量797 收藏
点赞数
分类专栏： topcoder 文章标签： topcoder SRM
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/merlininice/article/details/18719819
版权

好多天没有进度了。。就因为这个。。这个丧心病狂的题。。丧心病狂的lyrically

分四种情况

1、三个相邻

2、两个相邻，和第三个不相邻

3、三个都不相邻并且成Y型

4、三个都不相邻并且成一条直线

然后每种情况是形成一个大的结构，讨论讨论，都能找到其他点在这个大的结构上或者从大的结构某个点伸出去，而且位置唯一确定。。。。剩下的就是码代码啊。。。。


我觉得如果你不是闲得无聊或者有自虐倾向的，远离此题QAQ（或许有简单的写法呢呆逼up主

#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <cstring>

using namespace std;

class UnknownTree {
public:
	int getCount(vector<int> , vector<int> , vector<int> );
};
long long pmod = 1000000009;
int cc[100];
bool cmp(int a, int b) {
	return cc[a] < cc[b];
}
long long f(vector<int> p, int i) {
	if (p.size() == i)
		return 1;
	sort(p.begin(), p.end(), cmp);
	int v;
	long long ret, now, up;
	ret = 1;
	now = 1;
	up = 1;
	v = cc[p[i]];
	for (; i < p.size(); ++i) {
		if (v != cc[p[i]]) {
			v = cc[p[i]];
			now = up;
		}
		ret *= now;
		ret %= pmod;
		up++;
	}
	return ret;
}
long long gao1(vector<int> a, vector<int> b, vector<int> c) {
	int i;
	int ca, cb, cc;
	int n = a.size();
	vector<int> va, vb, vc;
	int x = abs(a[0] - b[0]);
	int y = abs(c[0] - b[0]);
	ca = cb = cc = 0;
	for (i = 0; i < n; ++i) {
		if (abs(a[i] - b[i]) != x)
			return 0;
		if (abs(c[i] - b[i]) != y)
			return 0;
		if (a[i] < b[i] && b[i] < c[i])
			va.push_back(i);
		else if (a[i] > b[i] && b[i] > c[i])
			vc.push_back(i);
		else if (a[i] > b[i] && c[i] > b[i])
			vb.push_back(i);
		else
			return 0;
	}
	return f(va, 0) * f(vb, 0) % pmod * f(vc, 0) % pmod;
}
long long gao2(vector<int> a, vector<int> b, vector<int> c) {
	int i, j;
	int n = a.size();
	int lab = abs(a[0] - b[0]);
	if (lab == 0)
		return 0;
	for (i = 0; i < n; ++i) {
		if (abs(a[i] - b[i]) != lab)
			return 0;
	}
	int lbc = b[0] + c[0];
	for (i = 0; i < n; ++i) {
		if (lbc > b[i] + c[i])
			lbc = b[i] + c[i];
	}
	int ed = 1;
	map<int, int> m;
	vector<int> v[55];
	vector<int> va;
	for (i = 0; i < n; ++i) {
		if (a[i] < b[i]) {
			if (c[i] - b[i] != lbc)
				return 0;
			va.push_back(i);
		} else {
			if (abs(c[i] - b[i]) > lbc)
				return 0;
			int p = m[c[i] - b[i]];
			if (p == 0) {
				p = ed++;
				m[c[i] - b[i]] = p;
			}
			v[p].push_back(i);
		}
	}
	long long ret = f(va, 0);
	for (i = 1; i < ed; ++i) {
		int rr = 0;
		for (j = 0; j < v[i].size(); ++j) {
			if (b[v[i][j]] + c[v[i][j]] == lbc)
				rr++;
		}
		if (abs(c[v[i][0]] - b[v[i][0]]) != lbc && rr != 1)
			return 0;
		if (abs(c[v[i][0]] - b[v[i][0]]) != lbc)
			ret *= f(v[i], 1);
		else
			ret *= f(v[i], 0);
		ret %= pmod;
	}
	return ret;
}
long long gao3(vector<int> a, vector<int> b, vector<int> c) {
	int i, j;
	int n = a.size();
	int lab = a[0] + b[0];
	for (i = 0; i < n; ++i) {
		if (lab > a[i] + b[i])
			lab = a[i] + b[i];
	}
	int lac = a[0] + c[0];
	for (i = 0; i < n; ++i) {
		if (lac > a[i] + c[i])
			lac = a[i] + c[i];
	}
	int lbc = b[0] + c[0];
	for (i = 0; i < n; ++i) {
		if (lbc > b[i] + c[i])
			lbc = b[i] + c[i];
	}
	int mid = -1;
	for (i = 0; i < n; ++i) {
		if (a[i] + b[i] == lab && a[i] + c[i] == lac && b[i] + c[i] == lbc) {
			if (mid != -1) {
				return 0;
			}
			mid = i;
		}
	}
	if (mid == -1)
		return 0;

	int an = 1, bn = 1, cn = 1;
	map<int, int> am, bm, cm;
	vector<int> av[55], bv[55], cv[55], vo;
	for (i = 0; i < n; ++i) {
		if (i == mid)
			continue;
		if (b[i] - a[i] == b[mid] - a[mid] && c[i] - a[i] == c[mid] - a[mid]) {
			vo.push_back(i);
		} else if (b[i] - a[i] > b[mid] - a[mid] && c[i] - a[i] > c[mid]
				- a[mid]) {
			int pos = b[i] - a[i] - b[mid];
			if (pos != c[i] - a[i] - c[mid])
				return 0;
			int p = am[pos];
			if (p == 0) {
				p = an++;
				am[pos] = p;
			}
			av[p].push_back(i);
		} else if (a[i] - b[i] > a[mid] - b[mid] && c[i] - b[i] > c[mid]
				- b[mid]) {
			int pos = a[i] - b[i] - a[mid];
			if (pos != c[i] - b[i] - c[mid])
				return 0;
			int p = bm[pos];
			if (p == 0) {
				p = bn++;
				bm[pos] = p;
			}
			bv[p].push_back(i);
		} else if (b[i] - c[i] > b[mid] - c[mid] && a[i] - c[i] > a[mid]
				- c[mid]) {
			int pos = b[i] - c[i] - b[mid];
			if (pos != a[i] - c[i] - a[mid])
				return 0;
			int p = cm[pos];
			if (p == 0) {
				p = cn++;
				cm[pos] = p;
			}
			cv[p].push_back(i);
		} else {
			return 0;
		}
	}
	long long ret = f(vo, 0);
	for (i = 1; i < an; ++i) {
		int rr = 0;
		for (j = 0; j < av[i].size(); ++j) {
			if (a[av[i][j]] + b[av[i][j]] == lab)
				rr++;
		}
		if (abs(a[av[i][0]] - b[av[i][0]]) != lab && rr != 1)
			return 0;
		if (abs(a[av[i][0]] - b[av[i][0]]) != lab)
			ret *= f(av[i], 1);
		else
			ret *= f(av[i], 0);
		ret %= pmod;
	}
	for (i = 1; i < bn; ++i) {
		int rr = 0;
		for (j = 0; j < bv[i].size(); ++j) {
			if (a[bv[i][j]] + b[bv[i][j]] == lab)
				rr++;
		}
		if (abs(a[bv[i][0]] - b[bv[i][0]]) != lab && rr != 1)
			return 0;
		if (abs(a[bv[i][0]] - b[bv[i][0]]) != lab)
			ret *= f(bv[i], 1);
		else
			ret *= f(bv[i], 0);
		ret %= pmod;
	}
	for (i = 1; i < cn; ++i) {
		int rr = 0;
		for (j = 0; j < cv[i].size(); ++j) {
			if (c[cv[i][j]] + b[cv[i][j]] == lbc)
				rr++;
		}
		if (abs(c[cv[i][0]] - b[cv[i][0]]) != lbc && rr != 1)
			return 0;
		if (abs(c[cv[i][0]] - b[cv[i][0]]) != lbc)
			ret *= f(cv[i], 1);
		else
			ret *= f(cv[i], 0);
		ret %= pmod;
	}
	return ret;
}
long long gao4(vector<int> a, vector<int> b, vector<int> c) {
	int i, j;
	int n = a.size();
	int lab = a[0] + b[0];
	for (i = 0; i < n; ++i) {
		if (lab > a[i] + b[i])
			lab = a[i] + b[i];
	}
	int lac = a[0] + c[0];
	for (i = 0; i < n; ++i) {
		if (lac > a[i] + c[i])
			lac = a[i] + c[i];
	}
	int lbc = b[0] + c[0];
	for (i = 0; i < n; ++i) {
		if (lbc > b[i] + c[i])
			lbc = b[i] + c[i];
	}
	if (lab + lbc != lac)
		return 0;
	int an = 1, cn = 1;
	map<int, int> am, cm;
	vector<int> av[55], bv, cv[55];
	for (i = 0; i < n; ++i) {
		if (c[i] - b[i] == lbc && a[i] - b[i] == lab)
			bv.push_back(i);
		else if (c[i] - b[i] == lbc) {
			if (abs(a[i] - b[i]) > lab)
				return 0;
			int pos = a[i] - b[i];
			int p = am[pos];
			if (p == 0) {
				p = an++;
				am[pos] = p;
			}
			av[p].push_back(i);
		} else if (a[i] - b[i] == lab) {
			if (abs(c[i] - b[i]) > lbc)
				return 0;
			int pos = c[i] - b[i];
			int p = cm[pos];
			if (p == 0) {
				p = cn++;
				cm[pos] = p;
			}
			cv[p].push_back(i);
		}
	}
	long long ret = f(bv, 0);
	for (i = 1; i < an; ++i) {
		int rr = 0;
		for (j = 0; j < av[i].size(); ++j) {
			if (a[av[i][j]] + b[av[i][j]] == lab)
				rr++;
		}
		if (abs(a[av[i][0]] - b[av[i][0]]) != lab && rr != 1)
			return 0;
		if (abs(a[av[i][0]] - b[av[i][0]]) != lab)
			ret *= f(av[i], 1);
		else
			ret *= f(av[i], 0);
		ret %= pmod;
	}
	for (i = 1; i < cn; ++i) {
		int rr = 0;
		for (j = 0; j < cv[i].size(); ++j) {
			if (c[cv[i][j]] + b[cv[i][j]] == lbc)
				rr++;
		}
		if (abs(c[cv[i][0]] - b[cv[i][0]]) != lbc && rr != 1)
			return 0;
		if (abs(c[cv[i][0]] - b[cv[i][0]]) != lbc)
			ret *= f(cv[i], 1);
		else
			ret *= f(cv[i], 0);
		ret %= pmod;
	}
	return ret;
}
int UnknownTree::getCount(vector<int> a, vector<int> b, vector<int> c) {
	long long ans = 0;
	int i;
	for (i = 0; i < a.size(); ++i)
		cc[i] = a[i];
	ans += gao1(a, b, c);
	ans += gao1(a, c, b);
	ans += gao1(b, a, c);
	ans += gao2(a, b, c);
	ans += gao2(b, a, c);
	ans += gao2(a, c, b);
	ans += gao2(c, a, b);
	ans += gao2(b, c, a);
	ans += gao2(c, b, a);
	ans += gao3(b, c, a);
	ans += gao4(a, b, c);
	ans += gao4(a, c, b);
	ans += gao4(b, a, c);
	return ans % pmod;
}


Merlininice
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
SRM系统框架
11-28
多年SRM实施经验总结，对希望从事SRM实施或规划的同学们有帮助
vmware srm 安装
12-28
可提供先进的复制、灾难恢复管理、无中断测试、自 动故障切换和故障恢复以及计划内迁移功能。它可以管理从生产数据中心到灾难恢复站点的故障切 换，以及两个具有活动工作负载的站点之间的故障切换，甚至还可以将多个...
SRM 565 div1
make_it_for_good的博客
11-10 274
250pts直接dp，f[i][j]表示到i花费为j怪的最大价值和。#include <bits/stdc++.h> using namespace std; #define ll long long int n; ll f[51][110],ans; void upd(ll &x,ll y){x=max(x,y);} class MonstersValley { public: int m
【SRM 565 UnknownTree】计数 分类讨论
Owen_hzt的专栏
01-12 737
一个有N + 3个点的树，告诉你123号点到其他点的距离，求合法的边权为正整数的树个数。 #include #include #include #include #include #include #define Rep(i, x, y) for (int i = x; i <= y; i ++) #define Dwn(i, x, y) for (int i = x; i >= y
topcoder srm 565 div1
weixin_34361881的博客
01-08 118
problem1 link $f[i][j]$表示经过前$i$个怪物之后，花费$j$个硬币可以得到的最大值。 problem2 link 设$nim[i]$表示数字$i$的nim值。那么题目就是求有多少个区间$[L,R]$满足这个区间内所有数字的nim值的抑或不是０．通过记录前缀和的个数，可以$O(n)$计算得到。现在问题是计算每个值的nim值。 可以用数学归纳法证明$i$的nim值等于$...
Topcoder/SRM565
wangyongcs的专栏
12-31 393
SRM565 DIV II: Problem Statement      A divisible sequence with starting value N and length H is a sequence of positive integers with the following properties: The sequence has H elements, labeled
srm 565 div 2
JWZH_CS的专栏
12-23 515
250 水题 500 枚举每个怪打还是不打，dfs水过 1000  题意：定义一个序列S[0] = N，length(S) = H,且 A[i+1] | A[i]，给定N,H,问有多少个数列满足要求。 将N质因子分解后，发现A[i+1]的质因子集是A[i]的子集，，然后就可以想象有n个盒子，每个盒子里有ai个球，每次选任意个（可以是0个），求有多少种选法。每个盒子都可以看成是独立的。然
Topcoder SRM 565 Div.2
Canon of Code
12-22 739
Pt.250: 水题，略过 #include #include #include #include #include using namespace std; class ValueHistogram{ public: vector build(vector fuck) { int n,i,j,max = 0,sum[10] = {0}; string xx; v
SRM565 DIV2
MADNESS
01-03 584
http://community.topcoder.com/stat?c=round_overview&er=5&rd=15187 250-point problem Problem Statement   A histogram of a given collection of numbers graphically represents the frequ
srm 565 div2
weixin_30522095的博客
12-28 66
1、水题，简单的画图，略过。 2、用到了贪心算法，需要把当前n个金币可以买到的最强能力保存下来，然后在可以达到终点的路径中取最小的即可。 1 #include <iostream> 2 #include <string> 3 #include <vector> 4 #include <cstdlib> 5 #include ...
srm
04-01
标题中的“srm”可能指的是“Statistical Regularization Method”（统计正则化方法）或者“Software Requirements Management”（软件需求管理）。在这个上下文中，由于没有更具体的信息，我们可以假设这是关于使用...
SRM SRM 平台功能介绍.pdf
01-04
SAP SRM 介绍
VSAN和SRM.rar
07-05
【标题】"VSAN与SRM"涉及到的是VMware虚拟化环境中的两个关键组件：Virtual SAN（VSAN）和Site Recovery Manager（SRM）。这两个工具在企业级数据中心中发挥着至关重要的作用，确保业务连续性和灾难恢复能力。 VSAN...
<-这是目录-> SRM div1 hard solutions(52)
Merlininice的专栏
12-19 2316
刚从百度空间搬过来。。百度空间没有置顶功能差评啊。。。 从599向后，如果有做那场SRM，会写前两题的简单题解。 （前提有做出来第三题233333，因为此帖本意是第三题练习记录做过哪些，tc你们懂得。。。） practice room练习的话不准备做前两题。 599 598 597 596 595 594 593 592 591 590 5
SRM 601
Merlininice的专栏
12-23 1652
唔。。。做悲剧了- -  250 一开始以为苹果的包和桔子的包不是同一个包。。。500 根本没思路。。950 看的晚了（早看了也写不完。。） 250pt 枚举从每包里那几个，然后找到能拿的最小的苹果数和最大的苹果树（中间那些情况都可以达到） ans += maxapple-minapple+1; 500pt 枚举哪一位开始不同（比如 r），然后dp dp[2005][3
SRM 613
Merlininice的专栏
03-25 1627
因为博客被奇怪的封掉了两天，于是。。。现在才写 easy： 排序，然后枚举分界位置，分界左边向右走，分界右边向左走 medium： 因为The difference high - low will be less than or equal to 100,000. 所以如果选了两个不同的，gcd肯定小于100000 枚举1-100000，算出gcd是
SRM 614
Merlininice的专栏
05-05 1360
这个矩阵很好构造出来，然后算期望要高斯消元，直接算肯定超时，有神奇的方法
SRM 619
Merlininice的专栏
05-06 1254
easy：  如果每堆石头不全为1，那么每次我们总能取一堆石头分给另外两堆，堆数-1，并且新的局面肯定有一堆的个数大于1。 于是，如果每堆石头数都为1 -> lose，否则的话判断堆数奇偶即可（堆数==1要特判） medium： 树形dp，dp[i][j]，表示以i为根的子树分配完成，并且i节点会j技能的最小花费。 转移：枚举i点在自己的department所用的技能(k)，然后
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值