模板

原创 2013年12月05日 21:30:53
//FFT
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
using namespace std;

const int mxn = (1 << 19) + 100;
const double PI = acos(-1.0);
int N = 1 << 19, id = 0, n;
int cnt[mxn];
pair<double,double> P[mxn], PB[mxn], tmp[mxn];

void fill(int s[], int m, int d)
{
    if (m == N)
        P[d] = make_pair(s[id++],0.);
    else 
    {
        fill(s, m<<1, d);
        fill(s, m<<1, d+m);
    }
}

void fill2(int m, int d)
{
    if (m == N) P[d] = tmp[id++];
    else {
        fill2(m<<1, d);
        fill2(m<<1, d+m);
    }
}

inline void FFT(double oper)
{
    for (int d = 0;(1 << d) < N;++d) {
        int m = (1 << d);
        double p0 = PI / m * oper;
        double sinp0 = sin(p0);
        double cosp0 = cos(p0);
        for (int i = 0;i < N;i += (m << 1)) {
            double sinp = 0;
            double cosp = 1;
            for (int j = 0;j < m;++j) {
                double ta = cosp * P[i+j+m].first - sinp*P[i+j+m].second;
                double tb = cosp * P[i+j+m].second + sinp * P[i+j+m].first;
                P[i+j+m].first = P[i+j].first - ta;
                P[i+j+m].second = P[i+j].second - tb;
                P[i+j].first += ta;
                P[i+j].second += tb;
                double tsinp = sinp;
                sinp = sinp * cosp0 + cosp * sinp0;
                cosp = cosp * cosp0 -tsinp * sinp0;
            }
        }
    }
}

void Mul(int p1[],int p2[])
{
    id = 0;
    fill(p1, 1, 0);
    FFT(1.0);
    for (int i = 0;i < N;++i)
        PB[i] = P[i];
    
	id = 0;
    fill(p2, 1, 0);
    FFT(1.0);

    for (int i = 0;i < N;++i)
    {
        tmp[i].first = P[i].first * PB[i].first - P[i].second * PB[i].second;
        tmp[i].second = P[i].first * PB[i].second + P[i].second * PB[i].first;
    }
    id = 0;
    fill2(1, 0);
    FFT(-1.0);
    for (int i = 0;i < N;++i)
    cnt[i] = P[i].first / N + 0.1;
}

int a[250010],b[250010];

int p1[mxn],p2[mxn],ans[mxn];

void solve(int n,int m)
{
	N=n+m;
	while(N!=(N&-N))	N+=(N&-N);

	for(int i=0;i<N;++i)	p1[i]=p2[i]=0;
	for(int i=0;i<n;++i)	p1[i]=(a[i]==0);
	for(int i=1;i<=m;++i)	p2[i-1]=(b[m-i]==0);
	Mul(p1,p2);
	for(int i=m;i<=n;++i)	ans[i-m]=cnt[i-1];

	for(int i=0;i<N;++i)	p1[i]=p2[i]=0;
	for(int i=0;i<n;++i)	p1[i]=(a[i]==1);
	for(int i=1;i<=m;++i)	p2[i-1]=(b[m-i]==1);
	Mul(p1,p2);
	for(int i=m;i<=n;++i)	ans[i-m]+=cnt[i-1];
}

int nxt[250010];

vector<int>V;

void init(int n,int m)
{
	int i,j;
	for(i=1,nxt[0]=j=-1;i<m;i++)
	{
		while(j+1&&b[i]/2-b[j+1]/2)
			j=nxt[j];
		if(b[i]/2==b[j+1]/2)
			j++;
		nxt[i]=j;
	}
	V.clear();
	for(i=0,j=-1;i<n;i++)
	{
		while(j+1&&a[i]/2-b[j+1]/2)
			j=nxt[j];
		if(a[i]/2==b[j+1]/2)
			j++;
		if(j==m-1)
			V.push_back(i-m+1);
	}
}
int getnum(char *s)
{
	int ans=0;
	for(int i=0;s[i];i++)
		ans=ans*2+s[i]-'0';
	return ans;
}
int main()
{
	int n,m,i,j,k;
	char s[20];

	scanf("%d%d",&n,&m);
	for(i=0;i<n;i++)
		scanf("%s",s),a[i]=getnum(s);
	for(i=0;i<m;i++)
		scanf("%s",s),b[i]=getnum(s);

	init(n,m);

	for(i=0;i<n;i++)
		a[i]%=2;
	for(i=0;i<m;i++)
		b[i]%=2;
	if(V.size()==0)
		puts("No");
	else
	{
		solve(n,m);
		puts("Yes");
		int res=-1,pos=-1;
		for(int i=0;i<V.size();++i)
		{
			if(ans[V[i]]<=res)	continue;
			pos=V[i];
			res=ans[pos];
		}
		printf("%d %d\n",m-res,pos+1);
	}
	return 0;
	//a,b里面存现在的数,0或者1,长度分别位n,m,下标从0开始,vector里面存匹配的开始下表
}
//C( n, m ) % p;
#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;

