NOIP2015考前模拟！！！

Problem 1 机器人(robot.cpp/c/pas)
【题目描述】
早苗入手了最新的Gundam模型。最新款自然有着与以往不同的功能，那就是它能够自动行走，厉害吧。
早苗的新模型可以按照输入的命令进行移动，命令包括‘E’、‘S’、‘W’、‘N’四种，分别对应东南西北。执行某个命令时，它会向对应方向移动一个单位。作为新型机器人，它可以执行命令串。对于输入的命令串，每一秒它会按命令行动一次。执行完命令串的最后一个命令后，会自动从头开始循环。在0时刻时机器人位于（0,0）。求T秒后机器人所在位置坐标。
【输入格式】
第1行：一个字符串，表示早苗输入的命令串，保证至少有1个命令
第2行：一个正整数T
【输出格式】
2个整数，表示T秒时，机器人的坐标。
【样例输入】
NSWWNSNEEWN
12
【样例输出】
-1 3
【数据范围】
对于60%的数据 T<=500,000 且命令串长度<=5,000
对于100%的数据 T<=2,000,000,000 且命令串长度<=5,000

【注意】
向东移动，坐标改变改变为(X+1,Y);
向南移动，坐标改变改变为(X,Y-1);
向西移动，坐标改变改变为(X-1,Y);
向北移动，坐标改变改变为(X,Y+1);

---

自己想吧

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define LL long long
#define INF 1000000000
#define eps 1e-10
#define sqr(x) (x)*(x)
#define pa pair<int,int>
#define cyc(i,x,y) for(int i=(x);i<=(y);i++)
#define cy2(i,x,y) for(int i=(x);i>=(y);i--)
using namespace std;
inline 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=10*x+ch-'0';ch=getchar();}
    return x*f;
}
char s[5010];
int t,x,y; 
int main()
{
    freopen("robot.in","r",stdin);
    freopen("robot.out","w",stdout);
    scanf("%s",s);
    int len=strlen(s);
    t=read();
    int ti=t/len,it=t%len;
    if(!ti)
    {
        cyc(i,0,t-1)
        {
            if(s[i]=='W') y++;else
            if(s[i]=='S') y--;else
            if(s[i]=='E') x++;else
            if(s[i]=='W') x--;
        }
        printf("%d %d",x,y);
    }else
    {
        cyc(i,0,len-1)
        {
            if(s[i]=='N') y++;else
            if(s[i]=='S') y--;else
            if(s[i]=='E') x++;else
            if(s[i]=='W') x--;
        }
        x=x*ti,y=y*ti;
        cyc(i,0,it-1)
        {
            if(s[i]=='N') y++;else
            if(s[i]=='S') y--;else
            if(s[i]=='E') x++;else
            if(s[i]=='W') x--;
        }
        printf("%d %d",x,y);
    }

    return 0;
}

Problem 2 数列(seq.cpp/c/pas)
【题目描述】
a[1]=a[2]=a[3]=1
a[x]=a[x-3]+a[x-1] (x>3)
求a数列的第n项对1000000007（10^9+7）取余的值。
【输入格式】
第一行一个整数T，表示询问个数。
以下T行，每行一个正整数n。
【输出格式】
每行输出一个非负整数表示答案。
【样例输入】
3
6
8
10
【样例输出】
4
9
19
【数据范围】
对于30%的数据 n<=100；
对于60%的数据 n<=2*10^7；
对于100%的数据 T<=100，n<=2*10^9；

---

