SRM 622

easy：

先floyd求出两两直接最短路，然后n^2枚举每条路，n^2判断他是不是某两点最短路可以路过的路径即可。

medium：

贪心。

从某点开始dfs。

如果孩子返回的分支长度大于maxdist，那把孩子单独分一块。

剩下孩子返回的分支长度排序，最大+次大>maxdist，最大分支单独分一块，这样一直算知道最大+次大<=maxdist，把最大分支长度返回给父亲

hard

数位dp。

Fibonacci base的数，相邻两位不能同时是1.

于是我们dp表示成dp[n][2][2]，后两维表示有没有到达上限和上一位是不是1，上一位是不是1，dp记录是这种方案的个数。

求dp数组和普通数位dp一样。那么求ans的时候。

需要知道到达这种状态的是奇数还是偶数个方案、后面可以填的全部方案是奇数还是偶数，如果都是奇数，那么答案这一位就是1。

代码：

#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 FibonacciXor {
public:
	int find(long long, long long);
};
long long base[105];
int a[105];
void pre() {
	int i;
	base[0] = 1;
	base[1] = 2;
	for (i = 2; i < 105; ++i)
		base[i] = base[i - 1] + base[i - 2];
}
long long pmod = 1000000007;
int dp[105][2][2];
long long b[105];
int ans[105];
void gao(long long x) {
	int n, i, j, k;
	n = 0;
	memset(a, 0, sizeof(a));
	for (i = 75; i >= 0; --i) {
		if (x >= base[i]) {
			a[i] = 1;
			x -= base[i];
		} else
			a[i] = 0;
		b[i] = x;
	}
	memset(dp, 0, sizeof(dp));
	dp[76][0][0] = 1;
	for (i = 75; i >= 0; --i) {
		if (a[i] == 1) {
			//放1
			dp[i][0][1] += dp[i + 1][0][0];
			dp[i][1][1] += dp[i + 1][1][0];
			if (base[i + 2] % 2 == 1 && dp[i + 1][1][0] % 2 == 1)
				ans[i] ^= 1;
			if ((b[i] + 1) % 2 == 1 && dp[i + 1][0][0] % 2 == 1)
				ans[i] ^= 1;

			//放0
			dp[i][1][0] += dp[i + 1][1][1];
			dp[i][1][0] += dp[i + 1][0][1];
			dp[i][1][0] += dp[i + 1][1][0];
			dp[i][1][0] += dp[i + 1][0][0];

		} else {
			//放1
			dp[i][1][1] += dp[i + 1][1][0];
			if (base[i + 2] % 2 == 1 && dp[i + 1][1][0] % 2 == 1)
				ans[i] ^= 1;
			//放0
			dp[i][0][0] += dp[i + 1][0][1];
			dp[i][1][0] += dp[i + 1][1][1];
			dp[i][0][0] += dp[i + 1][0][0];
			dp[i][1][0] += dp[i + 1][1][0];
		}
		dp[i][0][0] %= 2;
		dp[i][0][1] %= 2;
		dp[i][1][0] %= 2;
		dp[i][1][1] %= 2;
	}
}
int FibonacciXor::find(long long A, long long B) {
	pre();
	memset(ans, 0, sizeof(ans));
	gao(A - 1);
	gao(B);
	int ret = 0;
	for (int i = 75; i >= 0; --i)
		ret = (ret * 2 + ans[i]) % pmod;
	return ret;
}