BestCoder No.14 总结

A题：一个纯物理题，推一下吧。

代码：

#include
#include
#include
using namespace std;
int v,v0;
int main()
{
    while(scanf("%d%d",&v,&v0)!=EOF)
    {
        printf("%d\n",(v<<1)-v0);
    }
    return 0;
}

B题：这个题是说要把地图中所有的石子（值大于0的运到左上角），最开始起点就在左上角。再加上题中给了一次可以运无限个石子并且要求总时间最小，那么这个题就等于求这不到10堆的石子的先后顺序，其实也就是一个TSP问题了。

思路：

      首先由于只有不到10堆的石头，那么很明显容易想到状态压缩。然后如果用DFS的方法只需要枚举这10堆石头的拿的先后顺序即可，不过这里需要一个剪枝。如果用DP的话也比较简单，动态转移方程式dp[i|(1<

代码（DFS）：

#include
#include
#include
#include
#include
using namespace std;
const int maxn=51;
int n,m,cnt,ans;
vector > s;
void DFS(int x,int y,int v,int sum)
{
    if(sum>=ans)
        return;
    if(v==(1<>i)&1))
            DFS(s[i].first,s[i].second,v|(1<

代码（DP）：

#include
#include
#include
#include
#include
using namespace std;
const int inf=1<<29;
const int maxn=11;
int n,m,dp[1< > s;
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        s.clear();
        ans=1<<28;
        for(int i=0;i>j)&1))
                    for(int k=0;k

C题：

      这个题重点在于转成矩阵乘法吧，然后需要用线段树进行维护。相关信息看题解吧，我给出一个不太好的（时间耗费比较高）的版本。

代码：

#include
#include
#include
#include
using namespace std;
const int maxn=5e4+100;
const int mod=1e9+7;
struct Matrix
{
    long long v[2][2];
    Matrix()
    {
        memset(v,0,sizeof(v));
    }
};
struct Node
{
    int l;
    int r;
    Matrix v;
}t[maxn*4];
int n,m;
void Cal(Matrix &ans,Matrix a,Matrix b)
{
    memset(ans.v,0,sizeof(ans.v));
    for(int k=0;k<2;k++)
        for(int i=0;i<2;i++)
            for(int j=0;j<2;j++)
                ans.v[i][j]=(ans.v[i][j]+a.v[i][k]*b.v[k][j]%mod)%mod;
}
void Calc(int index,int s1,int s2)
{
    Cal(t[index].v,t[s1].v,t[s2].v);
}
void Build(int l,int r,int index)
{
    t[index].l=l;
    t[index].r=r;
    if(l==r)
    {
        for(int i=0;i<2;i++)
            for(int j=0;j<2;j++)
                t[index].v.v[i][j]=1;
        return;
    }
    int mid((l+r)/2);
    Build(l,mid,index<<1);
    Build(mid+1,r,index<<1|1);
    Calc(index,index<<1,index<<1|1);
}
void Query(int l,int r,int index,Matrix &ans)
{
    if(t[index].l==l&&t[index].r==r)
    {
        ans=t[index].v;
        return;
    }
    int mid=(t[index].l+t[index].r)>>1;
    if(r<=mid)
        return Query(l,r,index<<1,ans);
    if(l>mid)
        return Query(l,r,index<<1|1,ans);
    Matrix a,b;
    Query(l,mid,index<<1,a);
    Query(mid+1,r,index<<1|1,b);
    Cal(ans,a,b);
    //Cal(ans,Query(l,mid,index<<1,a),Query(mid+1,r,index<<1|1),b);
}
void Update(int l,int x,int y,int index)
{
    if(t[index].l==t[index].r)
    {
        t[index].v.v[x][y]^=1;
        return ;
    }
    int mid((t[index].l+t[index].r)>>1);
    if(l<=mid)
        Update(l,x,y,index<<1);
    else 
        Update(l,x,y,index<<1|1);
    Matrix ans;
    Calc(index,index<<1,index<<1|1);
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        Build(1,n-1,1);
        for(int i=0;i