Loi模拟题之Loi水题汇总 MST+DP+DP+树剖

最新推荐文章于 2017-10-24 19:00:46 发布

cgoldsss

最新推荐文章于 2017-10-24 19:00:46 发布

阅读量261

点赞数 1

分类专栏：考试

本文链接：https://blog.csdn.net/qq_36312502/article/details/78312579

版权

考试专栏收录该内容

14 篇文章 0 订阅

订阅专栏

T1：

题目：
http://www.lydsy.com/JudgeOnline/problem.php?id=1821

方法1：
求一遍最小生成树，第n-k+1条边就是最近距离;
找不到代码了……

方法2:
二分+并查集;
注意精度；
由于保留小数点后两位，所以r-l>0.001;

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;

const int MAXN=2000001;
int fst[MAXN],nxt[MAXN],fa[MAXN];
double l,r;
int n,k,tot,cnt;
struct sh {int x,y;}mp[MAXN];
struct hh {int from,to;double cost;}ma[MAXN];

void build(int f,int t,double c)
{
    ma[++tot]=(hh){f,t,c};
    nxt[tot]=fst[f];
    fst[f]=tot;
    return;
}

double get(int i,int j)
{
    return sqrt( (mp[i].x-mp[j].x) * (mp[i].x-mp[j].x) + (mp[i].y-mp[j].y) * (mp[i].y-mp[j].y) );
}

int find(int x)
{
    int r=x,t;
    while(r!=fa[r]) r=fa[r];
    while(x!=r) t=fa[x],fa[x]=r,x=t;
    return r;
}

bool check(double x)
{
    for(int i=1;i<=n;i++) fa[i]=i;
    cnt=0;
    for(int i=1;i<=tot;i++)
    {
        if(ma[i].cost<x)
        {
            int fx=find(ma[i].from),fy=find(ma[i].to);
            if(fx!=fy)  fa[fx]=fy;
        }
    }
    for(int i=1;i<=n;i++) 
        if(fa[i]==i) cnt++;
    if(cnt>=k) return true;
    else return false;
}

void solve()
{
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++) 
        scanf("%d%d",&mp[i].x,&mp[i].y);

    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
            build(i,j,get(i,j));

    for(int i=1;i<=tot;i++) 
        r=max(ma[i].cost,r);

    while(r - l >0.001)
    {
        double mid=(l+r)/2;
        if(check(mid)) l=mid;
        else r=mid;
    }

    printf("%.2lf",l);
    return;
}

int main()
{
    solve();
    return 0;
}

T2:

http://codevs.cn/problem/5251/

s[i]:底数不超过i的方案总数；
f[i]：以i为底数的方案数;

显然：
f[i]=s[i/k];
s[i]=s[i-1]+f[i];

因为模数可能大于n，所以s[0]=1；

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll MAXN=1000001;
ll T,k,p,n;
ll f[MAXN],s[MAXN];

void solve()
{
    scanf("%lld%lld",&T,&p);
    while(T--)
    {
        ll cnt=1;
        s[0]=s[1]=f[1]=1;
        scanf("%lld%lld",&n,&k);

        for(ll i=2;i<=n/k;i++) 
            f[i]=s[i/k],s[i]=(s[i-1]+f[i])%p;

        for(ll i=n;i/=k;cnt++);

        printf("%lld %lld\n",s[n/k]%p,cnt%p);
    }
    return;
}

int main()
{
    solve();
    return 0;
}

T3:

题目：
http://codevs.cn/problem/1684/

dp[i]：卡门走了i的高度最多存活的时间；
加了包装的01背包；
但其实高度就是”体积”，时间就是”价值”；
滚动数组优化；
详见：http://blog.csdn.net/qq_36312502/article/details/77196224

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;

const int MAXN=20001;
int dp[MAXN];
int D,G;
struct hh
{
    int T,F,H;
}ma[MAXN];

bool cmp(hh a,hh b)
{
    return a.T < b.T;
}
void solve()
{
    scanf("%d%d",&D,&G);
    for(int i=1;i<=G;i++) 
        scanf("%d%d%d",&ma[i].T,&ma[i].F,&ma[i].H);
    sort(ma+1,ma+G+1,cmp);
    dp[0]=10;

    for(int i=1;i<=G;i++)
    {
        for(int j=D;j>=0;j--)
        {
            if(ma[i].T<=dp[j])
            {
                if(dp[D]>=ma[i].T)
                {
                    cout<<dp[D]<<'\n';
                    return;
                }
                dp[j+ma[i].H]=max(dp[j],dp[j+ma[i].H]);
                dp[j]+=ma[i].F;
            }   
        }
    }

    cout<<dp[0]<<'\n';
    return;
}

