3.4

T1
美丽的序列
做法很多 单调栈 线段树…
T2
货币套利
较为复杂
比起一般的图论问题，这道题的边权是相乘而不是相加，于是我们先将所有的边权取对数，利用对数的性质，就变成了一个等价的、边权只用相加的问题。于是原问题变成了一个在图中寻找一个正环的问题。对于这个问题，使用基于dfs的spfa去寻找环。首先把所有点的距离都赋为0，如果某次访问到的顶点已经在栈中，那么就找到了一个正环。（虽然复杂度最坏是O(nm)的，但这种算法实际上还是很快的，尤其是重赋了距离之后。具体可以参见 it
当然本题还有坑点。首先，最后拿在手上的必须是人民币，因此找到正环后还要判断它能否由1号点到达、以及它能否回到1号点。另外，如果不取对数而直接相乘的话，当环的长度很长时，有可能会爆double。
70p

#include<bits/stdc++.h>
using namespace std;
#define eps 1e-10
const int N=400010;
int to[N],ne[N],h[N],v1[N],v2[N],qq[N];
double w[N],d[N],sta[N],cur[N];
int t,num,n,m; 
bool fd,iq[N];
inline int rd(){
	int x=0,k=1; char ch=getchar();
	while(ch<'0'||ch>'9'){ if(ch=='-') k=-1; ch=getchar(); }
	while(ch>='0'&&ch<='9'){ x=(x<<1)+(x<<3)+ch-'0'; ch=getchar(); }
	return x*k;
}
inline double rdd(){
	char c=getchar(); int x=0,y=1,k=1;	bool v=false;
	while(c<'0' || c>'9') { if(c=='-') k=-1; c=getchar(); }
	while((c>='0' && c<='9') || c=='.'){ 
		if(c=='.') v=true;
		else{
			if(v) y=(y<<1)+(y<<3);  
			x=(x<<1)+(x<<3)+c-'0'; 
		}
		c=getchar();
	}
	return x/(double)y;
}
inline void add(int x,int y,double z){
	to[++num]=y;
	ne[num]=h[x];
	h[x]=num;
	w[num]=z;
}
void czh(int u){
	int tp=1; sta[1]=u;
	while(tp && !fd)
	{
		int now=sta[tp];
		double ndist=d[now];
		bool dfsed=false;
      	if(!cur[now]) cur[now]=h[now];
		for(int i=h[now]; i; i=ne[i]) 
               if((i & 1) && v1[to[i]]&&v2[to[i]])
			if(ndist+w[i]-d[to[i]]<-eps)
			{
			    d[to[i]]=ndist+w[i];
				if(!iq[to[i]])
				{
					sta[++tp]=to[i];
					iq[to[i]]=dfsed=true;
				}
				else
					fd=true;
          		break;
			}
		if(!dfsed) 
		{
			tp--;//st.pop();
			iq[now]=false;
		}
	}
}

int main()
{
	//freopen("exchange.in","r",stdin);
	//freopen("exchange.out","w",stdout);
	t=rd(); 
	while(t--)
	{
		num=0; fd=false;
		memset(h,0,sizeof(h));
		memset(w,0,sizeof(w));
		memset(ne,0,sizeof(ne));
		memset(to,0,sizeof(to));
		memset(d,0,sizeof(d));
		memset(cur,0,sizeof(cur));
		memset(v1,0,sizeof(v1));
		memset(v2,0,sizeof(v2));

		n=rd(); m=rd();
		for(int i=1;i<=m;++i)
		{
			int x,y; double v;
			x=rd(); y=rd(); v=rdd();
			v=log(v);
			add(x,y,-v); add(y,x,v);
		}
		
		v1[1]=1; qq[1]=1; int ha=1,ta=1;
		while(ha<=ta)
		{
			int u=qq[ha++];
			for(int i=h[u];i;i=ne[i])
			    if(i&1)
					if(!v1[to[i]]){
				    v1[to[i]]=1;
					qq[++ta]=to[i];
				}
		}
		v2[1]=1; qq[1]=1;  ha=1,ta=1;
		while(ha<=ta)
		{
			int u=qq[ha++];
			for(int i=h[u];i;i=ne[i])
			    if(!(i&1))
					if(!v2[to[i]]){
				    v2[to[i]]=1;
					qq[++ta]=to[i];
				}
		}
		for(int i=1;i<=n;++i)
		if(v1[i]&&v2[i]) czh(i);
	
		if(fd) printf("YES\n");
		else printf("NO\n");
	}
	return 0;
}

std

