20161006模拟

评测数据下载：https://yunpan.cn/cvVysDxL7F3KC （提取码：b07f）

分析：

T1 KMP+矩阵乘法（参考 BZOJ 1009 GT考试）

T2 斐波那契数列变形（数据很大，要用矩阵乘法），注意判一下

T3 不是匈牙利，而是树形dp，自己慢慢悟吧

T4 蒟蒻只会30分的做法：求图的直径，再/2，就是ans

更正：第三组：不存在相同的字符|str|=26，26<=n<=100

更正：输出的顺序保证a<b

更正：输出样例：0 1000000006

更正：模数1000000007

T1

60分代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=1e5+7;
const int mod=1e9+7;
int n,cnt,len;
char str[N],pat[N];
ll f[N];
void dfs(int now){
    if(now==n){
        if(++cnt>=mod) cnt%=mod;
        return ;
    }
    bool flag=(now<len-1)|(!equal(str+now-len+1,str+now,pat));
    for(char x=‘a‘;x<=‘z‘;x++){
        if(!flag&&x==pat[len-1]) continue;
        str[now]=x;
        dfs(now+1);
    }
}
ll fpow(ll a,ll p){
    ll res=1;
    for(;p;p>>=1,a=a*a%mod) if(p&1) res=res*a%mod;
    return res;
}
int main(){
    freopen("helloworld.in","r",stdin);
    freopen("helloworld.out","w",stdout);
    while(scanf("%d",&n)==1){
        scanf("%s",pat);
        len=strlen(pat);
        cnt=0;
        if(!n){puts("0");continue;}
        if(n<=5){
            dfs(0);
            printf("%d\n",cnt);
        }
        else if(len==1){
            printf("%d\n",fpow(25,n));
        }
        else{
            f[0]=1;
            for(int i=1;i<=len;i++) f[i]=f[i-1]*26%mod;
            for(int i=len;i<=n;i++) f[i]=(f[i-1]*26-f[i-len]+mod)%mod;
            printf("%d\n",(int)f[n]);
        }
    }
    return 0;
}

Std

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 10100
#define MAXM 110
#define MOD 1000000007
int dp[MAXN][MAXM];
int fail[MAXM];
int trs[MAXN][26];
char str[MAXN];
inline void deal(int &x,int y){
    x+=y;
    if(x>=MOD) x-=MOD;
}
int main(){
    freopen("helloworld.in","r",stdin);
    freopen("helloworld.out","w",stdout);
    int n;
    while(~scanf("%d%s",&n,str+1)){
        memset(dp,0,sizeof(dp));
        int m=strlen(str+1);
        fail[1]=0;
        for(int i=2;i<=m;i++){
            int p=fail[i-1];
            while(p&&str[p+1]!=str[i]) p=fail[p];
            if(str[p+1]==str[i]) fail[i]=p+1;
        }
        memset(trs,-1,sizeof(trs));
        memset(trs[0],0,sizeof(trs[0]));
        for(int i=0;i<m;i++) trs[i][str[i+1]-‘a‘]=i+1;
        for(int i=0;i<=m;i++)
            for(int j=0;j<26;j++)
                if(trs[i][j]==-1)
                    trs[i][j]=trs[fail[i]][j];
        dp[0][0]=1;
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
                for(int k=0;k<26;k++)
                    deal(dp[i+1][trs[j][k]],dp[i][j]);
        long long ans=0;
        long long x=1;
        for(int i=n;i>=0;i--)
        {
            ans = (ans+x*dp[i][m])%MOD;
            if(i)x=x*26%MOD;
        }
        printf("%d\n",(int)((x-ans)%MOD+MOD)%MOD);
    }
    return 0;
}

T2

