USTCOJ 1266 Sequence Again

题目链接：http://acm.ustc.edu.cn/ustcoj/problem.php?id=1266

题意理解起来比较简单给一串数字，在其中找最长的一段连续串使得它是1--x的一个排列，比如5 1 2 3 4，这个串本身就是1--5的一个排列，所以输出5，3 1 2 3 3 4

则最长为3 1 2，输出3。至今没有想到比N^2快的算法，所以贴出来给大家探讨一下，如果有好方法，请赐教。

我的N^2算法是这样的，首先找到最小的没有在序列中出现的数字，这是答案的上届，然后我们从上界开始不断减小直到找到第一个满足要求的答案。对于判断某一个x是不是可以使解，我是这样做的，先求一段长为x的序列的和，判断一下是否为（x + 1）*x / 2,如果是继续判断是否出现了1--x，否则整段序列向某个方向移动一位，这样这段和可以O（N）求出，每次到达某一个边界（序列某一端）可以把序列另一端减去首位，得到的和作为下一个判断的初始和。这样也能省一些求和时间，（怕TLE）。这个作法最坏时最小的未出现的数字是N，而答案是1，比如这个序列：n-1， n-1， n-2，n-2， 。。。4， 4， 2， 2，3， 3， 1， 3，3，2，2，4，4。。。n - 2， n - 2， n - 1， n - 1。所以退化了之后是个N^2,算法。水水就过了。。。求更好算法。（Ps，按我的理解这题是不能二分答案来查找的，如果有想法，请留言）

代码如下：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<algorithm>
#include<stack>
#include<deque>
#include<list>
#include<set>
#include<vector>
#include<iomanip>
#include<cctype>
#include<string>
#include<memory>
#include<map>
#include<sstream>
#pragma warning (disable : 4996)
#define mem(a) memset(a, 0, sizeof(a))
#define For(i, n) for(int i = 0; i < (n); ++i)
#define sl(a) strlen(a)
typedef long long LL;
typedef double dou;
const int Mod = 1000000007;
const double eps = 1e-8;
const LL inf = 1e18;
const int inf1 = 1e9;
const int N = 20005;
using namespace std;
#define s(x) (((x) + 1) * (x) >> 1)
int a[N], vis[N];

int main()
{
		int n, i, j, k, t, s, e, sum, flag, up;
		cin >> t;
		while (t--)
		{
				cin >> n;
				sum = flag = 0, mem(vis);
				for (i = 0; i < n; ++i) cin >> a[i], vis[a[i]] = 1;
				for (i = 1; i <= n && vis[i]; ++i);
				up = i;
				for (i = 0; i < up; ++i) sum += a[i];
				if (sum == s(up))
				{
						for (i = 0; i < n; ++i) vis[a[i]] = 1;
						for (i = 1; i <= up && vis[i]; ++i);
						if (i > up) 
						{
								cout << up << endl;
								flag = 1;
								goto sym;
						}
				}
				s = 0, e = up - 1;
				for (i = up - 1; i >= 1; --i)
				{
						if (( i & 1) == ((up - 1) & 1)) 
						{
								sum -= a[e--];
								for (j = 0; j <= n - i; ++j)
								{
										if (sum == s(i))
										{
												for (k = s; k <= e; ++k) vis[a[k]] = 1;
												for (k = 1; k <= i && vis[k]; ++k);
												if (k > i) 
												{
														cout << i << endl;
														flag = 1;
														goto sym;
												}
												else mem(vis);
										}
										else if (e + 1 < n) sum += a[e++ + 1] - a[s++];
								}
						}
						else 
						{
								sum -= a[s++];
								for (j = n - i; j >= 0; --j)
								{
										if (sum == s(i))
										{
												for (k = s; k <= e; ++k) vis[a[k]] = 1;
												for (k = 1; k <= i && vis[k]; ++k);
												if (k > i) 
												{
														cout << i << endl;
														flag = 1;
														goto sym;
												}
												else mem(vis);
										}
										else if (s - 1 >= 0)sum += a[s-- - 1] - a[e--];
								}
						}
				}
				sym : 
				if (flag == 0) cout << 0 << endl;
		}
		return 0;
}