int main()
{
    solve();
    return 0;
}

T4：

详见：http://blog.csdn.net/qq_36312502/article/details/78121945

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;

const int MAXN=600000+6;
int fst[MAXN],nxt[MAXN],cdis[MAXN],dis[MAXN],tmp[MAXN];//建图; 
int son[MAXN],dfn[MAXN],siz[MAXN],fa[MAXN],deep[MAXN],top[MAXN];//树剖 
int n,m,tot,tim,l=-1,r;
//dfn[i]:第i个被访问的点; 
//cdis[i]:fa[i]到i的边权;
//tmp[]:差分数组; 
struct hh {int from,to,cost;}ma[MAXN];
struct sh {int from,to,lca,dis;}mp[MAXN];

void build(int f,int t,int c)
{
    ma[++tot]=(hh){f,t,c};
    nxt[tot]=fst[f];
    fst[f]=tot;
    return;
}

void dfs1(int x,int f,int c)//处理树上信息
{
    fa[x]=f,siz[x]=1;
    deep[x]=deep[f]+1;
    dfn[++tim]=x,dis[x]=c;
    for(int i=fst[x];i;i=nxt[i])
    {
        int v=ma[i].to;
        if(v==f) continue;
        cdis[v]=ma[i].cost;//边权映射为点权; 
        dfs1(v,x,c+ma[i].cost),siz[x]+=siz[v];
        if(!son[x] || siz[v]>siz[son[x]]) son[x]=v;//找重儿子; 
    }
    return;
}

void dfs2(int x,int st)//树剖;
{
    top[x]=st;
    if(!son[x]) return;
    dfs2(son[x],st);
    for(int i=fst[x];i;i=nxt[i])
    {
        int v=ma[i].to;
        if(v==fa[x] || v==son[x]) continue;
        dfs2(v,v);
    }
}

int get_lca(int x,int y)//树剖求x和y的lca 
{
    int fx=top[x],fy=top[y];
    while(fx!=fy)
    {
        if(deep[fx]<deep[fy]) swap(fx,fy),swap(x,y);
        x=fa[fx],fx=top[x];
    }
    return deep[x]>deep[y]?y:x;
}

bool check(int x)//二分; 
{
    int cnt=0,limi=-1;
    memset(tmp,0,sizeof(tmp));
    for(int i=1;i<=m;i++)
    {
        if(mp[i].dis>x)//找出所有长度大于x的路径； 
        {
            limi=max(mp[i].dis-x,limi);//找最大差值; 
            cnt++;
            tmp[mp[i].from]++,tmp[mp[i].to]++;//差分; 
            tmp[mp[i].lca]-=2;
        }
    }

    if(!cnt) return true;
    for(int i=n;i>=1;i--) tmp[fa[dfn[i]]]+=tmp[dfn[i]];

    for(int i=2;i<=n;i++)
        if(tmp[i]==cnt && cdis[i]>=limi) return true;
//存在一条边是所有大于x的路径的交,并且边权大于最大差值, 
//删掉此边,所有路径长度都会小于或等于x; 
    return false;
}

void solve()
{

    int f,t,c;

    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
    {
        scanf("%d%d%d",&f,&t,&c);
        build(f,t,c),build(t,f,c);
        r+=c;
    }
    dfs1(1,0,0),dfs2(1,1);

    for(int i=1;i<=m;i++)
    {//处理路径信息,dis为路径长度; 
        scanf("%d%d",&mp[i].from,&mp[i].to);
        mp[i].lca=get_lca(mp[i].from,mp[i].to);
        mp[i].dis=dis[mp[i].from]+dis[mp[i].to]-2*dis[mp[i].lca];
    }

    while(r - l > 1)//二分; 
    {
        int mid=(l+r)>>1;
        if(check(mid)) r=mid;
        else l=mid;
    }

    printf("%d",r);
    return;

}

int main()
{
    solve();
    return 0;
}