AC代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#ifdef unix
#define LL "%lld"
#else
#define LL "%I64d"
#endif
using namespace std;
#define ll long long
#define N 3
const ll mod=1e9+7;
ll a[N][N],b[N][N],c[N][N];
ll p,q,a1,a2,n;
inline const ll read(){
    register ll x=0,f=1;
    register char ch=getchar();
    while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)f=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}
    return x*f;
}
void work(){
    a[1][1]=0;a[1][2]=q;a[2][1]=1;a[2][2]=p;
    b[1][1]=0;b[1][2]=q;b[2][1]=1;b[2][2]=p;
    while(n){
        if(n&1){
            for(int i=1;i<=2;i++){
                for(int j=1;j<=2;j++){
                    for(int k=1;k<=2;k++){
                        c[i][j]=(c[i][j]+a[i][k]*b[k][j]%mod)%mod;
                    }
                }
            }
            memcpy(a,c,sizeof c);
            memset(c,0,sizeof c);
        }
        for(int i=1;i<=2;i++){
            for(int j=1;j<=2;j++){
                for(int k=1;k<=2;k++){
                    c[i][j]=(c[i][j]+b[i][k]*b[k][j]%mod)%mod;
                }
            }
        }
        memcpy(b,c,sizeof c);
        memset(c,0,sizeof c);
        n>>=1;
    }
}
int main(){
    freopen("gcd.in","r",stdin);
    freopen("gcd.out","w",stdout);
    p=1;q=1;a1=1;a2=2;
    n=read();ll bf=n;
    if(n==1){printf("1 1");return 0;}
    if(n%mod==0){printf("1 ");printf(LL,n-1);return 0;}
    n-=1;
    work();
    ll ans1=(a[1][1]%mod+a[2][1]%mod)%mod;
    n=bf;
    work();
    ll ans2=(a[1][1]%mod+a[2][1]%mod)%mod;
    if(ans1>ans2) swap(ans1,ans2);
    printf(LL,ans1);putchar(‘ ‘);printf(LL,ans2);
    return 0;
}

T3

AC代码1：

#include<cstdio>
#include<cstring>
#include<vector>
#define ll long long
using namespace std;
const int N=1e5+10;
const int mod=1e9+7;
const int inf=0x3f3f3f3f;
struct node{
    int v,next;
}e[N<<1];
int T,opt,n,tot,head[N],q[N],fa[N],mx[N][2],dp[N][2];
inline const int read(){
    register int x=0,f=1;
    register char ch=getchar();
    while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)f=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}
    return x*f;
}
inline void add(int x,int y){
    e[++tot].v=y,e[tot].next=head[x],head[x]=tot;
}
void work(){
    int h=0,t=1;
    q[1]=1;
    while(h<t){
        int now=q[++h];
        for(int i=head[now];i;i=e[i].next){
            int v=e[i].v;
            if(v==fa[now]) continue;
            fa[v]=now;
            q[++t]=v;
        }
    }
    for(int i=t;i;i--){
        int now=q[i];
        dp[now][0]=0;
        int mx0=0,mx1=0,dp0=1,dp1=0;
        for(int i=head[now];i;i=e[i].next){
            int v=e[i].v;
            if(v==fa[now]) continue;
            mx0+=mx[v][1];
            dp0=(ll)dp0*dp[v][1]%mod;
        }
        dp[now][0]=dp0;
        mx[now][0]=mx0;
        vector<int> vec1,vec2;
        int delta=inf;
        for(int i=head[now];i;i=e[i].next){
            int v=e[i].v;
            if(v==fa[now]) continue;
            vec1.push_back(dp[v][1]);
            vec2.push_back(dp[v][1]);
            delta=min(delta,mx[v][1]-mx[v][0]);
        }
        int vsize=vec1.size();
        for(int i=1;i<vsize;i++) vec1[i]=(ll)vec1[i]*vec1[i-1]%mod;
        for(int i=vsize-2;i>=0;i--) vec2[i]=(ll)vec2[i]*vec2[i+1]%mod;
        int cnt=0;
        for(int i=head[now];i;i=e[i].next){
            int v=e[i].v;
            if(v==fa[now]) continue;
            if(mx[v][1]-mx[v][0]==delta)
                dp1=(dp1+(ll)(cnt?vec1[cnt-1]:1)*(cnt==vsize-1?1:vec2[cnt+1])%mod*dp[v][0])%mod;
            cnt++;
        }
        mx[now][1]=mx0-delta+1;
        dp[now][1]=dp1;
        if(mx[now][0]==mx[now][1]) dp[now][1]=(dp[now][0]+dp[now][1])%mod;
        if(mx[now][0]>mx[now][1]) dp[now][1]=dp[now][0],mx[now][1]=mx[now][0];
    }
}
int main(){
    freopen("hungary.in","r",stdin);
    freopen("hungary.out","w",stdout);
    T=read();opt=read();
    while(T--){
        memset(dp,0,sizeof dp);
        memset(mx,0,sizeof mx);
        memset(fa,0,sizeof fa);
        memset(head,0,sizeof head);
        tot=0;
        n=read();
        for(int i=1,x,y;i<n;i++) x=read(),y=read(),add(x,y),add(y,x);
        work();
        opt==1?printf("%d\n",mx[1][1]):printf("%d %d\n",mx[1][1],dp[1][1]);
    }
    return 0;
}