#define N 100010

long long mod_pow(int a,int n,int p)
{
    long long ret=1;
    long long A=a;
    while(n)
    {
        if (n & 1)
            ret=(ret*A)%p;
        A=(A*A)%p;
        n>>=1;
    }
    return ret;
}

long long factorial[N];

void init(long long p)
{
    factorial[0] = 1;
    for(int i = 1;i <= p;i++)
        factorial[i] = factorial[i-1]*i%p;
    //for(int i = 0;i < p;i++)
        //ni[i] = mod_pow(factorial[i],p-2,p);
}

long long Lucas(long long a,long long k,long long p) //求C(n,m)%p p最大为10^5。a,b可以很大!
{
    long long re = 1;
    while(a && k)
    {
        long long aa = a%p;long long bb = k%p;
        if(aa < bb) return 0; //这个是最后的改动!
        re = re*factorial[aa]*mod_pow(factorial[bb]*factorial[aa-bb]%p,p-2,p)%p;//这儿的求逆不可先处理
        a /= p;
        k /= p;
    }
    return re;
}

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        long long n,m,p;
        cin >> n >> m >> p;
        init(p);
        cout << Lucas(n+m,m,p) << "\n";
    }
    return 0;
}

//插头dp
#include<iostream>
#include<cstdio>
#include<vector>
#include<string>
#include<queue>
#include<cstring>
#include<map>
#include<cmath>
#include<algorithm>
using namespace std;
#define  ll long long
#define inf 0x3fffffff
#define N 13
#define MOD 60007
#define STA 1000010
  
int n,m,S;
struct HashMap
{
    int hd[MOD],nxt[STA],sta[STA],f[STA],size;
    void init()
    {
        size = 0;
        memset(hd,-1,sizeof(hd));
    }
    void add(int s,int a)
    {
        int h = s%MOD;
        for(int i=hd[h];i+1;i=nxt[i])
        {
            if(sta[i]==s)
            {
                f[i] = min(f[i],a);
                return;
            }
        }
        nxt[size] = hd[h];
        sta[size] = s;
        f[size] = a;
        hd[h] = size++;
        return;
    }
}hs[2];
  
int mmap[N][N];
int solve()
{
    S = (1<<((m+1)<<1))-1;
    int cur = 0;
    hs[0].init();
    hs[0].add(0,0);
      
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<m;j++ )
        {
            hs[cur=!cur].init();
            for(int k=0;k<hs[cur^1].size;k++)
            {
                int s = hs[cur^1].sta[k];
                int fk = hs[cur^1].f[k]+1;
                int x = (3<<(j<<1))&s;
                int y = (3<<((j+1)<<1))&s;
                int left = x>>(j<<1);
                int up = y>>((j+1)<<1);
                if(mmap[i][j]==1)
                {
                    if(left==0&&up==0)
                    {
                        hs[cur].add(s,fk-1);
                        if(mmap[i+1][j]&&mmap[i][j+1]) 
                        {
                            hs[cur].add(s|(2<<(j<<1))|(2<<((j+1)<<1)),fk);
                            hs[cur].add(s|(3<<(j<<1))|(3<<((j+1)<<1)),fk);
                        }
                    }
                    else if(left==up&&left)
                    {
                        hs[cur].add(s^x^y,fk);
                    }
                    else if((left==0&&up==2)||(left==2&&up==0))
                    {
                        if(mmap[i+1][j]) hs[cur].add(s^y|(2<<(j<<1)),fk);
                        if(mmap[i][j+1]) hs[cur].add(s^x|(2<<((j+1)<<1)),fk);
  
                    }
                    else if((left==0&&up==3)||(left==3&&up==0))
                    {
                        if(mmap[i+1][j]) hs[cur].add(s^y|(3<<(j<<1)),fk);
                        if(mmap[i][j+1]) hs[cur].add(s^x|(3<<((j+1)<<1)),fk);
  
                    }
                }
                else if(mmap[i][j]==2)
                {
                    if(left==2&&up==0)
                        hs[cur].add(s^x,fk);
                    else if(left==0&&up==2)
                        hs[cur].add(s^y,fk);
                    else if(left==0&&up==0)
                    {
                        if(mmap[i+1][j]==1||mmap[i+1][j]==2) hs[cur].add(s|(2<<(j<<1)),fk);
                        if(mmap[i][j+1]==1||mmap[i][j+1]==2) hs[cur].add(s|(2<<((j+1)<<1)),fk);
                    }
  
                }
                else if(mmap[i][j]==3)
                {
                    if(left==3&&up==0)
                        hs[cur].add(s^x,fk);
                    else if(left==0&&up==3)
                        hs[cur].add(s^y,fk);
                    else if(left==0&&up==0)
                    {
                        if(mmap[i+1][j]==1||mmap[i+1][j]==3) hs[cur].add(s|(3<<(j<<1)),fk);
                        if(mmap[i][j+1]==1||mmap[i][j+1]==3) hs[cur].add(s|(3<<((j+1)<<1)),fk);
                    }
                }
                else
                {
                    if(left==0&&up==0) hs[cur].add(s,fk-1);
                }
  
            }
        }
        hs[cur^1].init();
        for(int i=0;i<hs[cur].size;i++)
            hs[cur^1].add((hs[cur].sta[i]<<2)&S,hs[cur].f[i]);
        cur ^= 1;
          
  
    }
    cur ^= 1;
    for(int i=0;i<hs[cur].size;i++)
        if(hs[cur].sta[i]==0) return hs[cur].f[i]-2;
    return 0;
}
  