矩阵乘法快速幂！！！
第一个代码用法有点高端高端个P
就是说在dev-cpp4.4上编译错了！
但是我写的是对的

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define LL long long
#define INF 1000000000
#define eps 1e-10
#define sqr(x) (x)*(x)
#define pa pair<int,int>
#define cyc(i,x,y) for(int i=(x);i<=(y);i++)
#define cy2(i,x,y) for(int i=(x);i>=(y);i--)
using namespace std;
inline 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=10*x+ch-'0';ch=getchar();}
    return x*f;
}
#define mod (1000000000+7)
struct matrix{
    LL d[3][3];
};
const matrix dw=(matrix){1,0,0,0,1,0,0,0,1};//单位矩阵 
const matrix zy=(matrix){0,0,1,1,0,0,0,1,1};//转移矩阵
const matrix zero=(matrix){0,0,0,0,0,0,0,0,0};//零矩阵 
matrix mul(matrix A,matrix B)
{
    matrix ans=zero;
    cyc(i,0,2) cyc(j,0,2) cyc(k,0,2)
        ans.d[i][j]=(ans.d[i][j]+A.d[i][k]*B.d[k][j])%mod;
    return ans;
} 
matrix matrix_pow(int x)
{
    matrix r=dw,tmp=zy;
    while(x)
    {
        if(x&1) r=mul(r,tmp);
        tmp=mul(tmp,tmp);
        x>>=1;
    }
    return r;
}
int solve(int x)
{
    if(x<=3) return 1;
    matrix ans=matrix_pow(x-3);
    return (ans.d[2][0]+ans.d[2][1]+ans.d[2][2])%mod;
}
int main()
{
    freopen("seq.in","r",stdin);
    freopen("seq.out","w",stdout);
    int T,x;
    T=read();
    while(T--)
    {
        x=read();
        printf("%d\n",solve(x));
    }
    return 0;
}

第二份代码
改成数组，不用结构体
编译过了

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define LL long long
#define INF 1000000000
#define eps 1e-10
#define sqr(x) (x)*(x)
#define pa pair<int,int>
#define cyc(i,x,y) for(int i=(x);i<=(y);i++)
#define cy2(i,x,y) for(int i=(x);i>=(y);i--)
using namespace std;
inline 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=10*x+ch-'0';ch=getchar();}
    return x*f;
}
#define mod (1000000000+7)
struct matrix{
    LL d[3][3];
};
//const matrix dw=(matrix){1,0,0,0,1,0,0,0,1};//单位矩阵 
//const matrix zy=(matrix){0,0,1,1,0,0,0,1,1};//转移矩阵
//const matrix zero=(matrix){0,0,0,0,0,0,0,0,0};//零矩阵 
matrix mul(matrix A,matrix B)
{
    matrix ans;
    ans.d[0][0]=0;
    ans.d[0][1]=0;
    ans.d[0][2]=0;
    ans.d[1][0]=0;
    ans.d[1][1]=0;
    ans.d[1][2]=0;
    ans.d[2][0]=0;
    ans.d[2][1]=0;
    ans.d[2][2]=0;
    cyc(i,0,2) cyc(j,0,2) cyc(k,0,2)
        ans.d[i][j]=(ans.d[i][j]+A.d[i][k]*B.d[k][j])%mod;
    return ans;
} 
matrix matrix_pow(int x)
{
    matrix r,tmp;
    r.d[0][0]=1;
    r.d[0][1]=0;
    r.d[0][2]=0;
    r.d[1][0]=0;
    r.d[1][1]=1;
    r.d[1][2]=0;
    r.d[2][0]=0;
    r.d[2][1]=0;
    r.d[2][2]=1;
    tmp.d[0][0]=0;
    tmp.d[0][1]=0;
    tmp.d[0][2]=1;
    tmp.d[1][0]=1;
    tmp.d[1][1]=0;
    tmp.d[1][2]=0;
    tmp.d[2][0]=0;
    tmp.d[2][1]=1;
    tmp.d[2][2]=1;
    while(x)
    {
        if(x&1) r=mul(r,tmp);
        tmp=mul(tmp,tmp);
        x>>=1;
    }
    return r;
}
int solve(int x)
{
    if(x<=3) return 1;
    matrix ans=matrix_pow(x-3);
    return (ans.d[2][0]+ans.d[2][1]+ans.d[2][2])%mod;
}
int main()
{
    freopen("seq.in","r",stdin);
    freopen("seq.out","w",stdout);
    int T,x;
    T=read();
    while(T--)
    {
        x=read();
        printf("%d\n",solve(x));
    }
    return 0;
}