ylf‘sAC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#define maxn 100010
#define mod 1000000007
#define ll long long
using namespace std;
ll T,P,n,head[maxn],num,f[maxn][2],g[maxn][2],L,R[maxn],l,r[maxn],son[maxn];
struct node{
    ll v,pre;
}e[maxn*2];
ll init(){
    ll x=0,f=1;char s=getchar();
    while(s<‘0‘||s>‘9‘){if(s==‘0‘)f=-1;s=getchar();}
    while(s>=‘0‘&&s<=‘9‘){x=x*10+s-‘0‘;s=getchar();}
    return x*f;
}
void Add(ll from,ll to){
    num++;e[num].v=to;
    e[num].pre=head[from];
    head[from]=num;
}
void Clear(){
    num=0;
    memset(f,0,sizeof(f));
    memset(g,0,sizeof(g));
    memset(head,0,sizeof(head));
}
void DP(ll now,ll from){
    g[now][0]=1;
    ll mx,sum;
    for(int i=head[now];i;i=e[i].pre){
        ll v=e[i].v;
        if(v==from)continue;
        DP(v,now);//x不连儿子 儿子们可连可不连
        mx=max(f[v][1],f[v][0]);sum=0;
        if(mx==f[v][1])sum+=g[v][1];
        if(mx==f[v][0])sum+=g[v][0];
        g[now][0]=g[now][0]*sum%mod;
        f[now][0]+=mx;
    }
    //x连某个儿子 这个不选 其他的连或者不连
    L=0;l=1;ll S=0;
    for(int i=head[now];i;i=e[i].pre)
        if(e[i].v!=from)son[++S]=e[i].v;
    R[S+1]=0;r[S+1]=1;
    for(int i=S;i>=1;i--){
        ll v=son[i];sum=0;
        mx=max(f[v][1],f[v][0]);
        if(mx==f[v][1])sum+=g[v][1];
        if(mx==f[v][0])sum+=g[v][0];
        R[i]=R[i+1]+mx;
        r[i]=r[i+1]*sum%mod;
    }
    for(int i=1;i<=S;i++){
        ll v=son[i];
        mx=L+f[v][0]+R[i+1]+1;
        if(mx>f[now][1]){
            f[now][1]=mx;
            g[now][1]=l*g[v][0]%mod*r[i+1]%mod;
        }
        else if(mx==f[now][1])
            g[now][1]=(g[now][1]+l*g[v][0]%mod*r[i+1]%mod)%mod;
        sum=0;
        mx=max(f[v][1],f[v][0]);
        if(mx==f[v][1])sum+=g[v][1];
        if(mx==f[v][0])sum+=g[v][0];
        l=l*sum%mod;L+=mx;
    }
}
int main()
{
    freopen("hungary.in","r",stdin);
    freopen("hungary.out","w",stdout);
    T=init();P=init();
    while(T--){
        n=init();
        ll u,v;Clear();
        for(int i=1;i<n;i++){
            u=init();v=init();
            Add(u,v);Add(v,u);
        }
        DP(1,0);ll sum,mx;
        mx=max(f[1][0],f[1][1]);sum=0;
        if(mx==f[1][0])sum+=g[1][0],sum%=mod;
        if(mx==f[1][1])sum+=g[1][1],sum%=mod;
        if(P==1)cout<<mx<<endl;
        if(P==2)cout<<mx<<" "<<sum<<endl;
    }
    return 0;
}

T4

30分代码：

#include<cstdio>
using namespace std;
const int N=205;
const int inf=0x3f3f3f3f;
inline int max(int a,int b){return a>b?a:b;}
inline const int read(){
    register int x=0,f=1;
    register char ch=getchar();
    while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)f=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}
    return x*f;
}
int n,m;
int f[N][N];
int main(){
    freopen("radius.in","r",stdin);
    freopen("radius.out","w",stdout);
    n=read();m=read();
    for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) f[i][j]=inf;
    for(int i=1,x,y,z;i<=m;i++){
        x=read();y=read();z=read();
        if(f[x][y]>z) f[y][x]=f[x][y]=z;
    }
    for(int k=1;k<=n;k++){
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(i!=j&&j!=k&&i!=k){
                    if(f[i][j]>f[i][k]+f[k][j]){
                        f[i][j]=f[i][k]+f[k][j];
                    }
                }
            }
        }
    }
    int d=0;
    for(int i=1;i<n;i++){
        for(int j=i+1;j<=n;j++){
            if(f[i][j]!=inf) d=max(d,f[i][j]);
        }
    }
    double ans=d/2.0;
    printf("%.2lf",ans);
    return 0;
}

