ZROI 提高十连测 DAY3

由于我不太会写 觉得从比赛开始就冷静分析。然后看完三道题心态有点爆炸没有紧扣题目的性质。

这个心态是不可取的尽量不要有畏难心理 不要草草的写暴力。

LINK：[最长01子序列](http://zhengruioi.com/contest/399/problem/960)

对于一个序列要求最长01子序列 显然不太能写感觉无从下手的样子 不妨简化一下问题 先推出来一些性质。

如果求连续的最长01子序列且满足题目中的性质我们显然是根据每一个1进行统计答案。我们直接扫一遍即可。

如果是序列呢 我们发现这个序列两个1之间相差的x是关键 因为没有这个x我们不知道自己算的是哪一种。

这个比较暴力的做法是暴力枚举x 然后统计答案 如何统计答案？显然的有贪心前面的0能选就选了这样显然对后面会更优的。

不那么暴力的做法 事实上我有一个假做法 当时以为这是一个单峰函数 但是结果并非如此因为我无法证明这是单峰函数。

实际上也并非单峰函数 因为这个某个x与其相关的y的值和1和0的个数以及其分布位置是有关的 所以这并非单峰函数。(脑残写了一个三分

害怕超时所以把1抽了出来然后预处理了两个数组 暴力枚举x 然后贪心的判断 然后A掉了这道题。

说起来 这个也算是极大的优化吧。考虑正解：正解和上述做法有异曲同工之妙但是复杂度是有保证的。

都是暴力枚举 然后 匹配的话为了为了每次匹配到下一个1决定二分寻找下一个1 显然对于我们枚举的x 1的个数为n/x

那么我们要二分1个个数次 所以 1个数为n/1+n/2+..n/n显然根据调和级数我们1的个数最多有lnn个所以总复杂度为nlnnlogn

还算很好写..
```
//#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define INF 2147483646
#define ll long long
#define R register
using namespace std;
char buf[1'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch>1;
if(p[mid]-p[w]>=x)r=mid;
else l=mid+1;
}
return r;
}
inline int solve(int x)
{
int cnt=0,st=0;mark=0;
while(1)
{
int s=calc(st,x);
if(s==top+1)
{
if(w[st]>=x)
{
cnt+=x;
break;
}
else
{
--mark;
--cnt;
break;
}
}
st=s;
++mark;
cnt+=x+1;
}
if(!mark)return 0;
return cnt;
}
int main()
{
freopen("1.in","r",stdin);
scanf("%s",a+1);n=strlen(a+1);
last=0;
for(int i=1;i=1;--i)
if(a[i]=='1')
{
w[mark]=last-i-1+w[mark+1];
last=i;
--mark;
}
for(int i=0;iT2 绝对是一个比较神仙的题目 我是很难想到的解法 把到每一个点的 路径长度压缩到log1.1的级别 这个还是要推出一些性质 但是性质推出来之后那么一切也就很好写了。

对于每次询问我们直接二分即可 对于合并 采用归并可以做到O(n)而sort就显得比较满了复杂度也显得比较高因为复杂度是近乎mlognlogn这个的复杂度甚至会更高一点 T了也...我不清楚别问我

关于性质这个性质我是没有网这个方面去思考 可能思考也思考不出来什么 看题解也看了半天可能是数学没学好 对对数函数一点也不敏感。

性质是这样的对于三个边权 x y z来说 如果满足 $\frac{1}{1.1z}为什么?为什么？对于一个询问dis考虑三种情况：

$dis$\frac{1}{1.1z}$dis>y$ 那么y显然是不符合条件的。综上y真没用。

所以说log1.1(maxdis)是距离总数。当然写完了还是很虚 觉得这个东西不是我自己推出来的非常不真切 也对log1.1的值没有大概的把握所以比较自闭。

然而我 思考了一晚上还是觉得不真切 因为这个式子真的很难列出来 但是这个方向的优化还是可以写的。
```
const ll MAXN=200010;
ll n,m,Q,t,h,len;
ll ru[MAXN],q[MAXN],tmp[MAXN];
ll f[MAXN][500],flag[MAXN];
ll lin[MAXN],ver[MAXN],nex[MAXN],e[MAXN];
inline void add(ll x,ll y,ll z)
{
ver[++len]=y;
nex[len]=lin[x];
lin[x]=len;
e[len]=z;
}
inline ll get_x(ll x){return x+x/10ll;}
inline void modify(ll x,ll y,ll z)
{
if(!flag[x])return;
ll i=1,j=1,cnt=0;
for(ll k=1;kflag[y])
{
tmp[++cnt]=f[x][i]+z;
++i;
}
else tmp[++cnt]=f[y][j],++j;
}
f[y][flag[y]=1]=tmp[1];f[y][0]=-200;
for(ll i=2;i>1;
if(f[x][mid]>=y)r=mid;
else l=mid+1;
}
ll w=get_x(y);
if(r!=flag[x]+1&&f[x][r]