Problem 3 虫洞(holes.cpp/c/pas)
【题目描述】
N个虫洞，M条单向跃迁路径。从一个虫洞沿跃迁路径到另一个虫洞需要消耗一定量的燃料和1单位时间。虫洞有白洞和黑洞之分。设一条跃迁路径两端的虫洞质量差为delta。
1.从白洞跃迁到黑洞，消耗的燃料值减少delta，若该条路径消耗的燃料值变为负数的话，取为0。
2.从黑洞跃迁到白洞，消耗的燃料值增加delta。
3.路径两端均为黑洞或白洞，消耗的燃料值不变化。
作为压轴题，自然不会是如此简单的最短路问题，所以每过1单位时间黑洞变为白洞，白洞变为黑洞。在飞行过程中，可以选择在一个虫洞停留1个单位时间，如果当前为白洞，则不消耗燃料，否则消耗s[i]的燃料。现在请你求出从虫洞1到N最少的燃料消耗，保证一定存在1到N的路线。
【输入格式】
第1行：2个正整数N,M
第2行：N个整数，第i个为0表示虫洞i开始时为白洞，1表示黑洞。
第3行：N个整数，第i个数表示虫洞i的质量w[i]。
第4行：N个整数，第i个数表示在虫洞i停留消耗的燃料s[i]。
第5..M+4行：每行3个整数，u,v,k，表示在没有影响的情况下，从虫洞u到虫洞v需要消耗燃料k。
【输出格式】
一个整数，表示最少的燃料消耗。
【样例输入】
4 5
1 0 1 0
10 10 100 10
5 20 15 10
1 2 30
2 3 40
1 3 20
1 4 200
3 4 200
【样例输出】
130
【数据范围】
对于30%的数据： 1<=N<=100,1<=M<=500
对于60%的数据： 1<=N<=1000,1<=M<=5000
对于100%的数据： 1<=N<=5000，1<=M<=30000
其中20%的数据为1<=N<=3000的链
1<=u,v<=N, 1<=k,w[i],s[i]<=200
【样例说明】
按照1->3->4的路线。

---

神奇拆点最短路

#include<cstdio>
#include<cstring>
#define rep(i,n) for(int i=1;i<=n;i++)
#define abs(x) (x > 0 ? x : -(x))
const int maxn=10010,maxm=70010,maxq=16383;
int n,m,u,v,k,sz,l,r,x,delta,
    p[maxn],w[maxn],d[maxn],q[maxq+1],en[maxn],pre[maxm],g[maxm],len[maxm];
bool Inq[maxn];
inline void Ins(int u,int v,int d)
{
    pre[++sz]=en[u];en[u]=sz;g[sz]=v;len[sz]=d;
}
void Init_And_Set_Map()
{
    scanf("%d%d",&n,&m);
    rep(i,n) scanf("%d",&p[i]);
    rep(i,n) scanf("%d",&w[i]);
    rep(i,n)
    {
        scanf("%d",&k);
        Ins(i+i-1,i+i,0);
        Ins(i+i,i+i-1,k);
    }
    rep(i,m)
    {
        scanf("%d%d%d",&u,&v,&k);
        if (p[u]==p[v])
        {
            Ins(u+u-1,v+v,k);
            Ins(u+u,v+v-1,k);
        }
        else
        {
            delta=abs(w[u]-w[v]);
            Ins(u+u-1,v+v-1, (k-delta>0 ? k-delta : 0));
            Ins(u+u,v+v, k+delta);
        }
    }

}
void SPFA()
{
    memset(d,0x7F,sizeof(d));
    q[0]=p[1]+1;d[q[0]]=l=r=0;
    while ((r+1&maxq) != l)
    {
        x=q[l];l=(l+1)&maxq;Inq[x]=0;
        for(int i=en[x];i;i=pre[i])
        {
            v=g[i];
            if (d[x]+len[i]<d[v])
            {
                d[v]=d[x]+len[i];
                if (Inq[v]) continue;
                Inq[v]=1;
                if (d[v]<d[q[l]]) q[l=(l-1)&maxq]=v;
                             else q[r=(r+1)&maxq]=v;
            }
        }
    }

}
int main()
{
    freopen("holes.in","r",stdin);
    freopen("holes.out","w",stdout);
    Init_And_Set_Map();
    SPFA();
    printf("%d",d[n+n-1]<d[n+n] ? d[n+n-1] : d[n+n]);
    return 0;
}