Std

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXV 170000
#define MAXE 170000
#define MAXN 900
#define PROB "radius"
#define INF 0x3f3f3f3f
struct Edge{
    int np,val;
    Edge *next;
    bool flag;
}E[MAXE],*V[MAXV];
int m,n;
int map[MAXN][MAXN];
struct aaa
{
    int x,y,d;
}el[MAXE];
int tope=-1,topl=-1;
void addedge(int x,int y,int z){
    el[++topl].x=x;
    el[topl].y=y;
    el[topl].d=z;
    E[++tope].np=y;
    E[tope].val=z;
    E[tope].next=V[x];
    V[x]=&E[tope];
    E[++tope].np=x;
    E[tope].val=z;
    E[tope].next=V[y];
    V[x]=&E[tope];
}
void init(){
    int i,j,k;
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
            for(k=1;k<=n;k++){
                if(map[j][k]>map[j][i]+map[i][k])
                    map[j][k]=map[j][i]+map[i][k];
            }
        }
    }
}
pair<int,int> seg[MAXN];
int tops=-1,rg;
void set_range(int x){
    tops=-1;
    rg=x;
}
void add_seg(int x,int y){
    tops++;
    if(x>y)throw "E";
    seg[tops].first=x;
    seg[tops].second=y;
}
bool full(){
    int i,x=0;
    sort(seg,&seg[tops+1]);
    for(i=0;i<=tops;i++){
        if(x<seg[i].first)return false;
        if(x>seg[i].second)continue;
        x=seg[i].second+1;
    }
    if(x>rg)return true;
    return false;
}
int main(){
    freopen(PROB".in","r",stdin);
    freopen(PROB".out","w",stdout);
    int i,j,k,x,y,z;
    scanf("%d%d",&n,&m);
    memset(map,INF,sizeof(map));
    for(i=0;i<m;i++){
        scanf("%d%d%d",&x,&y,&z);
        addedge(y,x,z*2);
        map[x][y]=map[y][x]=min(map[x][y],z*2);
    }
    for(i=1;i<=n;i++) map[i][i]=0;
    init();
    int l,r,mid;
    l=0;r=10000006;
    int flag=-1;
    while (l+1<r){
        mid=(l+r)/2;
        flag=0;
        for(i=0;i<=topl;i++){
            set_range(el[i].d);
            for(j=1;j<=n;j++){
                x=mid-map[el[i].x][j];
                y=mid-map[el[i].y][j];
                if(x<0&&y<0){
                    add_seg(0,el[i].d);
                    break;
                }
                if(x>=el[i].d||y>=el[i].d) continue;

                if(x+y>=el[i].d) continue;
                add_seg(max(0,x+1),min(el[i].d,el[i].d-y-1));

            }
            if(!full())flag=1;
            if(flag==1) break;
        }
        if(flag==0){
            l=mid;
            continue;
        }
        if(flag==1){
            r=mid;
            continue;
        }
    }
    double ans=r/2.0;
    printf("%.2lf",ans);
    return 0;
}
时间： 2024-10-03 06:17:21

20161006模拟的相关文章

CentOS系统启动及内核大破坏模拟实验