int main()
{
      
    int ans;
    while(scanf("%d%d",&n,&m)!=EOF&&(m||n))
    {
          
        memset(mmap,0,sizeof(mmap));
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
            {
                scanf("%d",&mmap[i][j]);
                if(mmap[i][j]<2) mmap[i][j] ^= 1;
            }
  
        ans = solve();
        printf("%d\n",ans);
    }
}




C++模板:函数模板、类模板、模板与继承

C++模板:描述             C++提供一种模板的机制来减少代码重复。比如:对于同一样函数使用不同的数据类型,int,double,char等。C++模板属于“元编程”的范畴。 C++ 模...
  • xiaoding133
  • xiaoding133
  • 2013年09月13日 23:09
  • 37133

函数模板与模板函数及模板类与模板的特化

函数模板( Function templates) 模板(Templates)使得我们可以生成通用的函数,这些函数能够接受任意数据类型的参数,可返回任意类型的值,而不需要对所有可能的数据类型进行函数...
  • yc2zgh1314
  • yc2zgh1314
  • 2016年04月25日 00:07
  • 3426

C++类模板中的模板函数

在C++中类模板中可以嵌套函数模板,但是在使用时候有些注意事项如下: 一、如果类模板中参数为template,在函数模板中,函数的不同模板名字是可以重载的,如template void assign...
  • u010585135
  • u010585135
  • 2015年01月24日 10:05
  • 1180

毕业论文开题报告模板

毕 业 设 计(论文) 开 题 报 告   题   目:   基于SSH的楼盘销售系统                  的设计与实现          学      院:    计算机...
  • qq_36524291
  • qq_36524291
  • 2018年01月15日 14:21
  • 13

小程序开之——发送模板消息

1.下发条件说明 支付 当用户在小程序内完成过支付行为,可允许开发者向用户在7天内推送有限条数的模板消息(1次支付可下发3条,多次支付下发条数独立,互相不影响) 提交表单 当用户在小程...
  • qq_33616529
  • qq_33616529
  • 2018年01月16日 22:42
  • 0

静态成员函数 全局变量 局部变量 标识符 结构体 函数模板

点击打开链接 对象指针的指向原则是:小-->大,父类指向子类。 函数指针刚好相反,子类的函数指针可以指向父类,但是父类的函数指针不能指向子类。 在创建派生类对象时,构造函数的执行顺序:基...
  • yigezhanglu
  • yigezhanglu
  • 2018年01月16日 16:35
  • 11

使用模板设计模式解决问题

最近,在做一个系统,这个系统有部分数据是用户通过CSV文件进行导入的,我并不是很常操作这个,在网上搜索到一个操作CSV的文件类,很简单,具体的操作请看连接 http://blog.csdn.net/...
  • zxc_user
  • zxc_user
  • 2018年01月16日 16:14
  • 18

后缀自动机转后缀树模板

终于找到教程了,现在来写一发 题目:Jzoj4072(bzoj3998)弦论 第一问 其实非常好写,连空间都不需要多开一倍,只需要sa数组rank数组和height数组就好了,dfs一次可以求出 ...
  • JacaJava
  • JacaJava
  • 2018年01月12日 19:38
  • 67

【TINY4412】U-BOOT移植笔记:(2)拷贝模板

【TINY4412】U-BOOT移植笔记:(2)拷贝模板 宿主机 : 虚拟机 Ubuntu 16.04 LTS / X64 目标板[底板]: Tin...
  • sinat_20006769
  • sinat_20006769
  • 2018年01月12日 17:08
  • 9

普通平衡树模板

直接上代码 内容都在代码里 //不正经的splay注释版代码 //bzoj 3224 普通平衡树 /*平衡树的本质其实是二叉搜索树,所以很多操作是基于二叉搜索树的操作。 splay的本...
  • Michael_Byrant
  • Michael_Byrant
  • 2018年01月12日 15:12
  • 27
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:模板
举报原因:
原因补充:

(最多只允许输入30个字)