2021-09-30

聪明的质检员

额，这道题刚开始给我整蒙了，看这数据就……会T掉。所以：
前缀和+二分：
在二分W参数的函数中计算能算入的矿石的前缀和，在w[i]>=W时这个i矿石会在统计里（若<W就不管它了直接pre[i]=pre[i-1]
矿石价值和是：pre_v[i]=pre_v[i-1]+v[i]，前面的和加上当前这一个i矿石;
矿石数量和是：pre_n[i]=pre_n[i-1]+1。
就可以愉快的写出

for(int i=1;i<=n;i++){
		if(f[i].val>=x)pre_n[i]=pre_n[i-1]+1,pre_v[i]=pre_v[i-1]+f[i].w;
		else pre_n[i]=pre_n[i-1],pre_v[i]=pre_v[i-1];
	}

二分也没有什么好说的，就不说了。

但是

很重要的一点

sum的大小是LL_max ！！！

这个点恶心死了很重要！！！

#include<bits/stdc++.h>
#define N 200010
#define LL long long
using namespace std;
LL l[N],r[N],sum=2147483647483647,num,n,pre_n[N],pre_v[N],m,s;
LL minn=2147483647,maxx=-1;
struct node{
	LL w,val;
}f[N];
bool chk(LL x){
	LL ans=0;
	memset(pre_v,0,sizeof(pre_v));
	memset(pre_n,0,sizeof(pre_n));
	for(int i=1;i<=n;i++){
		if(f[i].val>=x)pre_n[i]=pre_n[i-1]+1,pre_v[i]=pre_v[i-1]+f[i].w;
		else pre_n[i]=pre_n[i-1],pre_v[i]=pre_v[i-1];
	}
	for(int i=1;i<=m;i++){
		ans+=(pre_n[r[i]]-pre_n[l[i]-1])*(pre_v[r[i]]-pre_v[l[i]-1]);
	}
	num=llabs(ans-s);
	if(ans>s)return true;
	else return false;
}
int main(){
	scanf("%lld%lld%lld",&n,&m,&s);
	for(int i=1;i<=n;i++){
		scanf("%lld%lld",&f[i].val,&f[i].w);
	}
	minn=0;
	maxx=1e6;
	for(int i=1;i<=m;i++)scanf("%lld%lld",&l[i],&r[i]);
	while(minn<=maxx){
		int mid=(minn+maxx)>>1;
		if(chk(mid))minn=mid+1;
		else maxx=mid-1;
		sum=min(sum,num);
	}
	printf("%lld",sum);
}

完事儿

---

城市交通路网

标准dp题：
以原点作为参考，向后dp，每一层作为状态，然后……
你就发现找不来路了……
所以应当：反过来想，逆序dp，然后就完事

#include<bits/stdc++.h>
using namespace std;
int a[1001][1001],f[1001],n,x,ma[1001],tot;
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)f[i]=1e9;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			scanf("%d",&x);
			a[i][j]=x;
		}
	}
	f[n]=0;
	for(int j=n;j>=1;j--)
		for(int i=1;i<=n;i++)
			if(a[i][j]){
				if(f[i]>f[j]+a[i][j]){
					f[i]=f[j]+a[i][j];
					ma[i]=j;
				}
			}
	printf("minlong=%d\n",f[1]);
	for(int i=1;i;i=ma[i])printf("%d ",i);
}

---

猪猪储蓄罐（piggy）

本来应该是个完全背包，但是又有点不一样，就是反过来的完全丢包
直接上代码

#include<bits/stdc++.h>
using namespace std;
int s[501][10001],e,w[501],v[501],n,m;
int main(){
	scanf("%d%d",&e,&m);
	memset(s,0x3f,sizeof(s));
	m-=e;
	scanf("%d",&n);
	s[0][0]=0;
	for(int i=1;i<=n;i++){
		scanf("%d%d",&v[i],&w[i]);
	}
	for(int i=1;i<=n;i++){
		for(int j=0;j<=m;j++){
			for(int k=0;k*w[i]<=j;k++){
				s[i][j]=min(s[i][j],s[i-1][j-k*w[i]]+k*v[i]);
			}
		}
	}
	int ans=s[n][m];
	if(ans==s[n+1][m+1])printf("This is impossible.");
	else printf("The minimum amount of money in the piggy-bank is %d.",ans);
}
/*
10 110
2
1 1
30 50
*/

---

金明的预算方案(budget)

