NOIP模拟题[数论}[暴力][拆点]

5 篇文章 0 订阅
3 篇文章 0 订阅

沉迷T3,不能自拔。
T1:
题意:求2的n次方modp,p为质数;
分析:
如果n小,就是快速幂。
如果n大,又是质数,就是费马小定理先把nmod了再直接暴力算。
有人写高精度快速幂,过了70,辛苦了。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<ctime>
#define ll long long 
#define inf 2e8
#define modd 1e9+7
#define clr(x) memset(x,0,sizeof(x))
#define maxen(x) memset(x,127,sizeof(x))
#define maxer(x) memset(x,31,sizeof(x))
#define minus(x) memset(x,-1,sizeof(x))
#define each(i,n) for(int i=1;i<=n;i++)
#define minn(a,b,c) min(a,min(b,c))
#define maxx(a,b,c) max(a,max(b,c))
#ifdef WIN32
#define lld "%I64d"
#else
#define lld "%lld"
#endif
#define PROC "cell"
//for(int i=1;i<=n;i++)
//(double) (ll) LL (int)
//(double)clock()/CLOCKS_PER_SEC
using namespace std;
const int Maxn=3e5+5;
int len,num,p,ans;
char a[Maxn];
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void init()
{
    ans=1;
    scanf("%s",a);
    len=strlen(a);
    scanf("%d",&p);
    if(len>4){
        for(int i=0;i<len;i++)
        num=(num*10+a[i]-'0')%(p-1);}
    else{
        for(int i=0;i<len;i++)
        num=num*10+a[i]-'0';
    }
}
void work()
{
    for(int i=1;i<=num;i++)
        ans=ans*2%p;
    printf("%d",ans);
}
void debug()
{
    //
}
int main()
{
    freopen(PROC".in","r",stdin);
    freopen(PROC".out","w",stdout);
    init();
    work();
    //debug();
    return 0;
}

T2:
题意:
在立方体里装管道,给定出发点和到达点,求最少使用管道数。
分析:
暴力,三维加上一个管道形态,加上一个流出方向。
推增量数组需要耐心。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<ctime>
#define ll long long 
#define inf 2e8
#define modd 1e9+7
#define clr(x) memset(x,0,sizeof(x))
#define maxen(x) memset(x,127,sizeof(x))
#define maxer(x) memset(x,31,sizeof(x))
#define minus(x) memset(x,-1,sizeof(x))
#define each(i,n,m) for(int i=n;i<m;i++)
#define minn(a,b,c) min(a,min(b,c))
#define maxx(a,b,c) max(a,max(b,c))
#ifdef WIN32
#define lld "%I64d"
#else
#define lld "%lld"
#endif
#define PROC "tube"
//for(int i=1;i<=n;i++)
//(double) (ll) LL (int)
//(double)clock()/CLOCKS_PER_SEC
using namespace std;
const int Maxn=25;
int zlx[10][10]={{2,2,3,3,2,2,3,3},
        {0,2,0,1,0,-2,0,-1},
        {2,0,1,0,-2,0,-1,0},
        {-2,-2,-3,-3,-2,-2,-3,-3},
        {0,2,0,1,0,-2,0,-1},
        {2,0,1,0,-2,0,-1,0}};

int zly[10][10]={{0,2,0,1,0,-2,0,-1},
         {2,2,3,3,2,2,3,3},
        {0,2,0,1,0,-2,0,-1},
        {0,2,0,1,0,-2,0,-1},
        {-2,-2,-3,-3,-2,-2,-3,-3},
        {0,2,0,1,0,-2,0,-1},
        };

int zlz[10][10]={{2,0,1,0,-2,0,-1,0},
        {2,0,1,0,-2,0,-1,0},
        {2,2,3,3,2,2,3,3},
        {2,0,1,0,-2,0,-1,0},
        {2,0,1,0,-2,0,-1,0},
            {-2,-2,-3,-3,-2,-2,-3,-3}};

int zls[10][10]={{3,2,3,2,6,5,6,5},
        {3,1,3,1,6,4,6,4},
        {1,2,1,2,4,5,4,5},
        {3,2,3,2,6,5,6,5},
        {3,1,3,1,6,4,6,4},
        {1,2,1,2,4,5,4,5}} ; 