#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
#define maxn 800010
#define eps 1e-10
int cur[maxn],num[maxn],n,m,nd[maxn],to[maxn],nxt[maxn],teg=0,sta[maxn];
double dist[maxn],w[maxn];
bool iq[maxn],v1[maxn],v2[maxn],ab[maxn],fd=false;
int ni()
{
	char c=getchar();
	while(c<'0' || c>'9') c=getchar();
	int d=0;
	while(c>='0' && c<='9')
	{
		d=d*10+c-48;
		c=getchar();
	}
	return d;
} 
double nr()
{
	char c=getchar();
	while(c<'0' || c>'9') c=getchar();
	int d=0,xi=1;
	bool v=false;
	while((c>='0' && c<='9') || c=='.')
	{
		if(c=='.') v=true;
		else
		{
			if(v) xi=xi*10;
			d=d*10+c-48;
		}
		c=getchar();
	}
	return d/(double)xi;
}
void dfs(int s)
{
	int tp=1; sta[1]=s;
	while(tp && !fd)
	{
		int now=sta[tp];
		double ndist=dist[now];
		bool dfsed=false;
      	if(!cur[now]) cur[now]=nd[now];//没搞懂
		for(int &i=cur[now]; i; i=nxt[i]) //“&”必须写，不然第七个点超时，然而我也不知道为什么
               if((i & 1) && ab[to[i]])
			if(ndist+w[i]-dist[to[i]]<-eps)
			{
				dist[to[i]]=ndist+w[i];
				if(!iq[to[i]])
				{
					sta[++tp]=to[i];//st.push(to[i]);
					iq[to[i]]=dfsed=true;
				}
				else
					fd=true;
          		break;
			}
		if(!dfsed) 
		{
			tp--;//st.pop();
			iq[now]=false;
		}
	}
}
void adeg(int p,int q,double r)
{
	to[++teg]=q;
	nxt[teg]=nd[p];
	nd[p]=teg;
	w[teg]=r;
}
int main()
{
	freopen("exchange.in","r",stdin);
	freopen("exchange.out","w",stdout);
	int t; scanf("%d",&t);
	while(t--)
	{
		memset(cur,0,sizeof(cur));
		teg=0; fd=false;
		memset(num,0,sizeof(num)); memset(nd,0,sizeof(nd));
		memset(to,0,sizeof(to)); memset(nxt,0,sizeof(nxt));
		memset(iq,0,sizeof(iq));
		memset(v1,0,sizeof(v1));
		memset(v2,0,sizeof(v2));
		memset(dist,0,sizeof(dist));
		scanf("%d%d",&n,&m);
		for(int i=1; i<=m; i++)
		{
			int p,q; double r;
			p=ni(),q=ni(),r=nr();
			r=-log(r);
			adeg(p,q,r);
			adeg(q,p,-r);//反向建边
		}
		v1[1]=true; int ta=1,fr=1; sta[1]=1;
		while(ta<=fr)
		{
			int np=sta[ta++];
			for(int i=nd[np]; i; i=nxt[i]) if(i & 1)正向边 1节点能到的点。
				if(!v1[to[i]])
				{
					v1[to[i]]=true;
					sta[++fr]=to[i];
				}
		}
		v2[1]=true; ta=1,fr=1; sta[1]=1;
		while(ta<=fr)
		{
			int np=sta[ta++];
			for(int i=nd[np]; i; i=nxt[i]) if(!(i & 1))//反向边 能回到1节点。
				if(!v2[to[i]])
				{
					v2[to[i]]=true;
					sta[++fr]=to[i];
				}
		}
		for(int i=1; i<=n; i++)
			if(v1[i] && v2[i]) ab[i]=true; else ab[i]=false;
		for(int i=1; i<=n; i++) if(ab[i]) dfs(i);
		if(fd) printf("YES\n");
		else printf("NO\n");
	}
	fclose(stdin);
	fclose(stdout);
	return 0;
}

T3
大海战
更复杂
这道题目的原型是某次CF DIV2 的D题，当时naïve的我去写了线段树。后来联想到《借教室》才发现我真的naïve了。于是魔改了一下（也就是套了个DP）就变成了这个样子。接下来我们分测试点分析。
第一个点没有询问，只需要判断初始棋盘能否放下就行了。

——直接
把所有ci* ti加起来,判是否小于等于n即可。但是大家似乎忘记了ci和ti的范围——乘起来会爆int，造成许多P党直接RE……
第2到第8个点，只有一次询问。这个询问把棋盘分成了两段，我们考虑其中较短的那部分，就变成了一个背包问题——要在这部分中放置尽量多的物品，空间浪费得越少越好。（至于剩下的那部分，直接用物品总体积减去前面能够放入的物品体积，与剩余的空间进行比较即可。）
每个战舰的价值即为它的长度，直接多重背包即可。第2、3个点直接是01背包，而后面的点由于c*n较大，朴素的背包会T，需要用单调队列优化得到O(cn)的算法,

可以参见 it
（我尝试卡过二进制拆分的多重背包，但写得优美的话也可以过）。
第9到14个点，只有一种物品。我们二分答案，那么询问序列把棋盘分成了多个部分。假设某一段的长度为L，那么就可以放下L/size个物品。于是我们扫一遍就可以统计最多放置的物品数量，然后再根据要求二分下去即可。
第15到20个点，有两种物品。在二分的基础上（设分为了mid段），我们用f[i][j]表示考虑到前i段，恰好放置j种1号物品的前提下，可以放置的最多2号物品的数量。设第i段的长度为Li，转移方程为f[i][j]=max(f[i-1][j-k]+(Li-k * size1)/size2)，其中0<=k<=j且k*size1<=Li。那么f[mid][num1]就是当前状态可以放置的最多2号物品的数量，与要求进行比较，再二分即可。看上去这个方程是O(n^3)的，
但实际上我们有k<=Li，因此对于确定的j，
考虑所有的i，总共的决策数量不会超过sigma(Li)，即不会超过O(n)个。所以DP的复杂度为O(n^2)，
整个算法的时间复杂度为O(n^2logq),
但远远达不到这个界，稍微优化一下就可以过了。