分组背包？？？
很像，但不是，又可以看出附件很少，所以考虑将其分为”单主，主+附1，主+附2，主+附1+附2，所以就是分组背包了

#include<bits/stdc++.h>
using namespace std;
const double eps(1e-8);
const int maxn=1100;
const long long maxint=-1u>>1;
const long long maxlong=maxint*maxint;
typedef long long lint;
int f[41000],n,m,v[maxn],p[maxn],q[maxn],ans;
void work()
{
    for (int i=1; i<=m; i++){
        bool flag1=false,flag2=false;
        int t1=0,t2=0;
        if (q[i]) continue;
        for (int j=i+1; j<=m; j++){
        	if (q[j]==i){
                flag1=true;
                t1=j;
                break;
            }
		}
        for (int j=t1+1; j<=m; j++){
        	if (q[j]==i){
                flag2=true;
                t2=j;
                break;
            }
		}
        for (int j=n; j>=v[i]; j--){
            f[j]=max(f[j],f[j-v[i]]+v[i]*p[i]);
            int x=j-v[i]-v[t1];
            if (x>=0 && flag1)f[j]=max(f[j],f[x]+v[i]*p[i]+v[t1]*p[t1]);
            x=j-v[i]-v[t2];
            if (x>=0 && flag2)f[j]=max(f[j],f[x]+v[i]*p[i]+v[t2]*p[t2]);
            x=j-v[i]-v[t1]-v[t2];
            if(x>=0 && flag1 && flag2)f[j]=max(f[j],f[x]+v[i]*p[i]+v[t1]*p[t1]+v[t2]*p[t2]);
        }
    }
    printf("%d",f[n]);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)scanf("%d%d%d",&v[i],&p[i],&q[i]);
    work();
    return 0;
}

---

热浪(heatwv)

SPFA模板，抬走

#include<cstdio>
#include<bits/stdc++.h>
using namespace std;
struct node{
	int x,y,z,next;
} a[1240100];
int hd[1240100],s=0,q[1240100],minn[1240100];
bool f[1240100];
void add(int xx,int yy,int zz)
{
	a[++s].x=xx;
	a[s].y=yy;
	a[s].z=zz;
	a[s].next=hd[xx];
	hd[xx]=s;
}
void spfa(int s)
{
	f[s]=1,q[1]=s,minn[s]=0;
	int h=0,t=1;
	do
	{
		h++;
		int tx=q[h];
		for(int i=hd[tx];i;i=a[i].next){
			if(minn[a[i].y]>minn[tx]+a[i].z){
				minn[a[i].y]=minn[tx]+a[i].z;
				if(f[a[i].y]==0){
					t++;
					f[a[i].y]=1;
					q[t]=a[i].y;
				}
			}
		}
		f[tx]=0;
	}while(h<=t);
}
int main()
{
	memset(f,false,sizeof(f));
	memset(minn,127/3,sizeof(minn));  //求最短，初始定最长
	int t,c,ts,te;
	cin>>t>>c>>ts>>te;
	for(int i=1;i<=c;i++)
	{
		int u,v,w;
		cin>>u>>v>>w;
		add(u,v,w);  //邻接表建立无向边
		add(v,u,w);
	} 
	spfa(ts);
	cout<<minn[te]<<endl;
	return 0;
}

---

Balloons

单调队列栈优化dp？
开始真的没想到，然后有一个结论：$d^2=(r1+r2{)}^2-(r1-r2{)}^2=4\times$r1$\times$r2，可以用单调栈实现以下操作：取出栈顶 j，如果 r_k<r_j不受之前的影响,加入栈顶，为了维护单调栈是递增的，如果r_k>r_j，弹出栈顶 j 保证递增，返回继续循环执行三个操作

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<stack>
#include<map>
#include<algorithm>
#include<map>
#include<cmath>
using namespace std;
const int N=2e5+9;
double r[N];// 
double x[N];//第i个气球与地面在x_i的位置相遇 
int stk[N];
int top,n;
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%lf%lf",&x[i],&r[i]);
    for(int i=1;i<=n;i++)
    {
        while(top>0)
        {
            int di=stk[top];
            double dis=(x[i]-x[di])*(x[i]-x[di])/(4.0*r[di]);
            r[i]=min(r[i],dis);
            if(r[i]<=r[di]) break;//必须保证单调递增 
            top--;
        }
        stk[++top]=i;
        printf("%.4lf\n",r[i]); 
    }
    return 0;
}