计蒜之道复赛 B D F

B题是一个简单的模拟

求一下两个点中间每过deta的整数点

然后更新一下每个点的最后一次经过就好了

#include<bits/stdc++.h>
#define cl(a,b) memset(a,b,sizeof(a))
#define debug(x) cerr<<#x<<"=="<<(x)<<endl
using namespace std;

const int maxn=300+10;

int gcd(int a,int b){return a?gcd(b%a,a):b;}
int abs(int a){if(a>0)return a;return -a;}

int a[maxn][maxn];

int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    cl(a,0);
    for(int i=1;i<=n;i++)
    {
        int x1,y1,x2,y2,dx,dy;
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        dx=x2-x1;
        dy=y2-y1;
        int k=gcd(abs(dx),abs(dy));
        dx/=k,dy/=k;
        x2+=dx;
        while(x1!=x2)
        {
//            debug(x1),debug(y1);
            a[x1][y1]=i;
            x1+=dx,y1+=dy;
        }
    }
    int q;
    scanf("%d",&q);
    while(q--)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        printf("%d
",a[x][y]);
    }
    return 0;
}/*

5 8
2 5 5 2
5 2 3 8
8 4 1 4
2 2 5 8
8 7 4 1
4
3 4
5 2
6 4
3 5

*/

D是一个最短路

而且用到了虚点的操作

每个城市群都是一个虚点

也就是对每个城市群需要加两个点 u'和u''

在每个群内的城市连自己的城市群是

u'->v

v->u''

然后每个城市群之间

u'->v''

v'->u''

直接跑一遍spfa就行

需要注意的是：

1.n和m没有规定谁大谁小 所以加点的时候n*2+i和n+i是会冲突的

2.因为最大的值是2e5*2e6>2^31-1 所以ans要用ll存

3.由于上一点的存在 inf也要改……

#include<bits/stdc++.h>
#define cl(a,b) memset(a,b,sizeof(a))
#define debug(x) cerr<<#x<<"=="<<(x)<<endl
using namespace std;
typedef long long ll;

const int MAXN=4e5+10;
const ll INF=2e12;

struct Edge
{
    int v;
    int cost;
    Edge(int _v=0,int _cost=0):v(_v),cost(_cost) {}
};

vector<Edge>E[MAXN];

void addedge(int u,int v,int w)
{
    E[u].push_back(Edge(v,w));
}

bool vis[MAXN];
int cnt[MAXN];
ll dist[MAXN];
int n,m;
int start,over;

ll SPFA(int start,int n)
{
    for(int i=0;i<n+m*2+10;i++)
    {
        dist[i]=INF;
        vis[i]=false;
    }
    vis[start]=true;
    dist[start]=0;
    queue<int>que;
    while(!que.empty())que.pop();
    que.push(start);
    memset(cnt,0,sizeof(cnt));
    cnt[start]=1;
    while(!que.empty())
    {
        int u=que.front();
        que.pop();
        vis[u]=false;
        for(int i=0; i<E[u].size(); i++)
        {
            int v=E[u][i].v;
            if(dist[v]>dist[u]+E[u][i].cost)
            {
                dist[v]=dist[u]+E[u][i].cost;
                if(!vis[v])
                {
                    vis[v]=true;
                    que.push(v);
                    if(++cnt[v]>n) return -1;
                }
            }
        }
    }
    return dist[over];
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int k;
        scanf("%d",&k);
        for(int j=0;j<k;j++)
        {
            int x;
            scanf("%d",&x);
            addedge(x,n+m+i,0);
            addedge(n+i,x,0);
        }
    }
    int m1;
    scanf("%d",&m1);
    for(int i=0;i<m1;i++)
    {
        int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        addedge(u,v,w);
        addedge(v,u,w);
    }
    int m2;
    scanf("%d",&m2);
    for(int i=0;i<m2;i++)
    {
        int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        addedge(n+m+u,n+v,w);
        addedge(n+m+v,n+u,w);
    }
    scanf("%d%d",&start,&over);
    ll ans=SPFA(start,n);
    printf("%lld
",ans>=INF?-1:ans);
    return 0;
}/*

5 4
2 5 1
2 2 4
1 3
2 3 4
2
1 2 9
1 5 18
2
1 2 6
1 3 10
1 5

*/

F是一个状压dp

（虽然没有16这么性感的数字 但是18也很明显了~

第一维表示当前状态 01表示是否删除

第二维表示已经删除的步数

用l r来枚举区间

数据范围不是很大(2^n)*n*n是可以容忍的

（也是调了好久啊……

#include<bits/stdc++.h>
#define cl(a,b) memset(a,b,sizeof(a))
#define debug(x) cerr<<#x<<"=="<<(x)<<endl
using namespace std;
typedef long long ll;

const int maxn=20;
const int mod=1000000007;

int gcd(int a,int b)
{
    return a?gcd(b%a,a):b;
}

bool check(int x,int y)
{
    return x&(1<<y);
}

int n,k;
int a[maxn];
int dp[(1<<maxn)][maxn];

void solve()
{
    for(int i=0; i<(1<<n); i++)
    {
        for(int l=0; l<n; l++)
        {
            int tmp=i;
            if(check(i,l))
            {
                int tmp1=a[l];
                for(int r=l; r<n; r++)
                {
                    if(check(i,r))
                    {
                        tmp^=(1<<r);
                        tmp1=gcd(tmp1,a[r]);
                        if(tmp1<k) break;
                        for(int j=1; j<=n; j++)
                        {
                            dp[i][j]=(dp[i][j]+dp[tmp][j-1])%mod;
                        }
                    }
                }
            }
        }
    }
}

int main()
{
    scanf("%d%d",&n,&k);
    cl(a,0),cl(dp,0);
    for(int i=0; i<n; i++)
    {
        scanf("%d",&a[i]);
    }
    dp[0][0]=1;
    solve();
    ll ans=0;
    for(int i=1; i<=n; i++)
    {
        ans+=1ll*i*dp[(1<<n)-1][i];
        ans%=mod;
    }
    printf("%lld
",ans);
    return 0;
}/*

4 1
1 1 1 1

*/