const int lim=Maxn*Maxn*Maxn*Maxn;
int stackx[lim],stacky[lim],stackz[lim],stacks[lim];
int xmax,ymax,zmax,stax,stay,staz,stas,finx,finy,finz,fins;
int vis[Maxn][Maxn][Maxn][Maxn],d[Maxn][Maxn][Maxn][Maxn];
int x,y,z,head,tail;
int qx[lim],qy[lim],qz[lim],qs[lim],xx,yy,zz,ss;
char sym,ano;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void calc(int& s1)
{
    sym=getchar();
    ano=getchar();
    if(ano=='x')s1=1;
    if(ano=='y')s1=2;
    else if(ano=='z')s1=3;
    if(sym=='-')s1+=3;
}
void init()
{
    xmax=read();ymax=read();zmax=read();
    stax=read();stay=read();staz=read();
    calc(stas);
    if(stas==1)stax--;
    else if(stas==2)stay--;
    else if(stas==3)staz--;
    else if(stas==4)stax++;
    else if(stas==5)stay++;
    else if(stas==6)staz++;
    finx=read();finy=read();finz=read();
    calc(fins);maxer(d);
}
void work()
{
    vis[stax][stay][staz][stas]=1;
    d[stax][stay][staz][stas]=0;
    stackx[++tail]=stax;
    stacky[tail]=stay;
    stackz[tail]=staz;
    stacks[tail]=stas;
    head=1;
    while(head<=tail){
        int x=stackx[head];
        int y=stacky[head];
        int z=stackz[head];
        int s=stacks[head++];
        vis[x][y][z][s]=0;
        each(i,0,8){
        int tox=x+zlx[s-1][i];
        int toy=y+zly[s-1][i];
        int toz=z+zlz[s-1][i];
        int tos=zls[s-1][i];
        if(tox<=xmax&&tox>0&&toy<=ymax&&toy>0&&toz<=zmax&&toz>0&&
            d[tox][toy][toz][tos]>d[x][y][z][s]+1){

        d[tox][toy][toz][tos]=d[x][y][z][s]+1;
        if(!vis[tox][toy][toz][tos]){
            vis[tox][toy][toz][tos]=1;
            stackx[++tail]=tox;
            stacky[tail]=toy;
            stackz[tail]=toz;
            stacks[tail]=tos;
        }
        }
    }
    }
    if(d[finx][finy][finz][fins]<=6)
        printf("%d",d[finx][finy][finz][fins]);
    else printf("Impossible");
}
void debug()
{
    //
}
int main()
{
    freopen(PROC".in","r",stdin);
    freopen(PROC".out","w",stdout);
    init();
    work();
    //debug();
    return 0;
}

T3:
题意:
在路上走,转弯前后和出发到达都要时间加倍,求最短时间。
分析:
我当成记忆化搜索来做的,刘汝佳拆点,其实表示方法和我一模一样。就是存下来的方向和是否加倍。
不过有了拆点这个概念以后,我们就可以直接SPFA了,比我搜索要快很多还不容易错(在无向图上用记忆化搜索走单源最短真的很艰难啊!)
易错点就是到了终点可能还要走,但更新答案是不再走的(看程序意会)

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<ctime>
#define ll long long 
#define inf 2e8
#define modd 1e9+7
#define clr(x) memset(x,0,sizeof(x))
#define maxen(x) memset(x,127,sizeof(x))
#define maxer(x) memset(x,31,sizeof(x))
#define minus(x) memset(x,-1,sizeof(x))
#define each(i,n) for(int i=1;i<=n;i++)
#define minn(a,b,c) min(a,min(b,c))
#define maxx(a,b,c) max(a,max(b,c))
#ifdef WIN32
#define lld "%I64d"
#else
#define lld "%lld"
#endif
#define PROC "roller"
//for(int i=1;i<=n;i++)
//(double) (ll) LL (int)
//(double)clock()/CLOCKS_PER_SEC
using namespace std;
const int Maxn=105;
const int zl[4][2]={{0,1},{0,-1},{1,0},{-1,0}};
int tov[Maxn][Maxn][4],f[Maxn][Maxn][5][2];
int n,m,stan,stam,finn,finm,ans=2e8;
queue<int>qx,qy,qd,qp,qf;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void init()
{
    n=read();m=read();
    stan=read();stam=read();
    finn=read();finm=read();
    each(i,2*n-1){
        if(i%2){
            each(j,m-1)
                tov[i/2+1][j][0]=tov[i/2+1][j+1][1]=read();
        }
        else{
            each(j,m)
                tov[i/2][j][2]=tov[i/2+1][j][3]=read();
        }
    }
    maxer(f);
}
void work()
{
    qx.push(stan);
    qy.push(stam);
    qd.push(4);
    qp.push(1);
    qf.push(0);
    f[stan][stam][4][1]=0;
    while(!qx.empty()){
        int x=qx.front();qx.pop();
        int y=qy.front();qy.pop();
        int d=qd.front();qd.pop();
        int p=qp.front();qp.pop();
        int fr=qf.front();qf.pop();
        for(int i=0;i<4;i++)
            if(tov[x][y][i]){
            int tox=x+zl[i][0];
            int toy=y+zl[i][1];
            int pdd;
            if(i==d)pdd=0;
            else pdd=1;
            if(f[tox][toy][i][pdd]>
                    f[x][y][d][p]+(pdd+1)*tov[x][y][i]+((pdd==1&&p==0)?fr:0)){
                f[tox][toy][i][pdd]=
                    f[x][y][d][p]+(pdd+1)*tov[x][y][i]+((pdd==1&&p==0)?fr:0);
                    qx.push(tox);
                    qy.push(toy);
                    qd.push(i);
                    qp.push(pdd);
                    qf.push(tov[x][y][i]);
            }
            if(tox==finn&&toy==finm){
            f[tox][toy][i][pdd]=min(f[tox][toy][i][pdd],
                    f[x][y][d][p]+2*tov[x][y][i]+((p==0&&i!=d)?fr:0));
                ans=min(ans,f[x][y][d][p]+2*tov[x][y][i]+
                ((p==0&&i!=d)?fr:0));
            }
    }
    }
    if(ans==2e8)
    printf("Impossible");
    else printf("%d",ans);
}
void debug()
{
    //
}
int main()
{
    freopen(PROC".in","r",stdin);
    freopen(PROC".out","w",stdout);
    init();
    work();
    //debug();
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值