讲过了centos的启动流程,此时是不是想来点破坏呢?那就尽情的玩耍吧,记得在实验之前拍个快照,万一哪个环节错误恢复不回来了呢,毕竟数据无价,话不多说,开始. 一.删除伪系统根.(ramdisk文件) (1)模拟误操作删除ramdisk文件. ①模拟误删除initramfs-3.10.0-514.el7.x86_64.img文件. ②为当前正在使用的内核重新制作ramdisk文件 格式为:mkinitrd /boot/initramfs-$(uname -r).img $(uname -r) (

NYOJ 2356: 哈希计划【模拟】

题目描述 众所周知,LLM的算法之所以菜,就是因为成天打游戏,最近LLM突然想玩<金庸群侠传X>,结果进去后各种被虐,LLM就开始研究这个游戏的代码,顺便还学会了一点点点点lua语言,然后就开始了伟大的改游戏代码之旅,然后LLM发现自己too young了,这个游戏把所有的文本都进行了哈希,如果自己改了代码或者剧情文本的话它哈希出来的值就会和原来的哈希值不一样......然后游戏就会打不开.....,现在LLM发现了文本的哈希函数,要求你写个程序,功能为: 输入一段字符串,输出一个哈希值 为了

爬虫——模拟点击动态页面

动态页面的模拟点击: 以斗鱼直播为例:http://www.douyu.com/directory/all 爬取每页的房间名.直播类型.主播名称.在线人数等数据,然后模拟点击下一页,继续爬取 #!/usr/bin/python3 # -*- conding:utf-8 -*- __author__ = 'mayi' """ 动态页面的模拟点击: 模拟点击斗鱼直播:http://www.douyu.com/directory/all 爬取每页房间名.直播类型.主播名称.在线人数

爬虫——网站模拟登录

使用Selenium与PhantomJS模拟登录豆瓣:https://www.douban.com/ #!/usr/bin/python3 # -*- conding:utf-8 -*- __author__ = 'mayi' """ 模拟登录豆瓣:https://www.douban.com/ """ from selenium import webdriver # 调用环境变量指定的PhantomJS浏览器创建浏览器对象,executable

python爬虫 模拟登陆校园网-初级

最近跟同学学习爬虫的时候看到网上有个帖子,好像是山大校园网不稳定,用py做了个模拟登陆很有趣,于是我走上了一条不归路..... 先上一张校园网截图 首先弄清一下模拟登陆的原理: 1:服务器判定浏览器登录使用浏览器标识,需要模拟登陆 2: 需要post账号,密码,以及学校id python走起,我用的2.7版本,用notepad++写的,绑定python可以直接运行 由于是模拟网页登陆,需要导入urllib urllib2 cookielib库,前两个有与网页直接的接口,cookielib就是用来

Android模拟位置信息

Android模拟位置程序,俗称GPS欺骗,只能修改采用GPS定位的软件. 手机定位方式目前有4种:基站定位,WIFI定位,GPS定位,AGPS定位 常见的修改手法: 1. 抓包欺骗法,抓包改包欺骗服务器端, 但是得专门去针对某款app,而且现在很多app数据包都加密了 2. hook java层经纬度获取函数, 这个方法以前可以用,现在不行了 3. hook native层经纬度获取函数 4. 使用允许模拟地址位置信息(不是很通用有版本限制) 为了修改微信朋友圈地理位置信息,为了好玩 试过了上

CSU 1112: 机器人的指令【模拟题】

1112: 机器人的指令 Time Limit: 1 Sec  Memory Limit: 128 MB Submit: 1858  Solved: 682 [Submit][Status][Web Board] Description 数轴原点有一个机器人.该机器人将执行一系列指令,你的任务是预测所有指令执行完毕之后它的位置. ·LEFT:往左移动一个单位 ·RIGHT: 往右移动一个单位 ·SAME AS i: 和第i 条执行相同的动作.输入保证i 是一个正整数,且不超过之前执行指令数 In

Regionals 2015 &gt;&gt; Europe - Central &gt;&gt;7325 - Book Borders【模拟】

Europe - Central >>7325 - Book Borders 题目链接:7325 题目大意:给你一个字符串(含空格),每行x个字符,将单词排列进去,单词不能断开,问每行第一个单词的长度时多少,注意加空格 题目思路:直接模拟.第一个for遍历[a,b],第二个大致为n/a.复杂度大概为nlogn. 开两个数组,v[i]记录i这个位置所属的单词开始位置,e[v[i]]记录第i个位置所属的单词的结束位置. 然后每次判断这一行结尾,所在位置.如果在两个单词中间,则将该单词视为下一行开始

GNS 3模拟防火墙ASA

[模拟环境] 所使用的GNS3版本为0.7.4,如果低于这个版本,有些版本会缺少些选项无法支持. [ASA]     ASA有2种模式的编译文件,分别为单模式和多模式,可选择使用.我使用的是单模式,我试用过多模式,不太好用. [配置] 打开GNS3,编辑→首选项→Qemu→ASA:     添加单模式:Identifier name:asa802-k8-sing(自己填名称,但不能是中文)RAM:256(使用默认的256)Number of NICs:6(网卡数量,默认是6)NIC model: